Работа с матрицами

#include<iostream>
#include "math.h"
 
using namespace std;
 
double eps;//погрешность
 
class Matrix{
    double **a;//матрица это двумерный массив, нам будет удобней работать с квадратными матрицами, поэтому делаем их такими
    int n;//размер матрицы
public:
    Matrix(int k){
        n=k;
        a=new double *[n];
        for(int i=0;i<n;i++){
            a[i]=new double [n];
            for(int j=0;j<n;j++){
                a[i][j]=i+(j*(i/2));
            }
        }
    };//конструктор
    ~Matrix(){//деконструктор
    };
    friend const Matrix operator+(const Matrix &x, const Matrix &y);//перегрузка для сложения двух матриц
    friend const Matrix operator*( const Matrix& x, const double z);//перегрузка для умножения матрицы на число
    friend const Matrix operator*(const Matrix &x, const Matrix &y);//перегрузка для умножения матрицы на матрицу
    void transp();//функция транспонирования
    double d();//функция которая считает дискрименант
    void obr();//функция обращающая матрицу
    void Show();//функция выводит матрицу на экран
};
void Matrix::Show(){
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            cout<<a[i][j]<<" ";
        }
        cout<<"\n";
    }
    cout<<"\n";
}
void Matrix::transp(){
    double **b;//создаем вспомогающую матрицу
    b=new double *[n];
    for(int i=0;i<n;i++){
        b[i]=new double [n];
        for(int j=0;j<n;j++){
            b[i][j]=a[j][i];//транспонируем матрицу
        }
    }
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            a[i][j]=b[i][j];//меняем элементы главную матрицу на транспонированую
        }
    }
    for(int i=0;i<n;i++){
        delete []b[i];
    }
    delete []b;
}
 
double Matrix::d(){
    double D=0;//сюда записываем результат
    double **b;//создаем вспомагательный массив
    b=new double *[n];
    for(int i=0;i<n;i++){
        b[i]=new double [n+n-1];//создаем удобную для просчета определителя матрицу
        for(int j=0;j<n+n-1;j++){
            b[i][j]=a[i][j%n];//заполняем матрицу b
        }
    }
    for(int i=0;i<n;i++){
        double temp=1;
        for(int j=0;j<n;j++){
            temp=temp*b[j%n][j+i];//просчитываем главные диагонали
        }
        D=D+temp;
    }
    for(int i=0;i<n;i++){
        double temp=1;
        for(int j=0;j<n;j++){
            temp=temp*b[n-1-(j%n)][j+i];//просчитываем побочные диагонали
        }
        D=D-temp;
    }
    for(int i=0;i<n;i++){
        delete []b[i];
    }
    delete []b;
    return D;
}
void Matrix::obr(){
    double r=d();
    if(abs(r)<eps){
        cout<<"obratnoi matrici net!\n";
    }else{
        Matrix M(n);//минор матрицы
        if(n>2){//если размер матрицы больше 2 то, чтобы найти минор нужно считать определители, для этого создаем новую вспомогательную матрицу
            Matrix temp(n-1);
            int ik=0,jk=0;
            for(int i=0;i<n;i++){//циклы для 
                for(int j=0;j<n;j++){//просчета минора
                    for(int k=0;k<n;k++){//циклы для
                        for(int l=0;l<n;l++){//занесения чисел в temp
                            if(k!=i&&l!=j){//если элемент не из вычеркнутого строки или столбца
                                temp.a[ik][jk]=a[k][l];//занести в temp
                                if(jk==n-2){//перенос места под элемент в temp-е
                                    ik++;
                                    jk=0;
                                }else jk++;
                                if(ik==n-1)break;
                            }
                        }
                    }
                    M.a[i][j]=temp.d()*((-1)*(i+j+2));//записываем минор для элемента [i][j]
                }
            }
        }else{
            M.a[0][0]=a[1][1];//иначе матрица
            M.a[1][0]=a[0][1]*(-1);//2Х2 и ее
            M.a[0][1]=a[1][0]*(-1);//можно сделать
            M.a[1][1]=a[0][0];//вручную
        }
        M.transp();//транспонируем
        Matrix v(n);
        v=M*(double(1)/r);
        for(int i=0;i<n;i++){
            for(int j=0;j<n;j++){
                a[i][j]=v.a[i][j];
            }
        }
    }
}
 
const Matrix operator+(const Matrix &x, const Matrix &y){
    if(x.n!=y.n)cout<<"matrici ne ravni!\n";
    else{
        Matrix b(x.n);
        for(int i=0;i<x.n;i++){
            for(int j=0;j<x.n;j++){
                b.a[i][j]=x.a[i][j]+y.a[i][j];
            }
        }
        return b;
    }
}
const Matrix operator*( const Matrix& x, const double z){
    Matrix b(x.n);
    for(int i=0;i<x.n;i++){
        for(int j=0;j<x.n;j++){
            b.a[i][j]=x.a[i][j]*z;
        }
    }
    return b;
}
const Matrix operator*(const Matrix &x, const Matrix &y){
    Matrix b(x.n);
    if(x.n!=y.n)cout<<"matrici ne ravni!\n";
    else{
        for(int i=0;i<x.n;i++){
            for(int j=0;j<x.n;j++){
                b.a[i][j]=0;
                for(int k=0;k<x.n;k++){
                    b.a[i][j]=b.a[i][j]+x.a[i][k]*y.a[k][j];
                }
            }
        }
        return b;
    }
}
int main(){
    int n1, n2;
    cout<<"vvedite razmer pervoi matrici: ";
    cin>>n1;
    cout<<"vvedite pogreshnost: ";
    cin>>eps;
    Matrix Mat1(n1);
    cout<<"vvedite razmer vtoroi matrici: ";
    cin>>n2;
    Matrix Mat2(n2);
    Mat1.Show();
    Mat2.Show();
    Mat2=Mat1;
    Mat1.obr();
    Mat1.Show();
    Mat1.transp();
    Mat1.Show();
    Matrix v(n1);
    v=Mat1*Mat2;
    v.Show();
    v=Mat1+Mat2;
    v.Show();
    v=Mat1*0.1;
    v.Show();
    return 0;
}