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

void main()
{
    //program d4r25
    //driver for routine bessi
    int i,n;
	char text[20];
    double nval,value,x; 
    const double pi = 3.1415926;
    fstream fin;
    fin.open("d:\\vc常用數值算法集\\data\\fncval.dat",ios::in);
    while ( strcmp(text,"In,")!=0 ) 
	{
      fin>>text;
	}
	fin>>text;
    fin>>nval; 
    fin>>text;
	cout<<"Modified Bessel Function In "<<endl;
    cout<<endl;
    cout<<"      n      x          actual          bessi(n,x)"<<endl;
    for( i = 1; i<=nval; i++)
	{
		fin>>n;
        fin>>x;
		fin>>value;  
		cout<<setprecision(2)<<setw(6)<<n;
		cout<<setprecision(2)<<setw(7)<<x;
        cout<<setprecision(6)<<setw(18)<<value;
		cout<<setprecision(6)<<setw(18)<<bessi(n,x)<<endl;			          
    }
}