CPP   20

xbessi0

Guest on 20th July 2022 04:22:21 PM

  1. #include <string>
  2. #include <fstream>
  3. #include <iostream>
  4. #include <iomanip>
  5. #include "nr.h"
  6. using namespace std;
  7.  
  8. // Driver for routine bessi0
  9.  
  10. int main(void)
  11. {
  12.         string txt;
  13.         int i,nval;
  14.         DP val,x;
  15.         ifstream fp("fncval.dat");
  16.  
  17.         if (fp.fail())
  18.           NR::nrerror("Data file fncval.dat not found");
  19.         getline(fp,txt);
  20.         while (txt.find("Modified Bessel Function I0")) {
  21.           getline(fp,txt);
  22.           if (fp.eof()) NR::nrerror("Data not found in fncval.dat");
  23.         }
  24.         fp >> nval;
  25.         getline(fp,txt);
  26.         cout << endl << "Modified Bessel Function I0" << endl;
  27.         cout << setw(5) << "x" << setw(13) << "actual";
  28.         cout << setw(14) << "bessi0(x)" << endl << endl;
  29.         for (i=0;i < nval;i++) {
  30.           fp >> x >> val;
  31.           cout << fixed << setprecision(2);
  32.           cout << setw(6) << x;
  33.           cout << fixed << setprecision(7);
  34.           cout << setw(13) << val << setw(13) << NR::bessi0(x) << endl;
  35.         }
  36.         fp.close();
  37.         return 0;
  38. }

Raw Paste


Login or Register to edit or fork this paste. It's free.