void main()
{
// Розв'язування системи лінійних рівнянь методом Гауса.
// Нисхідний прорахунок
int i;
int j;
int N;
float r;
float M[20+1][20+1];
//open(unit=1,file='Gaus.txt',status='Old')
//read(1,*)n
//read(1,*)((M(i,j),j=1,N+1),i=1,N)
N=3;
M[1][1]=1;
M[1][2]=1;
M[1][3]=1;
M[1][4]=6;
M[2][1]=1;
M[2][2]=-1;
M[2][3]=1;
M[2][4]=2;
M[3][1]=1;
M[3][2]=1;
M[3][3]=-1;
M[3][4]=0;
l: r=M[1][1];
// do 1 j=1,N+1
// M(1,j)=M(1,j)/r
//1 continue
for(j=1; j<=N+1; j++)
{
M[1][j]=M[1][j]/r;
};
// do 2 i=2,N
// r=M(i,1)
// if (r.EQ.0) goto 9
// do 3 j=1,N+1
// M(i,j)=M(i,j)/r
// M(i,j)=M(i,j)-M(1,j)
//3 continue
//9 continue
//2 continue
for(i=2; i<=N; i++)
{
r=M[i][1];
if(r!=0)
{
for(j=1; j<=N+1; j++)
{
M[i][j]=M[i][j]/r;
M[i][j]=M[i][j]-M[1][j];
};
}
else
{
};
};
// do 4 i=2,N
// do 5 j=1,N+1
// M(i-1,j)=M(i,j)
//5 continue
//4 continue
for(i=2; i<=N; i++)
{
for(j=1; j<=N+1; j++)
{
M[i-1][j]=M[i][j];
};
};
// do 7 i=1,N
// do 8 j=2,N+1
// M(i,j-1)=M(i,j)
//8 continue
//7 continue
for(i=1; i<=N; i++)
{
for(j=2; j<=N+1; j++)
{
M[i][j-1]=M[i][j];
};
};
N--;
if(N>=2)
goto l;
r=M[1][2]/M[1][1];
// r-останній корінь системи
// close(1)
}
