Произведение матрицы на вектор

#include <iostream> 
using namespace std;
 
struct Vector
{
    float V1[100];
    float V2[100];
    float V3[100];
    float M[100][100];
    
};
/******************************************
Функция ReadM вводит с консольной клавиатуры элементы матрицы
******************************************/
void ReadM(Vector &x, int n)
{
 
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++) 
            cin >> x.M[i][j];
    return;
}
/******************************************
Функция ReadV вводит с консольной клавиатуры элементы вектора
******************************************/
void ReadV(Vector &x, int n)
{
    for(int j=0;j<n;j++)
        x.V1[j]=0;
    for(int i=0;i<n;i++)
    {
        cin >> x.V1[i];
    }
    return;
}
/*******************************************
Функция Proizved считает произведение матрицы на вектор
******************************************/
/*Vector Proizved(Vector x, int n)
{
    //Vector z;
    for(int j=0;j<n;j++)
    {
        x.V2[j]=0;
        for(int i=0;i<n;i++)
            x.V2[j]+=x.M[j][i]*x.V1[i];
    }
    //return x.V2;
}
*/
/***********************************
Функция Lenght считает длину вектора
**********************************/
float Lenght(float x[100], int n)
{
    float dlin=0;
    for(int i=0;i<n;i++)
        dlin = dlin + pow (x[i],2);
    return sqrt(dlin);
} 
/***********************************
Функция MaxL считает максимальную длину вектора
**********************************/
float MaxL(float a,float b, float c)
{
    float max;
    if(a>b)
        max=a;
    else
        max=b;
    if (c>max)
        max=c;
    return max;
} 
int main() 
{ 
    setlocale(0,""); //подключение русского языка
    int n,i,j; 
    float A[100][100]; //матрица А
    float x[100];      //вектор x
    float Ax[100];     //вектор Ax
    float Axx[100];    //вектор Ax-x
    Vector s,p,k;
    cout<<"Введите размерность матрицы : "<<endl;
    cin>>n;
    cout<<"Введите элементы матрицы : "<<endl;
    ReadM(s,n);   
    cout<<endl;
    cout<<"Введите координаты вектора : "<<endl;
    ReadV(p,n);
    //k.V2=Proizved(p,n);
    cout<<k.V1[0]<<endl;
    cout<<k.V1[1]<<endl;
    cout<<k.V1[2]<<endl;
 
    //считаем вектор Ax
    for(j=0;j<n;j++)
    {
        k.V2[j]=0;
        for(i=0;i<n;i++)
            k.V2[j]+=k.M[j][i]*k.V2[i];
    }
    //считаем вектор Ax-x
    for(i=0;i<n;i++)
        k.V3[i]=k.V2[i]-k.V1[i];
    cout<<endl;
    
    float a,b,c;//длины векторов
    //находим длины векторов x, Ax, Ax-x
    a=Lenght(k.V1,n);
    b=Lenght(k.V2,n);
    c=Lenght(k.V3,n);
    
    //выводим длины на экран
    cout<<"Длина вектора x = "<<a<<endl;
    cout<<"Длина вектора Ax = "<<b<<endl;
    cout<<"Длина вектора Ax-x = "<<c<<endl;
    //находим вектор максимальной длины
    float max=MaxL(a,b,c);
    cout<<"Максимальная длина = "<<max<<endl;
 
    system("pause"); 
    return 0;
}