《编写用追赶法解三对角线性方程组的程序,并解下列方程组.docx》由会员分享,可在线阅读,更多相关《编写用追赶法解三对角线性方程组的程序,并解下列方程组.docx(5页珍藏版)》请在三一文库上搜索。
1、-15 .计算方法与实习上机实验(二) 实验名称:编写用追赶法解三对角线性方程组的程序,并解下列方程组:(1)2 x -x =5,1 2-x +2 x -x =-12, 1 2 3-x +2 x -x =11, 2 3 4-x +2 x =-13 4(2)Ax=b,其中-4 1 -27A1010=1 -4 11 -4 1. . . -15 , b =101 程序代码:#includeusing namespace std;#includeint main()1 -4 11 -4 -15 -15float a100,b100,c100,x100;int i,k,N;while(1)int abi
2、lity=1;break;else ak+1=ak+1/bk;ability=0;bk+1=bk+1-ak+1*ck;xk+1=xk+1-ak+1*xk;/这个过程执行的是消元过程(即追赶法的追):对应于书上的i=bi-lic(i-1),yi=di-liy(i-1)if(ability)xN-1=xN-1/bN-1; /回代法的第一项for(i=N-2;i=0;i-) / 下标从大到小变化 ,是 赶的过程xi=(xi-ci*xi+1)/bi;cout此方程的解为:endl;for(i=0;iN;i+)coutsetiosflags(ios:showpoint);coutxi+1=setiosflags(ios:fixed)setprecision(1)xiendl; /保留一位有效数字return 0;运行结果: