Матрица,нужен исходник

class square_matrix{
    private:
            double** matrix;
            int dimm;
            
            void create(){
                matrix = new double* [dimm];
                for (int i=0; i<dimm; i++)
                    matrix[i]= new double [dimm];
            }
    
    public:
            square_matrix(int n): dimm(n) {create();}
            
            ~square_matrix(){
                for (int i=0; i<dimm; i++)
                    delete[] matrix[i];
                delete[] matrix;
            }
            
            double &E(int i, int j){
            return matrix[i][j]; // E(lement)
            }
            
            int get_dimm()  { return dimm; }
                
};
 
 
#include<iostream>
void show(square_matrix &obj);          // Отображает матрицу obj 
void generate(square_matrix &obj);      // Заполняет матрицу obj случайными числами
 
void get_max(int &i, int &j , square_matrix &obj,int n);    
                                        // Кладёт в i,j координаты максимального из 
                                        // элементов матрицы obj не стоящих на первых n
                                        // позициях главной диагонали. Ключевая функция.
 
void swap(double &a, double &b);        // обменивает значениями переменные а и b.
                                        
int main()
{
    int dimm = 3;
    square_matrix M(dimm);
    
    generate(M);
    show(M);
    
    for(int n=0; n<dimm;n++){
        int row,col;
        get_max(row, col, M, n);
        swap(M.E(row,col),M.E(n,n));
    }
    
    
    show(M);
    
    
}
 
 
void show(square_matrix &matrix){
    int dimm = matrix.get_dimm();
    
    for(int i=0;i<dimm;i++){
        for(int j=0;j<dimm;j++)
            std::cout<<matrix.E(i,j)<<"\t";
        std::cout<<std::endl;
    }
    std::cout<<std::endl;
}
 
void generate(square_matrix &matrix){
    int dimm = matrix.get_dimm();
    
    srand(time(NULL));
    for (int i=0; i<dimm;i++)
        for (int j=0;j<dimm;j++)
            matrix.E(i,j) = rand()%100+1;
    
}
 
void swap(double &a, double &b){
    double tmp = a;
    a = b;
    b = tmp;
}
 
void get_max(int &row, int &col ,square_matrix &matrix,int n){
// Кладёт в row,col координаты максимального из 
// элементов матрицы matrix не стоящих на первых n
// позициях главной диагонали. Ключевая функция.
    int dimm = matrix.get_dimm();
    double max = matrix.E(dimm-1,dimm-1);
    
    for(int i = 0; i<dimm; i++){
        for(int j = 0; j<dimm; j++){
            
            if((i!=j)||(i>=n))
                if(matrix.E(i,j)>max){
                    row = i;
                    col = j;
                    max = matrix.E(i,j);
                }
        }
    }
}