Даны две матрицы. Необходимо их сложить,умножить на значение и см.тз

#include <iostream>
#include <math.h>
using namespace std;
// Создание матрицы
//
double ** Create( size_t n, size_t m ) {
    double ** M = new double * [n];
    for ( size_t i = 0; i < n; ++i ) {
        M[i] = new double [m];
    }
    return M;
}
// Удаление матрицы
//
void Free( double ** M, size_t n ) {
    for ( size_t i = 0; i < n; ++i ) {
        delete [] M[i];
    }
    delete [] M;
}
//---- ввод матрицы--------------------------------------
//
void Input( double ** M, size_t n, size_t m ) {
    for ( size_t i = 0; i < n; ++i ) {
        for ( size_t j = 0; j < m; ++j ) {
            cout << "M[" << i << "][" << j << "] = ";
            cin >> M[i][j];
        }
    }
}
//-------- Печать матрицы ------------------------------------------------
//
void Print( double ** M, size_t n, size_t m ) {
    for ( size_t i = 0; i < n; ++i ) {
        for ( size_t j = 0; j < m; ++j ) {
            cout<<M[i][j]<<' ';
        }
        cout<<endl;
    }
}
//вычисление определителя методом гаусса
double gaus_det(double **mass,int n)
{
int i,j,k;
 double det=1;
//прямой ход
for(i=0;i<n;i++)
{
for(j=i+1;j<n;j++)
{
if(mass[i][i]==0)
{
   for(int z=0;z<n;z++)
   {
       if(z!=i and mass[z][i]!=0)
       {
           for(int w=0;w<n;w++)
           {
               mass[i][w]+=mass[z][w];
           }
           continue;
       }
   }
}
 
double b=mass[j][i]/mass[i][i];
for(k=i;k<n;k++)
mass[j][k]=mass[j][k]-mass[i][k]*b;
}
det*=mass[i][i];//вычисление определителя
}
return det;
}
 
 
int main()
{
 
 
     size_t m;
    cout << "m=?" << endl;
   cin >> m;
   double** M1=Create(m,m);
   double** M2=Create(m,m);
   double** Result=Create(m,m);
   cout << "Vvedite pervuyu matrizu" << endl;
   Input(M1,m,m);
    cout << "Vvedite vtoruyu matrizu" << endl;
   Input(M2,m,m);
   for(int i=0;i<m;i++)
  for(int j=0;j<m;j++)
    Result[i][j] = M1[i][j] + M2[i][j];
 
    cout << "Summa\n";
    Print(Result,m,m);
 
    double x=1.5,a=4,p=1.6;
    double y=(1/exp(x)*sqrt(pow(x,2)-pow(a,2)))-exp(x)*sqrt(pow(x,2)-pow(a,2))+14.9*pow(p,2);
 
      for(int i=0;i<m;i++)
  for(int j=0;j<m;j++)
    Result[i][j] *= y;
 
     cout << "Matriza posle umnojeniya\n";
    Print(Result,m,m);
    int count=0;
    bool NullElements;
  for(int i=0;i<m;i++)
{
    NullElements=false;
  for(int j=0;j<m;j++)
  {
   if( Result[i][j] == 0)
   {
       NullElements=true;
   }
  }
  if(NullElements==false)
  {
      count++;
  }
}
 
cout << "Kolichestvo Ryadov bez nulevih elementov=" << count << endl;
 
double det=gaus_det(M1,m);
cout << "Opredelitel=" << det << endl;
 
    return 0;
}

#include <iostream>
#include <cmath>
 
using namespace std;
class matrix
{
public:
    float **a;
    int columns;
    int lines;
public:
    matrix()
    {
        a=NULL;
        columns=0;
        lines=0;
    }
 
    void summ(matrix &obj1,matrix &obj2,float y)
    {
        matrix f_summ;
        matrix mul;
        cout<<"Сумма матриц"<<endl;
        f_summ.columns=obj1.columns;
        f_summ.lines=obj1.lines;
        f_summ.a=new float*[f_summ.lines];
        for(int i=0;i<f_summ.lines;i++)
        {
            f_summ.a[i]=new float[f_summ.columns];
            for(int j=0;j<f_summ.columns;j++)
            {
                f_summ.a[i][j]=obj1.a[i][j]+obj2.a[i][j];
                cout<<f_summ.a[i][j]<<"  ";
            }
            cout<<endl;
        }
        cout<<"Умножение на У"<<endl;
        mul.columns=f_summ.columns;
        mul.lines=f_summ.lines;
        mul.a=new float*[mul.lines];
        for(int i=0;i<mul.lines;i++)
        {
            mul.a[i]=new float[mul.columns];
            for(int j=0;j<mul.columns;j++)
            {
                mul.a[i][j]=f_summ.a[i][j]*y;
                cout<<mul.a[i][j]<<"  ";
            }
            cout<<endl;
        }
 
    }
 
 
    friend  istream &operator >>(istream &in,matrix &obj);
    friend  ostream &operator <<(ostream &out,matrix &obj);
};
 
istream &operator >>(istream &in,matrix &obj)
{
    cout<<"Введите данные"<<endl;
    cout<<"Введите число столбцов: ";
    cin>>obj.columns;
    cout<<"Введите число линий: ";
    cin>>obj.lines;
    obj.a=new float*[obj.lines];
    for(int i=0;i<obj.lines;i++)
    {
        obj.a[i]=new float[obj.columns];
        for(int j=0;j<obj.columns;j++)
        {
            cout<<"a["<<i+1<<"]["<<j+1<<"]=";
            in>>obj.a[i][j];
        }
    }
    return in;
}
 
ostream &operator <<(ostream &out,matrix &obj)
{
 
    cout<<"Вывод данных"<<endl;
    for(int i=0;i<obj.lines;i++)
    {
        for(int j=0;j<obj.columns;j++)
        {
            out<<"a["<<i+1<<"]["<<j+1<<"]="<<obj.a[i][j]<<" ";
        }
        out<<endl;
    }
    return out;
}
 
float formula()
{
    float y;
    y=(1/(exp(1.5)/*sqrt(pow(1.5,2)-pow(4,2))*/))-exp(1.5)/*sqrt(pow(1.5,2)-pow(4,2))*/+14.9*pow(1.6,2);
    return y;
}
 
void notziro(matrix obj)
{
    int count=obj.lines;
    for(int i=0;i<obj.lines;i++)
    {
        for(int j=0;j<obj.columns;j++)
        {
            if(obj.a[i][j]==0)
            {
                count--;
            }
            break;
        }
    }
    cout<<"Кол-во рядов не содержащих ноль "<<count;
    cout<<endl;
}
 
 
int main()
{
    matrix obj1,obj2,s;
    float y;
    cin>>obj1;
    cout<<endl;
    cin>>obj2;
    y=formula();
    s.summ(obj1,obj2,y);
    notziro(obj1);
 
    return 0;
}