java to c++

#include "StdAfx.h"
#include "RealMatrix.h"
 
 
RealMatrix::RealMatrix(void)
{
    col=1;
    row=1;
    matrix=new double*[col];
    for(int i=0;i<col;i++)
    {
        matrix[i]=new double[row];
        for(int j=0;j<row;j++)
            matrix[i][j]=0;
    }
}
 
RealMatrix::RealMatrix(int c,int r,double** mass)
{
    col=c;
    row=r;
    matrix=new double*[row];
    for(int i=0;i<row;i++)
    {
        matrix[i]=new double[col];
        for(int j=0;j<col;j++)
        {
            matrix[i][j]=mass[i][j];
        }
        
    }
        ///for(int i=0;i<row;i++)
    //{
        //for(int j=0;j<col;j++)
            //printf("%d \t",matrix[i][j]);
        //printf("\n"); printf("\n");
        //}
}
 
RealMatrix::~RealMatrix(void)
{
}
 
bool RealMatrix::MatrixMultiply(RealMatrix in)
{
    if(row=in.col)
    {
    
 
double** newMatrix=new double*[row];
    for(int i=0;i<row;i++)
        newMatrix[i]=new double[row];
    
      
      for (int i = 0; i < row; i++) { 
          for (int j = 0; j < in.col; j++) {
            newMatrix[i][j]=0;
           for (int inner = 0; inner < col; inner++) {
                newMatrix[i][j] += matrix[i][inner] * in.matrix[inner][j];
                    //printf("%d * %d + ", matrix[i][inner] , in.matrix[inner][j]);
            }
            //printf("%d ",newMatrix[i][j]);
        }
       // printf("\n");
      }
      matrix=newMatrix;
     return true;
    }
    return false;
}
void RealMatrix::printMatrix()
{
      for (int i = 0; i < row; i++) { 
          for (int j = 0; j <col; j++) {
            printf("%f ",matrix[i][j]);
        }
      printf("\n");
      }
}
int RealMatrix::getCol()
{
    return col;
};
int RealMatrix::getRow()
{
    return row;
}
double** RealMatrix::getMatrix()
{
    return matrix;
}
 
RealMatrix RealMatrix::LUDecomposition(RealMatrix in)
{
    int n=in.getCol();
  double sum = 0;
 
  
  double** A=in.getMatrix();
  double** L ;
  double** U ;
 
    L=new double*[n];
    for(int i=0;i<n;i++)
    {
        L[i]=new double[n];
        for(int j=0;j<n;j++)
            L[i][j]=0;
    }   
 
     U=new double*[n];
    for(int i=0;i<n;i++)
    {
         U[i]=new double[n];
        for(int j=0;j<n;j++)
             U[i][j]=0;
    }
 
//находим первый столбец L[][] и первую строку U[][]
 
  for (int i = 0; i < n; i++)
  {
    L [i][0] = A [i][0];
    U [0][i] = A [0][i] / L [0][0];
  }
 
//дальше вычисляем L[][], U[][] по формуле
 
  for (int i = 1; i < n; i++)
  {
      for (int j = 1; j < n; j++)
      {
          if (i >= j) //нижний треугольник
          {
              sum = 0;
              for (int k = 0; k < j; k++)
                  sum += L [i][k] * U [k][j];
 
              L [i][j] = A [i][j] - sum;
          }
          else // верхний
          {
              sum = 0;
              for (int k = 0; k < i; k++)
                  sum += L [i][k] * U [k][j];
 
              U [i][j] = (A [i][j] - sum) / L [i][i];
          }
      }
   }
 
//====================================================
   printf("\n\n");
   for (int i = 0; i < n; i++)
   {
     for (int j = 0; j < n; j++)
       printf("  %f  ",L [i][j] );
      printf("\n\n");
   }
 
   printf("\n\n");
 
   for (int i  = 0; i < n; i++)
   {
     for (int j = 0; j < n; j++)
        printf("  %f  ",U [i][j] );
   printf("\n\n");
   }
   RealMatrix u =RealMatrix(n,n,U);
   RealMatrix l=RealMatrix(n,n,L);
 
   l.MatrixMultiply(u);
   printf("\n multi \n");
   l.printMatrix();
      printf("\n\n");
   return in;
}