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

　　实验名称：用追赶法解三对角线方程组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
　　３．分析与讨论