#include <iostream>


int solve(double** A, double *b, int n)
{
   int i,k;
   
   if (n<=0)
      return 0;

   for (i=0; i<n; i++)
   {
      for(k=0; k<i; k++)
         b[i] -= b[k]*A[i][k];

      if(A[i][i]==0)
         return -i;

      b[i] = b[i]/A[i][i];
   }

   return 1;
}

int main()
{
   double** A = new double*[5];
   double*  b = new double[5];

   for (int i=0; i<5; i++)
      A[i] = new double[i+1];

   for (int i=0; i<5; i++)
      for (int k=0; k<=i; k++)
         A[i][k]=i;

   b[0]=1;
   b[1]=3;
   b[2]=6;
   b[3]=10;
   b[4]=15;

   solve(A,b,5);
   
   for (int i=0; i<5; i++)
      std::cout << "b[" << i << "]=" << b[i] << std::endl;

   for (int i=0; i<5; i++)
      delete A[i];
   delete A;
   delete b;
}