Как составить алгоритм???

void main(void)
{double aa,bb;
ipMatr();
for(opMatr(),k=0;k<=n;k++)        //прямой ход метода Гаусса;
 {
     aa=fabs(a[k][k]);
     i=k;
    for(j=k+1;j<=n;j++)//поиск макс. ел-та;
    
    if(fabs(a[j][k])>aa)
    {
        i=j;
        aa=fabs(a[j][k]);
    }
    if(aa<0)
      {
          printf("Система не имеет единственного решения!!!");
       getch();
   
    }
 
  for(j=k;(j<=n+1)&&(i!=k);j++)       //перестановка i-ой строки;
      {bb=a[k][j];
       a[k][j]=a[i][j];
       a[i][j]=bb;
  }
  aa=a[k][k];
  a[k][k]=1;
  for(j=k+1;j<=n+1;j++)//преобразование ведущей строки;
  {
     a[k][j]=a[k][j]/aa;
  }
  for(i=k+1;i<=n;i++)                   //метод Гаусса;
  {
      bb=a[i][k],a[i][k]=0;
    for(j=k+1;(j<=n+1)&&(bb!=0);j++)
    {
        a[i][j]=a[i][j]-bb*a[k][j];
    }
  }
}
opMatr();//Обратный ход метода;
i=n;
while(i>=0)
{   aa=a[i][n+1];
 
    for(j=n;j>=i+1;j--)
    aa=aa-a[i][j]*c0[j];
 
    c0[i]=aa;
    i--;
 
 
}       
 
  i=0;
  while(i<=n)
   {
       printf("c%d=%10.15lf\n",i,c0[i]);
        i++;
   }    
pogrech();
system("PAUSE");
 
}