Вычислить косинус матрицы в производном классе с помощью перегруженных операций из базового

#include <iostream>
#include <ctime>
#include <cstdlib>
#include <cmath>
 
using namespace std;
 
class matrix{
    protected:
 
        int mtr[5][5];
 
    public:
        
 
        matrix(){
            for(int i=0; i<=4; i++)
                for(int j=0; j<=4; j++)
                    mtr[i][j]=rand()%10;
        }
        
        matrix(int a){
        for(int i=0; i<=4; i++)
                for(int j=0; j<=4; j++)
                    (i==j)?mtr[i][j]=1 : mtr[i][j]=0;
    }
        
        
        
    void show(){
        for(int i=0; i<=4; i++){
            cout<<endl;
            for(int j=0; j<=4; j++)
            cout<<mtr[i][j]<<" ";
        }
    }
    
 
    
     matrix operator*(int a){
            matrix C;
            for(int i=0; i<=4; i++)
                for(int j=0; j<=4; j++)
                    C.mtr[i][j]=mtr[i][j]*a;
        return C;           
    }
    
     matrix operator/(int a){
            matrix C;
            for(int i=0; i<=4; i++)
                for(int j=0; j<=4; j++)
                    C.mtr[i][j]=mtr[i][j]/a;
        return C;           
    }
    
    matrix operator+(matrix B){
        matrix C;
        for(int i=0; i<=4; i++)
                for(int j=0; j<=4; j++)
                    C.mtr[i][j]=mtr[i][j]+B.mtr[i][j];
        return C;               
    }
    
    matrix operator-(matrix B){
        matrix C;
        for(int i=0; i<=4; i++)
                for(int j=0; j<=4; j++)
                    C.mtr[i][j]=mtr[i][j]-B.mtr[i][j];
        return C;               
    }
    
    matrix operator*(matrix B){
        matrix C;
        for(int i=0; i<=4; i++)
                for(int j=0; j<=4; j++){
                int s=0;
                    for(int k=0; k<=4; k++)
                    s=s+mtr[i][k]*B.mtr[k][j];
                    C.mtr[i][j]=s;  
                }
                
        return C;
    }
};
 
class cosMatrix : public matrix{
public:
    cosMatrix(): matrix()
    {}
    cosMatrix(int a): matrix(a)
    {}
    cosMatrix operator=(matrix a) {
        cosMatrix C;
        for (int i = 0; i <= 4; i++)
            for (int j = 0; j <= 4; j++)
                mtr[i][j] = a.mtr[i][j];
        return *this;
    }
};
 
int main() {
    srand(time(NULL));
    cosMatrix m10(1), m11(1), m12;
    matrix m1(1),m2(1),m3=m1+m2,m4=m1-m2,m5;
    int c=2;
    m1.show();
    cout<<endl;
    m2.show();
    m1=m1*c;
    cout<<endl;
    m1.show();
    cout<<endl;
    m2=m2/2;
    m2.show();
    cout<<endl;
    m3.show();
    cout<<endl;
    m4.show();
    cout<<endl;
    m5.show();
    cout<<endl;
    m10.show();
    cout << endl;
    m11.show();
    cout << endl;
    m12.show();
    cout << endl;
    m12 = m10 + m11+m11;
    m12.show();
    cin.get();
    return 0;
}