c语言程序 用追赶法求解方程组,编写用追赶法解三对角线性方程组的程序,并解下列方程组(3页)-原创力文档...
計(jì)算方法與實(shí)習(xí)上機(jī)實(shí)驗(yàn)(二)
實(shí)驗(yàn)名稱:
編寫用追趕法解三對(duì)角線性方程組的程序,并解下列方程組:
(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];//這個(gè)過程執(zhí)行的是消元過程(即追趕法的追):對(duì)應(yīng)于書上的βi=bi-lic(i-1),yi=di-liy(i-1)
}
}
if(ability)
{
x[N-1]=x[N-1]/b[N-1]; //回代法的第一項(xiàng)
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;
}
運(yùn)行結(jié)果:
總結(jié)
以上是生活随笔為你收集整理的c语言程序 用追赶法求解方程组,编写用追赶法解三对角线性方程组的程序,并解下列方程组(3页)-原创力文档...的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 猴子选王c语言链表程序代码,c/c++编
- 下一篇: linux as86,记linux_ce