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

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