t2_4.cpp

本题采用的计算方法为：主要求解三对角阵方程组得解。采用的计算方法为“追赶法”。 算法思路为：求解方程Ly=d（追）——〉求解Ux=y（赶）

//  cpp  t2_4源代码  上机练习题2_4
#include <iostream.h>
#include <math.h>
const int N=20;
int main()
{
	int k;
	double a[N+1],b[N+1],c[N],d[N+1],u[N+1],y[N+1],l[N+1],x[N+1];
	for(k=1;k<=N;k++)         /////写入数据
	{
		if(k>1)a[k]=-1;
		b[k]=4;
		if(k<N)c[k]=-1;
		if(k==1||k==N)
			d[k]=3;
		else
			d[k]=2;
	}
	///////////追赶法求解
	u[1]=b[1];
	y[1]=d[1];
	for(k=2;k<=N;k++)
	{
		l[k]=a[k]/u[k-1];
		u[k]=b[k]-l[k]*c[k-1];
		y[k]=d[k]-l[k]*y[k-1];	
	}
	x[N]=y[N]/u[N];
	for(k=N-1;k>=1;k--)
	{
		x[k]=(y[k]-c[k]*x[k+1])/u[k];	
	}
	/////////////////输出结果
	cout<<"上机练习题2_4"<<endl;
	cout<<"三对角方程组的解为:"<<endl;
	for(k=1;k<=N;k++)
	{
		cout<<"x["<<k<<"]="<<x[k]<<"\t";
		if(k%5==0)cout<<endl;		
	}	
	return 1;
}