#include "Int.h"

Integrate::Integrate()
{
     ifstream infile("input.dat");
     if(!infile){
             cout<<"sorry!The input data can not be loaded!"<<endl;
             exit(-1);
	 }

     else{int i,j;
          for(i=0;i<n;i++)
             for(j=0;j<=i;j++)
                 infile >> Xknot[i][j];//輸入求積節(jié)點(diǎn)

          for (i=0;i<n;i++)
             for(j=0;j<=i;j++)
                 infile >> Aknot[i][j];//輸入求積系數(shù)
                 infile >> a >>b;    //輸入?yún)^(qū)間端點(diǎn)值
	 }
     infile.close();
}


void Integrate::Inte()
      {int i;
       double temp,result=0.0;//置result初值為零
       
      cout<<"Please input m（1<=m<=5）"<<endl<<"m=";
      cin>>m;
      m=m-1;//數(shù)組是從零開(kāi)始

      if(a==-1&&b==1)
	     {for(i=0;i<=m;i++)
	       {temp=Function(Xknot[m][i]);//若a=-1且b=1，直接使用式（2.2）

             result+=Aknot[m][i]*temp;
	       }
	      }

      else{
           for(i=0;i<=m;i++)
		      {
               temp=Function(0.5*(b-a)*Xknot[m][i]+0.5*(b+a));//否則，使用（2.3）式計(jì)算
               result+=Aknot[m][i]*temp;
		       }
               result*=0.5*(b-a);//結(jié)果再乘上(a+b)/2
	       }
       cout<<"The result is"<<result<<endl;

}

double Integrate::Function(double x)//要計(jì)算的函數(shù)sin(x^2)
        { double y;
           y=sin(x*x);
           return y;
        }


int main()
    { char cmd;
      cmd='n';

       while(cmd!='y')//當(dāng)從鍵盤(pán)讀入command為yes時(shí)，跳出。
	   {Integrate inter;
             inter.Inte();

             cout<<"Do you want to end?(y/n) "<<endl;
             cin >>cmd;
	   }

        return 1;
     }