#include <iostream.h>
#include <math.h>
#include <iomanip.h>
#include <stdlib.h>
#include <fstream.h>
#include <string>
#include <process.h>

void main()
{
    //program d4r26
    //driver for routine besjan
    int i,n[19],ih;
	char text[20];
    double nval,a[19],value[19],x[19],f[19],temp; 
    const double pi = 3.1415926;
	temp=4;
    fstream fin;
    fin.open("d:\\vc常用數值算法集\\data\\fncval.dat",ios::in);
    while ( strcmp(text,"Ja+n")!=0) 
	{
      fin>>text;
	}
    fin>>nval; 
    fin>>text;
	cout<<"Bessel Function Ja+n"<<endl;
    cout<<endl;
    cout<<"  n    a    ih    x     actual     besjan"<<endl;
    for (i = 0 ; i<=nval - 1; i++)
	{
        fin>>n[i];
        fin>>a[i];
        fin>>x[i];
        fin>>value[i];
    }
    ih = 1;
    besjan(x[0], a[0], temp, ih, f);
    for( i = 0;i<=4;i++)
	{
		cout<<setw(6)<<n[i];
        cout<<setw(6)<<a[i];
		cout<<setw(6)<<ih;
		cout<<setw(6)<<x[i];
		cout<<setw(16)<<value[i];
		cout<<setw(16)<<f[i]<<endl;
	}
    besjan(x[5], a[5], temp, ih, f);
    for( i = 5;i<=9;i++)
	{
		cout<<setw(6)<<n[i];
        cout<<setw(6)<<a[i];
		cout<<setw(6)<<ih;
		cout<<setw(6)<<x[i];
		cout<<setw(16)<<value[i];
		cout<<setw(16)<<f[i-5]<<endl;
    }
    ih = -1;
    besjan(x[10], a[10], temp, ih, f);
    for( i = 10;i<=13;i++)
	{
		cout<<setw(6)<<n[i];
        cout<<setw(6)<<a[i];
		cout<<setw(6)<<ih;
		cout<<setw(6)<<x[i];
		cout<<setw(16)<<value[i];
		cout<<setw(16)<<f[i-9]<<endl;
    }
    besjan(x[14], a[14], temp, ih, f);
    for( i = 14;i<=17;i++)
	{
		cout<<setw(6)<<n[i];
        cout<<setw(6)<<a[i];
		cout<<setw(6)<<ih;
		cout<<setw(6)<<x[i];
		cout<<setw(16)<<value[i];
		cout<<setw(16)<<f[i-13]<<endl;
    }
}