数值计算方法与算法:C语言实现利用Gauss消元法解方程组

源程序:

#include<stdio.h>

#include<math.h>

#define MAX_N 20

void main()

{  double a[MAX_N][MAX_N];

   double b[MAX_N];

   double x[MAX_N];

   double t;

   double s;

   int n;

   int i;

   int j;

   int k;

   printf("Please input the dimension of the matrix;");//输入方程组AX=B的维数

   scanf("%d",&n);

  

   printf("Now input the coefficient matrix a(i,j),i,j=1,...,%d:/n",n);//输入系数A矩阵

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

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

                 scanf("%lf",&a[i][j]);

          printf("Now input the matrix b(i),i=1,...,%d:/n",n);//输入b矩阵

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

                 scanf("%lf",&b[i]);

          for(k=1;k<=n-1;k++)//消元过程

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

          {t=a[i][k]/a[k][k];

          for(j=k+1;j<=n;j++)

                 {a[i][j]=a[i][j]-t*a[k][j];}

          b[i]=b[i]-t*b[k];

          }

          }

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

          {s=b[i];

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

                 s=s-a[i][j]*x[j];

                 x[i]=s/a[i][i];

          }

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

                 printf("x[%d]=%lf/n",i,x[i]);

                

输入提示:

在输入系数矩阵时,先输入维数

输入矩阵,一行一行地输入

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值