Вычислить вектор P=A^n*X, где A-матрица, X-вектор

#include <iostream>
#include <string>
#include <math.h>
 
using namespace std;
 
const int N=10;
 
void InMatr(double x[N][N], int k)
{
    int i,j;
    for(i=0;i<k;i++){
        for(j=0;j<k;j++){
            cout<<"("<<i+1<<")"<<"("<<j+1<<")"<<" "<<"элемент матрицы=";
            cin>>x[i][j];
        }
    }
    cout<<"\n";
}
void InVect(double x[N], int k)
{
    int i;
    for(i=0;i<k;i++){
        cout<<i+1<<" "<<"элемент вектора=";
        cin>>x[i];
    }
    cout<<"\n";
}
void CMatr(double x[N][N], int k, double y[N][N])
{
    int i, j;
    for(i=0;i<k;i++){
        for(i=0;i<k;i++){
            y[i][j]=x[i][j];
        }
    }
}
void UMatr(double x[N][N], double y[N][N], int k,double z[N][N])
{
    int i,j,l;
    double S;
    for(i=0;i<k;i++){
        for(j=0;j<k;j++){
            S=0;
            for(l=0;l<k;l++){
                S=S+x[i][l]*y[l][j];
                z[i][j]=S;
            }
        }
    }
}
void PowMatr(double x[N][N],int k, int l, double y[N][N])
{
    double yx[N][N];
    int i;
    CMatr(y,k,x);
    for(i=2;i<=l;i++){
        UMatr(y,x,k,yx);
        CMatr(y,k,yx);
    }
}
void UMatrVec(double x[N][N],double y[N], int k, double z[N][N])
{
    int i,j,l;
    double S;
    for(i=0;i<k;i++){
        for(j=0;j<k;j++){
            S=0;
            for(l=0;l<k;l++){
                S+=x[i][l]*y[l];
                z[i][j]=S;
            }
        }
    }
}
 
void OutMatr(double x[N][N], int k)
{
    int i,j;
    for(i=0;i<k;i++){
        cout<<"\n";
        for(j=0;j<k;j++)
            cout<<x[i][j]<<" ";
    }
    cout<<"\n";
}
 
void OutVec(double x[N], int k)
{
    int i,j;
    cout<<"Исходный вектор:"<<endl;
    for(i=0;i<k;i++){
        cout<<x[i]<<" ";
    }
    cout<<"\n";
}
 
int main()
{
    setlocale(LC_ALL, "Russian");
    double a[N][N],b[N],z[N][N],o[N][N],n,k,q,s;
    cout << "Введите размерность матриц квадратных матриц NxN\n"<<"Размер =";
    cin>>n;
    cout<<"\n";
    cout << "Введите длинну векторов\n"<<"Длинна=";
    cin>>k;
    cout<<"\n";
    cout << "Введите степень в которую нужно возвести матрицу\n"<<"Степень=";
    cin>>q;
    cout<<"\n";
    cout << "Введите элементы матрицы A:\n";
    InMatr(a, n);
    cout << "Введите элементы вектора B:\n";
    InVect(b, k);
    CMatr(a,n,z);
    PowMatr(a,n,q,z);
    UMatrVec(z,b,n,o);
    OutMatr(o,n);
    cout << "\nРезультат программы:\n";
    OutMatr(a,n);
}