void setarray(double** A, int n, int m)
{
  int i=0,j=0;
  for(i=0;i<n;i++)
  {
    for(j=0;j<m;j++)
    {
      if(i==j)
        A[0][i*m+j]=1;
      else if(i>j)
        A[0][i*m+j]=0.5;
      else
        A[0][i*m+j]=0.25;
    }
  }
}