Матричная экспонента

Известно, что экспонента может быть представлена в виде бесконечного ряда: $$exp\left(z \right)=\sum_{n=0}^{\infty }{\frac{z^{n}}{n!}}$$ если аргументом экспоненты считать матрицу А, подставить эту матрицу в разложение, то получится матричная экспонента. Кому это может понадобиться не известно, но здесь в листинге демонстрируется как матрица может быть представлена в виде класса, а потому код получается для такого сложного матричного ряда вполне нормальным. Ну или почти нормальным.

#include 
using namespace std;
//Количество членов ряда:
const int N=100;
//Класс для реализации матриц:
class Matrix{
public:
double a[2][2];
//Метод для умножения матриц:
Matrix mult(Matrix obj){
Matrix T;
int i,j,k;
for(i=0; i < 2; i++)
 for(j=0; j < 2; j++){
 T.a[i][j]=0;
 for(k=0; k < 2; k++)
 T.a[i][j]+=a[i][k]*obj.a[k][j];
 }
return T;}
//Метод для суммирования матриц:
Matrix mSum(Matrix obj){
Matrix T;
int i,j;
for(i=0;i<2;i++)
 for(j=0;j<2;j++)
 T.a[i][j]=a[i][j]+obj.a[i][j];
return T;}
//Деление матрицы на число:
Matrix div(double x){
int i,j;
for(i=0;i<2;i++)
 for(j=0; j < 2; j++)
 a[i][j]/=x;
return *this;}
//Отображение матрицы:
void mShow(){
printf("%8f%s%8f%s",a[0][0]," ",a[0][1],"\n");
printf("%8f%s%8f%s",a[1][0]," ",a[1][1],"\n");
}};
//Матричная экспонента:
Matrix mExp(Matrix obj){
int i;
Matrix E,T;
E.a[0][0]=1;
E.a[1][1]=1;
E.a[0][1]=0;
E.a[1][0]=0;
T=E;
for(i=1; i <= N; i++){
 T=T.mult(obj);
 T=T.div(i);
 E=E.mSum(T);}
return E;}
int main(){
Matrix A;
A.a[0][0]=1;
A.a[0][1]=-1;
A.a[1][0]=2;
A.a[1][1]=1;
A.mShow();
cout << endl;
mExp(A).mShow();
return 0;}