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

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

实验名称:

编写用追赶法解三对角线性方程组的程序,并解下列方程组: ?2x1?x2?5,??x?2x?x??12,?123(1)?

?x?2x?x?11,34?2???x3?2x4??1(2)Ax=b,其中

??41??-27??1?41??-15????????-15?1?41A10×10=?, b=10×1???

.........???...????15?1?41?????1?4?15????????程序代码：

#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; break; } else { a[k+1]=a[k+1]/b[k]; b[k+1]=b[k+1]-a[k+1]*c[k];

ability=0;

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

//保留一位有效数字

} } } return 0; }

运行结果：