chase.cpp

追赶法,很好 追赶法,很好 追赶法,很好

//追赶法
#include <iostream.h>

#define n 4
void main()
{
	int i;
	double a[n]={0,2,1,2};
	double b[n]={10,2,10,1};
	double c[n]={5,1,5,0};
	double f[n]={5,3,27,6};
	double u[n],y[n],x[n];

	u[0]=c[0]/b[0];
    y[0]=f[0]/b[0];
	for(i=1;i<n;i++){
		u[i]=c[i]/(b[i]-a[i]*u[i-1]);
		y[i]=(f[i]-a[i]*y[i-1])/(b[i]-a[i]*u[i-1]);
	}
    x[n-1]=y[n-1];
	for(i=n-2;i>=0;i--)x[i]=y[i]-u[i]*x[i+1];

	for(i=0;i<n;i++)cout<<x[i]<<endl;
}