Умножение + создание матриц

#include <stdlib.h>
#include <conio.h>
#include <iostream.h>
#include <stdio.h>
//---------------------------------------------------------------------------
enum {max_size=10};
struct Matrix {
        double coef[max_size][max_size];
        int sx, sy;
        };
 
void CreatMatrix(Matrix &m, int a, int b);
void PrintMatrix(const Matrix &);
 
void SummMatrix(
        const Matrix & m1,
        const Matrix & m2,
                Matrix & res);
void MultMatrix(
                const Matrix & m1,
                const Matrix & m2,
                Matrix & mres);
 
//---------------------------------------------------------------------------
void main() {
  clrscr();
  int n,qn,a,b,d,c;
  double tmp;
  cout <<"Enter size matrix A\n";
  cin >>a >> b;
  cout <<"Enter size matrix B\n";
  cin>> c >> d;
  clrscr();
  Matrix q1,q2,qSum,qMult;
/*  q1.sx=a; q1.sy=b;
  q2.sx=c; q2.sy=d;
   */
  cout << "Press Matrix A" << "\n";
  CreatMatrix(q1,a,b);
  cout << "Press matrix B" << "\n";
  CreatMatrix(q2,c,d);
  cout <<"\nMatrix A\n";
  PrintMatrix(q1);
  cout <<"\nMatrix B\n";
  PrintMatrix(q2);
  cout <<"\n\tMenu!\n 1.Sum of Matrix;\n 2.Multiplication of Matrix;\n 3.Exit\n";
  cout <<"=>: ";
  cin >>n;
  switch(n) {
    case 1: if (q1.sx!=q2.sx && q1.sy!=q2.sy) {cout <<"Sum - impossible, select 4!\n";
                         break;}
      SummMatrix(q1,q2,qSum);
      PrintMatrix(qSum);
      break;
  case 2: if (q1.sx!=q2.sy) {cout <<"Multiplication - impossible, select 4!\n";
                break;}
      MultMatrix(q1,q2,qMult);
      PrintMatrix(qMult);
      break;
    }
    getch();
}
 
void CreatMatrix(Matrix &m, int a, int b) {
 
m.sx=a; m.sy=b;//Creating matrix1...
for (int i=0; i<m.sy;i++)
for(int j=0;j<m.sx;j++) scanf("%d",&m.coef[i][j]);
}
 
void PrintMatrix(const Matrix &w)
{
for (int i=0;i<w.sy;i++)
    {
for (int j=0;j<w.sx;j++) printf("%d\t",w.coef[i][j]);
cout <<'\n';
     }
}
void SummMatrix( const Matrix &m1,
         const Matrix &m2,
         Matrix &sres) {
  sres.sx = m1.sx; sres.sy = m2.sy;
  cout<<"\nSummMatrix...\n";
  for (int i=0;i<m1.sy;i++)
    for (int j=0;j<m1.sx;j++) sres.coef[i][j]=m1.coef[i][j]+m2.coef[i][j];
}
 
void MultMatrix( const Matrix & m1,
         const Matrix & m2,
         Matrix & mres) {
  mres.sx=m1.sx; mres.sy = m2.sy;
  double M=0;
  cout <<"\nMultMatrix...\n";
  for (int i1=0;i1<m1.sx;i1++)
    for (int j1=0;j1<m2.sy;j1++)
    {
     for(int i2=0;i2<m2.sx;i2++)
     M+=(m1.coef[i1][i2]*m2.coef[i2][j1]);
     mres.coef[i1][j1]=M;
     //      M=0;
    }
 
}

#include <iostream>
#include <set>
#include <vector>
#include <limits>
#include <queue>
#include <string>
#include <map>
#include <algorithm>
#include "Windows.h"
#include <conio.h>
 
using namespace std;
 
const int max_size = 100;
 
struct Matrix {
    double coef[max_size][max_size];
    int n, m;
    Matrix(int _n, int _m) {
        n = _n;
        m = _m;
        for (int i = 0; i < n; i++) 
            for (int j = 0; j < m; j++)
                coef[i][j] = 0;
    }
    void ScanMatrix(int _n, int _m) {
        n = _n;
        m = _m;
        for (int i = 0; i < n; i++) 
            for (int j = 0; j < m; j++)
                scanf("%lf", &coef[i][j]);
    }
    void PrintMatrix() {
        for (int i = 0; i < n; i++) {
            for (int j = 0; j < m; j++) {
                printf("%.15lf ", coef[i][j]);
            }
            printf("\n");
        }
    }
};
 
Matrix MultMatrix(const Matrix & m1, const Matrix & m2) {
    Matrix ans(m1.n, m2.n);
    for (int i = 0; i < m1.n; i++) 
        for (int j = 0; j < m2.m; j++)
            for (int k = 0; k < m2.n; k++)
                ans.coef[i][j] += m1.coef[i][k] * m2.coef[k][j];
    return ans;
}
 
 
 
int main(){            
    freopen("input.txt", "rt", stdin);
    freopen("output.txt", "wt", stdout);
    Matrix a(100, 100), b(100, 100);
    a.ScanMatrix(3, 3);
    b.ScanMatrix(3, 3);
    Matrix c = MultMatrix(a, b);
    c.PrintMatrix();
    return 0;
}