追赶法.cpp

计算方法的追赶法

#include<iostream.h>
#include<math.h>
#include<stdlib.h>
void main()
{
	int N,i;
    cout<<" 输入主对角线维数N：\n";
    cin>>N;
    double *a= new double[N-1];
    double *b=new double[N];
    double *c= new double[N-1];
    double *f=new double[N];
    cout<<"输入数组a[]:\n";
    for(i=0;i<N-1;i++)
        cin>>a[i];
        cout<<"输入数组b[]:\n";
    for(i=0;i<N;i++)
        cin>>b[i];
        cout<<"输入数组c[]:\n";
    for(i=0;i<N-1;i++)
        cin>>c[i];
        cout<<"输入数组f[]:\n";
    for(i=0;i<N;i++)
        cin>>f[i];
    for(i=0;i<N-1;i++)
	{
        c[i]=c[i]/b[i];
        b[i+1]=b[i+1]-a[i]*c[i];
      if(i==0) f[i]=f[i]/b[i];
      else f[i]=(f[i]-a[i-1]*f[i-1])/b[i];
	}
    f[N-1]=(f[N-1]-a[N-2]*f[N-2])/b[N-1];
    for(i=N-2;i>=0;i--)
       f[i]=f[i]-c[i]*f[i+1];
    for(i=0;i<N;i++)
       cout<<"计算结果为："<<f[i]<<"　,";
    exit(0);
}