Код более понятней сделать

#include "stdafx.h"
#include <iostream>
using namespace std;
template <class Type> class Array
{
    friend ostream& operator<<(ostream& out, const Array& a)
    {
        
        for (int i = 0; i < a.nomberOfRows; i++)
        {
            out<<"\n";
            for(int j=0;j<a.nomberOfColumns;j++)
                out << a.pa[i][j] << ' ';
        }  
        return out;
 
    }
 
    friend istream& operator>>(istream& in, Array& a)
    {
        for (int i = 0; i < a.nomberOfRows; i++)
            for(int j=0;j<a.nomberOfColumns;j++)
                in >> a.pa[i][j];
        return in;
    }
 
public:
    Type **pa;
    int  nomberOfRows;
    int nomberOfColumns;
 
    class OutOfBounds
    {
        int raws;
        int columns;
    public:
        OutOfBounds(int a, int b) : raws(a) , columns(b) { }
        int getIndex() const { return raws; }
        int getIndex2() const {return columns;}
    };
 
    class IncorrectSize{
        int raws,raws1,columns,columns1;
    public:
        IncorrectSize(int a,int b, int c,int d) : raws(a),columns(b) ,raws1(c), columns1(d){}
        int getR() {return raws;}
        int getR1() {return raws1;}
        int getC() {return columns;}
        int getC1() {return columns1;}
    };
 
    Array(){pa=0;nomberOfRows=0;nomberOfColumns=0;}
    Array(int raws, int columns, Type **a);
    Array(const Array& arr);
    ~Array() { 
        if (pa) {
            for (int i=0;i<nomberOfRows;i++)
                delete [] pa[i];
            delete []pa;
        }
    }
    const Array& operator=(const Array& a);
    Type&        operator()(int raws, int columns);
 
    Array(int r,int c)
    {
        nomberOfRows=r;
        nomberOfColumns=c;
        pa = new Type * [ r];
        for(int i=0;i<r;i++)
            pa[i]=new Type[c];
        for(int i=0;i<r;i++)
            for (int j=0;j<c;j++)
                pa[i][j]=0;
    }
 
    int getRows(){
        return nomberOfRows;
    }
    int getColumns(){
        return nomberOfColumns;
    }
    Type** getPa(){return pa;}
    
 
    const Array<Type> operator-( const Array& a ){
        Array<Type> c(a.nomberOfRows, a.nomberOfColumns);
        if(a.nomberOfRows==nomberOfRows && a.nomberOfColumns==nomberOfColumns)
        {  
        for(int i=0;i<a.nomberOfRows;i++)
            for(int j=0;j<a.nomberOfColumns;j++)
                c.pa[i][j]=pa[i][j]-a.pa[i][j];
        return c;
        }
        else
        {cout <<"Matrix cant be minused, because of different sizes"<<'\t';
        return c;
        }
    }
 
    const Array<Type> operator*( const Array& a ){
        Array<Type> c(nomberOfRows, a.nomberOfColumns);
        if(nomberOfColumns!=a.nomberOfRows)
            throw Array<Type>::IncorrectSize(nomberOfRows,nomberOfColumns,a.nomberOfRows,a.nomberOfColumns);
 
            for (int i = 0; i < nomberOfRows; i++)
                for (int j = 0; j < a.nomberOfColumns; j++)
                    for (int k = 0; k < a.nomberOfRows; k++)
                        c.pa[i][j] += pa[i][k] * a.pa[k][j];
            return c;     
        
        
 
    }
}; 
 
template <class Type>
const Array<Type>& Array<Type>::operator=(const Array& a)
{
    if (&a != this)
    {
        if (pa)
        ~Array();
        nomberOfRows= a.nomberOfRows;
        nomberOfColumns=a.nomberOfColumns;
        for (int i = 0; i < nomberOfRows; i++)
            delete[] pa[i];
        delete[] pa;
        pa = new Type* [nomberOfRows];
        for (int i = 0; i < nomberOfRows; i++)
            pa[i]=new Type [nomberOfColumns];
 
        for(int i=0;i<nomberOfRows;i++)
            for(int j=0;j<nomberOfColumns;j++)
                pa[i][j]=a.pa[i][j];
    }
    return *this;
}
 
