//Gauss 消去法解线性方程组
//参考教材《计算方法教程》第二版,西安交通大学出版社
#include<stdio.h>
int main(void)
{
float A[7][7]={{3,-5,6,4,-2,-3,8},
{1,1,-9,15,1,-9 ,2},
{2,-1,7,5,-1,6,11},
{-1,1,3,2,7,-1,-2},
{4,3,1,-7,2,1,1},
{2,9,-8,11,-1,-4,-1},
{7,2,-1, 2,7,-1,9}};
float b[7]={11,2,29,9,5,8,25};
float x[7]={0};
float Aik,S;
int i,j,k;
int size=7;
printf("A[][]");
for(i=0;i<size;i++)
{
for(j=0;j<size;j++)
printf("%f ",A[i][j]);
printf("");
}
printf("b[]");
for(i=0;i<size;i++)
printf("%f ",b[i]);
printf("");
//消去过程
for(k=0;k<size-1;k++)
{
if(!A[k][k])
return -1;
for(i=k+1;i<size;i++)
{
Aik=A[i][k]/A[k][k];
for(j=k;j<size;j++)
{
A[i][j]=A[i][j]-Aik*A[k][j];
}
b[i]=b[i]-Aik*b[k];
}