欢迎来到天天文库
浏览记录
ID:55272727
大小:84.96 KB
页数:3页
时间:2020-05-08
《编写用追赶法解三对角线性方程组的程序,并解下列方程组.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、计算方法与实习上机实验(二)实验名称:编写用追赶法解三对角线性方程组的程序,并解下列方程组:(1)(2)Ax=b,其中A10×10=,b10×1=程序代码:#includeusingnamespacestd;#includeintmain(){floata[100],b[100],c[100],x[100];inti,k,N;while(1){intability=1;//ability用于判断可不可以执行追赶法的操作cout<<"输入三对角矩阵的维度:"<>N;cout<<"
2、输入三对角的数据:"<>b[0]>>c[0]>>x[0];for(i=1;i>a[i]>>b[i]>>c[i]>>x[i];}cin>>a[N-1]>>b[N-1]>>x[N-1];for(k=0;k3、;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<<"此方程的解为:"<4、ut<<"x["<
3、;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<<"此方程的解为:"<4、ut<<"x["<
4、ut<<"x["<
此文档下载收益归作者所有