#include<stdio.h>
#include<math.h>
#include<conio.h>
#define N 3
void Guass(double a[][N+1],double x[])
{
double l[N];
int i,j,k;
l[0]=1.0;
for(i=0;i<N-1;i++)
{
if(fabs(a[i][i])<0.000001)
{
printf("不能求解!/n");
return;
}
for(j=i+1;j<N;j++)
l[j]=a[j][i]/a[i][i];
for(j=i+1;j<N;j++)
for(k=i+1;k<N+1;k++)
a[j][k]=a[j][k]-l[j]*a[i][k];
}
#include<math.h>
#include<conio.h>
#define N 3
void Guass(double a[][N+1],double x[])
{
double l[N];
int i,j,k;
l[0]=1.0;
for(i=0;i<N-1;i++)
{
if(fabs(a[i][i])<0.000001)
{
printf("不能求解!/n");
return;
}
for(j=i+1;j<N;j++)
l[j]=a[j][i]/a[i][i];
for(j=i+1;j<N;j++)
for(k=i+1;k<N+1;k++)
a[j][k]=a[j][k]-l[j]*a[i][k];
}
for(i=N-1;i>=0;i--)
{
x[i]=a[i][N];
for(j=i+1;j<N;j++)x[i]=x[i]-a[i][j]*x[j];
x[i]=x[i]/a[i][i];
}
}
void main()
{
int i,j;
double x[N];
double s[N][N+1]=
{
x[i]=a[i][N];
for(j=i+1;j<N;j++)x[i]=x[i]-a[i][j]*x[j];
x[i]=x[i]/a[i][i];
}
}
void main()
{
int i,j;
double x[N];
double s[N][N+1]=