計算方法與實習(xí)上機實驗(二)

實驗名稱:

編寫用追趕法解三對角線性方程組的程序,并解下列方程組:

(1)

(2)Ax=b,其中

A10×10=, b10×1=

程序代碼：

#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用于判斷可不可以執(zhí)行追趕法的操作

cout<

cin>>N;

cout<

cin>>b[0]>>c[0]>>x[0];

for(i=1;i

{

cin>>a[i]>>b[i]>>c[i]>>x[i];

}

cin>>a[N-1]>>b[N-1]>>x[N-1];

for(k=0;k

{

if(b[k]==0)

{

cout<

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];//這個過程執(zhí)行的是消元過程(即追趕法的追):對應(yīng)于書上的β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--) //下標(biāo)從大到小變化,是趕的過程

{

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

}

cout<

for(i=0;i

{

cout<

cout<

}

}

}

return 0;

}

運行結(jié)果：

總結(jié)

以上是生活随笔為你收集整理的c语言程序 用追赶法求解方程组,编写用追赶法解三对角线性方程组的程序,并解下列方程组（3页）-原创力文档...的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。