template <class Type>
Array<Type>::Array(int raws, int columns , Type **a)
{
    nomberOfRows=raws;
    nomberOfColumns=columns;
    pa = new Type * [ raws];
    for(int i=0;i<raws;i++)
        pa[i]=new Type[columns];
 
    for (int i = 0; i < nomberOfRows; i++)
        for(int j=0;j<nomberOfColumns;j++)
            pa[i][j]=a[i][j];
}
 
template <class Type>
Array<Type>::Array(const Array& arr)
{
    nomberOfRows=arr.nomberOfRows;
    nomberOfColumns=arr.nomberOfColumns;
    pa = new Type* [nomberOfRows];
    for (int i = 0; i < nomberOfRows; i++)
        pa[i]=new Type [nomberOfColumns];
    for(int i=0;i<nomberOfRows;i++)
        for(int j=0;j<nomberOfColumns;j++)
            pa[i][j]=arr.pa[i][j];
}
 
template <class Type>
Type& Array<Type>::operator()(int raws,int columns)
{
    if (raws < 0 || raws >= nomberOfRows || columns<0 || columns >=nomberOfColumns)
        throw OutOfBounds(raws , columns);
    return pa[raws][columns];
}
 
template<class Type>
 Array<Type> operator+(const Array<Type>& a, const Array<Type> &b)
    {
        Array<Type> c(a.nomberOfRows,a.nomberOfColumns);
    if(a.nomberOfRows!= b.nomberOfRows || a.nomberOfColumns!=b.nomberOfColumns)
        throw Array<Type>::IncorrectSize(a.nomberOfRows,a.nomberOfColumns,b.nomberOfRows,b.nomberOfColumns);
        
        for(int i=0;i<a.nomberOfRows;i++)
            for(int j=0;j<a.nomberOfColumns;j++)
                c.pa[i][j]=a.pa[i][j]+b.pa[i][j];
 
        return c;
    
 }
 
template <typename Type>
Type max( Array<Type> &a ){
    Type **p=a.getPa();
    Type max=p[0][0];
    for(int i=0;i<a.getRows();i++)
        for(int j=0;j<a.getColumns();j++)
            if(max<p[i][j]) 
                max=p[i][j];
    return max;
}
 
 
int main(int argc, char* argv[])
{
    int raws;
    int columns;
    cout <<" Input nomber of raws"<< '\n';
    cin>> raws;
    cout <<"Input nomber of columns" <<'\n';
    cin>> columns;
    Array <double> a(raws,columns);
    cout << "Now input elements of Matrix by rows"<<'\n';
    cin>>a;
    int raws1,columns1;
    cout <<" Input nomber of raws of second matrix"<< '\n';
    cin>> raws1;
    cout <<"Input nomber of columns of second matrix" <<'\n';
    cin>> columns1;
    Array<double> b(raws1,columns1);
    cout << "Now input elements of  second Matrix by rows"<<'\n';
    cin>>b;
    cout << "Inputed data:";
    cout << a << endl;
    cout <<b<<endl;
    try{
    cout <<"a+b="<<endl <<a+b <<endl;
    cout <<"a-b="<<endl<<a-b<<endl;
    }
    catch(Array<double>::IncorrectSize e)
    {
        cout << "Operations + and - cant be perfomed, because of different sizes"<<endl<< e.getR()<<"x" << e.getC() <<endl <<e.getR1() <<"x" << e.getC1()<<endl;
    }
    try
    {
    cout <<"a*b=" <<endl<< a*b<<endl;
    }
    catch(Array<double>::IncorrectSize e)
    {
        cout << "Operation * cant be perfomed, because of different sizes"<<endl<< e.getR()<<"x" << e.getC() <<endl <<e.getR1() <<"x" << e.getC1()<<endl;
    }
    cout <<"max element of A(first) matrix =";
    double i =max <double>(a);
    cout <<i<<endl;
    try{
        cout << a(11,11);
    }
    catch(Array<double>::OutOfBounds e)
    {
        cout << "Wrong data, matrix out of bounds on such indexes"<<e.getIndex() << ' ' <<e.getIndex2();
    }
    system("pause");
    return 0;
}