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

void main()
{
    //program d4r13
    //driver for routine betai betacf
    int i;
	char text[20],text1[20];
    double nval,a,b,x,value; 
    const double pi = 3.1415926;
    fstream fin;
    fin.open("d:\\vc常用數值算法集\\data\\fncval.dat",ios::in);
    while ( (strcmp(text,"Incomplete")!=0)||(strcmp(text1,"Beta")!=0) )
	{
      fin>>text;
	  if(strcmp(text,"Incomplete")==0)
		  fin>>text1;
	}
	fin>>text;
    fin>>nval; 
    cout<<"Incomplete Beta Function "<<endl;
    fin>>text;
    cout<<endl;
    cout<<"      a       b       x        actual       betai(x)"<<endl;
    for( i = 1;i<=nval;i++)
	{
        fin>>a;
		fin>>b;
		fin>>x;
		fin>>value;
		cout<<setprecision(3)<<setw(8)<<a;
		cout<<setprecision(3)<<setw(8)<<b;
		cout<<setprecision(3)<<setw(8)<<x;
		cout<<setprecision(6)<<setw(14)<<value;
		cout<<setprecision(6)<<setw(14)<<betai(a, b, x)<<endl;			          
    }
}