#include <iostream>
#include <cstdio>
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("%8%s%8f%s",a[0][0]," ",a[0][1], "\n");
    printf("%8%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=0;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;
}