Помогите исправить ошибки в моей программе LU-разложение за Гауссом(Dev-C++)

#include <stdio.h>
#include <conio.h>
#include <iostream>
#include <stdlib.h>
#include <time.h>
#include <locale.h>
#include <iomanip>//setw()
#include <cstdlib>// rand (), srand ()
#include <windows.h>
using namespace std;
 
// пользовательские функции:
// меняем местами строки в массиве с коэфицентами уравнения
 
void exchange(double **mas, int n,int from, int to){
    for(int i=0; i < n; i++){
        mas[from][i] = mas[from][i] + mas[to][i];
        mas[to][i] = mas[from][i] - mas[to][i];
        mas[from][i] = mas[from][i] - mas[to][i];
    }
}
// определяем ближайшую строку с не 0 диоганалью вниз
int get_n_no_empty(double **MA, int n, int begin){
    for(int i=begin+1; i < n; i++){
        if(MA[i][i] != 0){
            return i;
        }
    }
    cout << "Данная система не решается методом Гаусcа" << endl;
    getch();
    exit(1);
}
 
// функция ждет пока пользователь прочтет и нажмет на кнопку, не обязательная функция
void pause(){
    cout << "Нажмите любую кнопку для продолжения" << endl;
    getch();
}
 
 
 
int main()
{
setlocale(LC_ALL, "Russian"); //чтоб буквы были русскими 
 
 int n,i,j; //размер матрицы
 double **MA,**U,**L // матрицы 
 
// переменной razmer присвоим размер матрицы */
cout<<"Введите размерность матрицы МА:\n";
cin>>razmer;
 
    // создание динамических переменных
    MA = new double*[n]; //память для матрицы
    U = new double*[n];
    L = new double*[razmer];
    
    for(int i=0; i < razmer; i++){
        MA[i] = new double[razmer];
  
//присвоим значение матрице MA 
            for (int j=0;j<n;j++)
            {time_t temp = time(NULL);//получаем время
             srand(temp); //инициализируем генератор случайныйх чисел полученным временем
             MA[i][j] = rand()%(-11)-5; //заполняем матрицу случайным образом от -5 до 5
        }
        }
        //вывод матрицы на экран
        cout<<"Сформированная матрица:\n";
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
                cout<<setw(5)<<*(*(MA+i)+j)<<ends;//setw(5)-ширина поля вывода з iomanip.h
            cout<<endl;
        }
 // начинаем решать
    // прямой проход
    
    for(int k=0; k < n-1; k++){//думаю,где-то здесь есть проблемы
        // вычитаем по строкам
        for(int m=k+1; m < n; m++){
            // если на диагонали элемент = 0, то поменяем местами строки
            if(MA[m][m] == 0) exchange(MA, n, m, get_n_no_empty(MA, n, m));
            
            double M[m][k]=MA[m][k]/MA[0][k];
            cout<<"Коэфициент M["m"]["k"]:"<<M[m][k]<<endl;
            //вычисление новых коэфицентов уравнения
           
            for(int z=0; z < n; z++){
                MA[m][z] = MA[m][k] - MA[0][k] * M[m][k];
             //вывод матрицы на экран
        cout<<"Сформированная эквивалентная матрица :\n";
        //Теперь  прировняем эквивалентную матрицу MA[m][z] к U
        for(int m=k+1; m < n; m++)
        {U[m] = new double[n];
            for(int z=0; z < n; z++){
                MA[m][z]=U[m][z];
            cout<<endl;
        }
            }
        }
    }  cout << "L берём как одиничную матрицу:\n";
        for ( i = 0; i < n; i ++ ){
            for ( j = 0; j < n; j ++ ){
                cout << matrixL[i][j] << " ";
            }
            cout << "\n";
        }
/*В нижнюю /\-ю матрицу L подставим те коэфициенты, которые которые использовались
для построения верхней треугольной матрицы U- как это сделать ещё не придумала(запуталась в коэф.*/
     
getch();
system("pause");
return 0;
}