# include # include # include # include using namespace std; # include "besselj_zero.h" int main ( ); void timestamp ( ); //****************************************************************************80 int main ( ) //****************************************************************************80 // // Purpose: // // besselj_zero_test() tests besselj_zero(). // // Licensing: // // This code is distributed under the MIT license. // // Modified: // // 13 June 2025 // // Author: // // John Burkardt // { int i; int n; int nt = 10; double *z0; double *z1; timestamp ( ); cout << "\n"; cout << "besselj_zero_test():\n"; cout << " C++ version\n"; cout << " Test besselj_zero()\n"; n = 0; z0 = jn_zeros ( n, nt ); n = 1; z1 = jn_zeros ( n, nt ); cout << "\n"; cout << " i J0(root i) J1(root i)\n"; cout << "\n"; for ( i = 0; i < nt; i++ ) { cout << " " << setw(2) << i << " " << setw(14) << z0[i] << " " << setw(14) << z1[i] << "\n"; } delete [] z0; delete [] z1; // // Terminate. // cout << "\n"; cout << "besselj_zero_test():\n"; cout << " Normal end of execution.\n"; timestamp ( ); return 0; } //****************************************************************************80 void timestamp ( ) //****************************************************************************80 // // Purpose: // // timestamp() prints the current YMDHMS date as a time stamp. // // Example: // // 31 May 2001 09:45:54 AM // // Licensing: // // This code is distributed under the MIT license. // // Modified: // // 19 March 2018 // // Author: // // John Burkardt // { # define TIME_SIZE 40 static char time_buffer[TIME_SIZE]; const struct std::tm *tm_ptr; std::time_t now; now = std::time ( NULL ); tm_ptr = std::localtime ( &now ); std::strftime ( time_buffer, TIME_SIZE, "%d %B %Y %I:%M:%S %p", tm_ptr ); std::cout << time_buffer << "\n"; return; # undef TIME_SIZE }