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

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

实验内容：

１．实验目的

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

2．实验步骤

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

2）软件： C语言

3）程序来源：自编

4）公式：

5）程序源代码：

#include <stdio.h>

#include <conio.h>

#include <dos.h>

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

(1

4);

clrscr();

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

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<n;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:rnn");

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

printf(" x%d=%fn",i,x[i]);

getche();

printf("n");

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

play=getche();

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

sound(5

5

5);

delay(10000);

nosound();

}

5）运行过程及结果：

先输入n的值，再分别输入a2---an,b1---bn,c1---cn-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

３．分析与讨论