FISH.


Nikita Beloshenko

Create Your Badge



Shop (under construction).





System of the linear equations.

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)
}




Click Here for Free Traffic!

FISH.