приведение матрицы к треугольному виду

#include <ctime>
#include <iostream>
#include <math.h>
#include <cstdlib>
using namespace std;
 
int main() {
    int mu=0;
    int razm1,razm2;
    cout<<"ENTER razm1 AND razm2"<<endl;
     while ((!(cin >> razm1>>razm2))||(razm1<0)||(razm2<0))
  {
        cin.clear();
        while (cin.get()!='\n');
        cout << "error! enter razm1 and razm2 again";
  }
  double **ptrarray = new double*[razm1];
  for (int i=0;i<razm2;i++) ptrarray[i]=new double(razm2);
  cout<<"enter matr";
  for(int i=0;i<razm1;i++)
      for (int j=0;j<razm2;j++)
      cin>>ptrarray[i][j]; 
 for (int j=0;j<razm2;j++)
   
      for (int i=j+1;i<razm1;i++)
      {
         mu=ptrarray[i][j]/ptrarray[j][j];
         for (int l=0;l<razm2;l++)
            ptrarray[i][l]-=ptrarray[j][i]*mu;
         }
         for (int i=0;i<razm1;i++)
         for (int j=0;i<razm2;j++)
         cout<<ptrarray[i][j]<<endl;
        
    
    system("pause");
    return 0;
}

#include <iostream>
 
using namespace std;
void print(int n, int **a)
{
    for(int i = 0; i < n; ++i)
    {
        for(int j = 0; j < n; ++j)
           cout << a[i][j] << "\t";
    cout << endl;
    }
}
int main()
{
setlocale(LC_ALL, "Russian");
cout<<"введите размер матрицы"<<endl;
int n;
cin>>n;
cout<<"введите элементы матрицы"<<endl;
    int **a=new int *[n];
    for(int i = 0; i < n; ++i)
       a[i]= new int [n];
for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
    {
      cout<<"A["<<i<<","<<j<<"]";
      cin>>a[i][j];
    }
cout << "исходная матрица" << endl;
print(n,a);
 for(int j=0;j<n;j++)
    for(int i=n-1;i>j;i--)
       for (int r=0;r<n;r++)
         a[i][r]=a[i][r]-((a[i][j]/a[i-1][j])*a[i-1][r]);
cout << "треугольная матрица" << endl;
print(n,a);
system("pause");
return 0;
}