追赶法源程序c语言,追赶法求三对角线方程组程序设计.doc

追赶法求三对角线方程组程序设计.doc

追赶法求三对角线方程组程序设计

实验名称：用追赶法解三对角线方程组C语言程序设计

实验内容：

１．实验目的

掌握追赶法解三对角线方程组求解方法的C语言程序设计思想；

2．实验步骤

1)解决的问题：用追赶法求解三对角线方程组

2)软件： C语言

3)程序：自编

4)公式：

5)程序源代码：

#include

#include

#include

void main()

{

char play;

int i,n;

float a[100],b[100],c[100],d[100];

float u[100],l[100],y[100],x[100];

do

{

textcolor(14);

clrscr();

cprintf("Please input the square distance set of piece n:\r\n");

scanf("%d",n);

clrscr();

printf("Please input a2--an:\n");

for(i=2;i<=n;i++)

{

printf("a%d=",i);

scanf("%f",a[i]);

}

clrscr();

printf("Please input b1--bn:\n");

for(i=1;i<=n;i++)

{

printf("b%d=",i);

scanf("%f",b[i]);

}

clrscr();

printf("Please input c1--c(n-1):\n");

for(i=1;i

{

printf("c%d=",i);

scanf("%f",c[i]);

}

clrscr();

printf("Please input d1--dn:\n");

for(i=1;i<=n;i++)

{

printf("d%d=",i);

scanf("%f",d[i]);

}

clrscr();

u[1]=b[1];

y[1]=d[1];

for(i=2;i<=n;i++)

{

l[i]=a[i]/u[i-1];

u[i]=b[i]-l[i]*c[i-1];

y[i]=d[i]-l[i]*y[i-1];

}

x[n]=y[n]/u[n];

for(i=n-1;i>0;i--)

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

cprintf("The result:\r\n\n");

for(i=n;i>=1;i--)

printf(" x%d=%f\n",i,x[i]);

getche();

printf("\n");

printf("continue?(y/n)");

play=getche();

}while(play=='y'||play=='Y');

sound(555);

delay(10000);

nosound();

}

5)运行过程及结果：

先输入n的值，再分别输入a2---an,b1---bn,c1----1,d1---dn的值即可。

如：用追赶法解三对角方程组

(转载自.ＮＳＥＡＣ.中国评价网)

过程：

Please input the square distance set of piece n:

4↙ (↙表示回车键)

Please input a2--an:

1↙

1↙

2↙

Please input b1--bn:

2↙

3↙

1↙

1↙

Please input c1--c(n-1):

1↙

1↙

1↙

Please input d1--dn:

1↙

2↙

2↙

0↙

结果：

The result:

X4=2.000000

X3=-1.000000

X2=1.000000

X1=0.000000

３．分析与讨论