多项式.cpp

在数值分析中多项式计算是一种比较复杂的运算

#include<iostream.h>
#include<math.h>
#include<iomanip.h>
void main()
{  
   int i=0;
   int k=0;
   int n=0;
   float x[4]={1,1,2,3};
   float y[4]={1,14,23,33};
   float a[4]={1,1,2,3};
   float u[4]={1,0,2,3};
   float v[4]={1,1,2,3};
   float t[4]={1,1,2,3};
   float sy=0;
   float u1=0;
   float s1=0;
   float f=0;
   float q=0;
   float p=0;
   float s0=0;
   cin>>n;
   cout<<n<<endl;
 /*for(i=1;i<=4;i++)
   {
	   cin>>x[i];
   }

   for(i=1;i<=4;i++)
   {
	   cin>>y[i];
   } */

   for(i=1;i<=4;i++)
   {
	   sy+=y[i];
	   u1+=x[i];
   }
   a[0]=sy/4;
   q=1;
   u[0]=u1/4;
   v[0]=0;
   s0=4;
   for(k=0;k<=n-1;k++)
   {
	   s1=0;
	   sy=0;
	   for(i=1;i<=4;i++)
	   {
		   f=(x[i]-u[k])*q-v[0]*p;
		   sy+=f*y[i];
           s1+=f*f;
           u[k+1]+=x[i]*f*f;
		   p=q;
		   q=f; 
		   
	   }
	   a[k+1]=(sy)/(s1);
	   u[k+1]=u[k+1]/(s1);
	   v[k+1]=(s1)/(s0);
	   (s0)=(s1);
	   cout<<a[k+1]<<" ";
   }
    cout<<endl;
    
   for(i=1;i<=4;i++)
   {
      cin>>t[i];

      s0=1;
	  s1=t[i]-u[0];
	  f=a[0]*(s0)+a[1]*(s1);
	  for(k=1;k<=n-1;k++)
	  {
          sy=(t[i]-u[k])*(s1)-v[k]*(s0);
		  f=f+a[k+1]*(sy);
		  s0=s1;
		  s1=f;

	  }
	  cout<<"f="<<f<<endl;


   }
}