Перемножить 2 матрицы

 //Код перемножения матриц:
for(i = 0; i < 4; i++)
{
for(j = 0; j < 4; j++)
{SMatrix[i][j] = 0;
for(k = 0; k < 4; k++)
{
SMatrix[i][j] += vhodMatrix[k][j]*sinapsMatrix[i][k];
}
}
}
 
//Вывод результирующей матрицы:
SMatrix[1][4] = 0;
for(i = 0; i < 4; i++)
{
for(j = 0; j < 4; j++)
{
cout << SMatrix[i][j] << ", ";
}
cout << endl;
}

for(i = 0; i < 4; i++)
{
   for(j = 0; j < 4; j++)
   {
        SMatrix[i][j] = vhodMatrix[i][j]*sinapsMatrix[i][j];
   }
}

#include <stdio.h>
#include <conio.h>
#include <math.h>
#include <windows.h> 
#include <iostream>
 
char bufRus[256];
 
char* Rus(const char* text) {
      CharToOem(text, bufRus);
      return bufRus;
      }
int i, j, k;      
using namespace std;
char h;
float n;
 
int main()
{
   float SMatrix[j][i]; 
   
   cout << Rus("Вывод значений матрицы синаптических весов нейронов W")<< endl;
float sinapsMatrix[4][4] = {{0.42, 1, 0.4, 0.5}, //Ввод матрицы W -синаптических весов нейронов
{0.09, 0, 0.3, 0.2},
{0.33, 0, 0.2, 0.2},
{0.17, 0, 0.1, 0.1}};
 
sinapsMatrix[1][4] = 0;
 
for(k = 0; k < 4; k++)
{
for(j = 0; j < 4; j++)
{
cout << sinapsMatrix[k][j] << ", ";
}
 
cout << endl;
}
 
cout << h << endl;
float vhodMatrix[4][4] = {{321.6, 380, -6.7, 80}, 
{381.6, 480, -7.8, 85},
{330.8, 550, -3.9, 80},
{302.3, 500, -5, 80}};
 
 
cout << Rus("Вывод значений матрицы входных параметров Х")<< endl;
 
 
vhodMatrix[1][4] = 0;
 
for(i = 0; i < 4; i++)
{
for(k = 0; k < 4; k++)
{
cout << vhodMatrix[i][k] << ", ";
}
 
cout << endl;
}
 
//Операция перемножения матрицы W и матрицы X 
 
for(i = 0; i < 4; i++)
{
for(j = 0; j < 4; j++)
{SMatrix[i][j] = 0;
for(k = 0; k < 4; k++)
{
SMatrix[i][j] += vhodMatrix[k][j]*sinapsMatrix[i][k];
}
}
}
 
//Вывод решения перемножения матрицы W и матрицы X - значений матрицы текущего состояния сети S=W*X "
cout << Rus("Расчет значений матрицы текущего состояния сети S=W*X ")<< endl;
SMatrix[i][j] = 0;
for(i = 0; i < 4; i++)
{
for(j = 0; j < 4; j++)
{
cout << SMatrix[i][j] << ", ";
}
cout << endl;
} 
 system("pause"); 
 return 0;
}

#include <iostream> 
#include <iomanip>
 
int main() 
{ 
    const int size = 4;  
    int i, j, z;
    double A[size][size] = { 0.42, 1, 0.4, 0.5, 
                            0.09, 0, 0.3, 0.2,
                            0.33, 0, 0.2, 0.2,
                            0.17, 0, 0.1, 0.1}; 
 
    double B[size][size] = {321.6, 380, -6.7, 80, 
                           381.6, 480, -7.8, 85,
                           330.8, 550, -3.9, 80,
                           302.3, 500, -5, 80}; 
     
 
    
    double C[size][size]; 
 
 
    for ( i = 0; i < size; i++) 
        for ( j = 0; j < size; j++) 
        { 
            C[i][j] = 0; 
            for ( z = 0; z < size; z++) 
                C[i][j] += A[i][z] * B[z][j]; 
        }  
         
        for ( i = 0; i < size; i++) 
            for ( j = 0; j < size; j++) 
            { 
                if( j % size == 0) 
                    std::cout << std::endl; 
                std::cout <<std::setw(8) <<  C[i][j]; 
            }
            std::cout << std::endl;
}