Классы для работами с матрицами

#include "stdafx.h"
#include <iostream>
using namespace std;
class Matrix{
private:
int n,m;
//размеры матрицы
double **a, **c;
// указатели
public:
Matrix::Matrix(int n, int m){
int i,j;
Matrix::n=n;
Matrix::m=m;
a=new double*[n]; // выделение памяти под указатели
for(i=0;i<n;i++)
a[i]=new double[m]; // выделение памяти под элемент
c=new double*[n];
for(int i=0; i<n; i++)
c[i]=new double[m];
// Заполнение массива элементами с клавиатуры
{puts("\nВведите элементы матрицы:");
for(i=0;i<n;i++)
for(j=0;j<m;j++)
{cout<<"a["<<i<<"]["<<j<<"]=";
cin>>a[i][j];}
cout<<endl;}
//};
Matrix::~Matrix()
{for(int i=0;i<n;i++)
delete []a[i];// освобождение памяти
delete []a;};
//Операция присваивания
Matrix &Matrix::operator =(const Matrix &A)
{for(int i=0; i<n; i++)
for(int j=0; j<m; j++)
a[i][j]=A.c[i][j];
return *this;}
//Сложение матриц
Matrix & Matrix::operator +(Matrix &A)
{for(int i=0; i<n; i++)
for(int j=0; j<m; j++)
c[i][j]=a[i][j]+A.a[i][j];
return *this;}
//Вычитание матриц
Matrix & Matrix::operator -(Matrix &A)
{for(int i=0; i<n; i++)
for(int j=0; j<m; j++)
c[i][j]=a[i][j]-A.a[i][j];
return *this;}
// Умножение матриц
Matrix & Matrix::operator *(Matrix &A)
{
double s=0;
for(int i=0;i<n;i++)
{
int z=0;
while(z<A.m)
{
c[i][z]=0;
for(int j=0;j<m;j++)
{
s=GetI(i,j)*A.GetI(j,z);
c[i][z]+=s;
}
z++;
}
}
return *this;
}
//Сравнение матриц
Matrix & Matrix::operator ==(Matrix &A)
{int f=0;
for(int i=0; i<n; i++)
{for(int j=0; j<m; j++)
if(a[i][j]!=A.a[i][j]){ cout<<"\nМатрицы не равны\n";f=1; break;}
if(f==1) break;}
if(f==0) cout<<"\nМатрицы равны\n";
return *this;}
//Транспонирование матрицы
void Matrix::Transp()
{c=new double*[m];
for(int i=0;i<m;i++)
c[i]=new double[n]; // захват памяти под элемент
for(int i=0; i<n; i++)
for(int j=0; j<m; j++)
c[j][i]=a[i][j];
for(int i=0;i<n;i++)
delete []a[i];// освобождение памяти
delete []a;
a=new double*[m];
for(int i=0;i<m;i++)
a[i]=new double[n]; // захват памяти под элемент
for(int i=0; i<m; i++)
for(int j=0; j<n; j++)
a[i][j]=c[i][j];
Matrix::n=m;
Matrix::m=n;
};
void main()
{
};