C语言 用Gauss列主元消元法解线性方程组

维他C

1. 
   
2. #include<math.h>
   
3. #include<stdio.h>
   
4. #define NUMBER 20
   
5. #define Esc   0x1b
   
6. #define Enter 0x0d
   
7. float A[NUMBER][NUMBER+1] ,ark;
   
8. int flag,n;
   
9. exchange(int r,int k);
   
10. float max(int k);
    
11. message();
    
12. main()
    
13. {
    
14.    float x[NUMBER];      /*此数组用于存放方程解*/
    
15.    int r,k,i,j;
    
16.    char celect;
    
17.    clrscr();
    
18.    
    
19.    printf("\n\n用Gauss列主元消元法解线性方程组");
    
20.    printf("\n\n1.解方程组请按Enter.");
    
21.    printf("\n\n2.退出程式请按Esc.");
    
22.    celect=getch();
    
23.    if(celect==Esc)
    
24.      exit(0);
    
25.    printf("\n\n 输入方程组的维数：n=");
    
26.    scanf("%d",&n);
    
27.      printf(" \n\n现在输入系数矩阵A和向量b:");
    
28.    for(i=1;i<=n;i++)
    
29.    {
    
30.     printf("\n\n请输入a%d1--a%d%d系数和向量b%d:",i,i,n,i);
    
31.          /*实现将每一行中的系数和向量一次性输入，数之间用空格格开，输完后回车确定*/
    
32.     for(j=1;j<=n+1;j++)     /*将刚才输入的数存入数组*/
    
33.     scanf("%f",&A[i][j]);
    
34.    }
    
35.    for(k=1;k<=n-1;k++)                       
    
36.    {
    
37.     ark=max(k);
    
38.     if(ark==0)                /*判断方程是否为线性方程，即是否合法*/
    
39.     {
    
40.       printf("\n\n此方程组不合法!");message();   
    
41.     }
    
42.     else if(flag!=k)
    
43.      exchange(flag,k);
    
44.      for(i=k+1;i<=n;i++)
    
45.      for(j=k+1;j<=n+1;j++)
    
46.      A[i][j]=A[i][j]-A[k][j]*A[i][k]/A[k][k];
    
47.    }
    
48.    x[n]=A[n][n+1]/A[n][n];
    
49.    for( k=n-1;k>=1;k--)
    
50.    {
    
51.      float me=0;
    
52.      for(j=k+1;j<=n;j++)
    
53.      {
    
54.        me=me+A[k][j]*x[j];
    
55.      }
    
56.        x[k]=(A[k][n+1]-me)/A[k][k];
    
57.    }
    
58.    for(i=1;i<=n;i++)
    
59.    {
    
60.      printf(" \n\nx%d=%f",i,x[i]);
    
61.    }
    
62.    message();
    
63. }
    
64. exchange(int r,int k)         /*交换行的矩函数*/
    
65. {
    
66.   int i;
    
67.   for(i=1;i<=n+1;i++)
    
68.     A[0][i]=A[r][i];
    
69.   for(i=1;i<=n+1;i++)
    
70.     A[r][i]=A[k][i];
    
71.   for(i=1;i<=n+1;i++)
    
72.     A[k][i]=A[0][i];
    
73. }
    
74. float max(int k)           /*比校系数大小的函数*/
    
75. {
    
76.   int i;
    
77.   float temp=0;
    
78.   for(i=k;i<=n;i++)
    
79.     if(fabs(A[i][k])>temp)
    
80.     {
    
81.       temp=fabs(A[i][k]);
    
82.       flag=i;
    
83.     }
    
84.   return temp;
    
85. }
    
86. message()                                       /*实现菜单选择的函数*/
    
87. {
    
88.   printf("\n\n 继续运算按 Enter ,退出程式按 Esc!");
    
89.   switch(getch())
    
90.   {
    
91.    case Enter: main();
    
92.    case Esc: exit(0);
    
93.    default:{printf("\n\n不合法的输入!");message();}
    
94.   }
    
95. }

复制代码

CCC

看看先~~~

wudijiangjuna

有错误啊