编写用追赶法解三对角线性方程组的程序,并解下列方程组

计算方法与实习上机实验(二)

实验名称:

编写用追赶法解三对角线性方程组的程序,并解下列方程组:

(1)⎪⎪⎩⎪⎪⎨⎧-=+-=-+--=-+-=-1

2,112,122,5243432

32121x x x x x x x x x x (2)Ax=b,其中

A 10×10=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----41141.........14114114, b 10×1=⎥⎥⎥⎥⎥⎥⎥⎥⎦

⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--1515...15-15-27- 程序代码：

#include

using namespace std;

#include

int main()

{

float a[100],b[100],c[100],x[100];

int i,k,N;

while(1)

{

int ability=1;

ability=0;

break;

}

else

{ a[k+1]=a[k+1]/b[k];

b[k+1]=b[k+1]-a[k+1]*c[k]; x[k+1]=x[k+1]-a[k+1]*x[k];//这个过程执行的是消元过程(即追赶法的追):对应于书上的βi=bi-lic(i-1),yi=di-liy(i-1)

}

}

if(ability)

{

x[N-1]=x[N-1]/b[N-1]; //回代法的第一项

for(i=N-2;i>=0;i--) //下标从大到小变化,是赶的过程

{

x[i]=(x[i]-c[i]*x[i+1])/b[i];

}

cout<<"此方程的解为："<

for(i=0;i

{

cout<

cout<<"x["<

}

}

}

return 0;

}

运行结果：