Умножение матриц

#include <iostream>
#include <cstdlib>  
#include <stdio.h>
#include <conio.h>  
 
using namespace std; 
void PrintMas(double *A,int n){
 
     for (int i = 0 ; i < n ; i++){
         for (int j = 0 ; j < n ; j++){
             cout<<*(A+i*n+j)<<"\t ";
         }
         cout<<endl;
     }
     cout<<endl;
     return ;
}
void EnterMas(double *A,int n){
     for (int i = 0 ; i < n ; i++){
         for (int j = 0 ; j < n ; j++){           
             cout<<"Введите ["<<j+1<<"] элемент ["<<i+1<<"] столбца = ";
             cin>>*(A+i*n+j);           
         }
     }
     return ;
}
void FillBMas(double *B,int n){ 
     for (int i = 0 ; i < n ; i++){
         for (int j = 0 ; j < n ; j++)
         {              
            *(B+i*n+j)=1/double((i)+j-1); 
             
         }
     }
     return ;
}
void MnojMartixMas(double *A,double *B, double *C,int n){
// 
    for(int irow = 0; irow < n; irow++) 
        for(int icol = 0; icol < n; icol++) {
            double sum = 0;
            for(int k = 0; k < n; k++) {
                sum += (*(A+irow*n+k)) * (*(B+k*n+icol));
            }
            *(C+irow*n+icol) = sum;
        }
     return ;
}
main ()
{
  char a;
  do
    {
     setlocale(LC_ALL,"Rus");
     int n; 
     do{
          cout<<Матрица\n Введите размерность матрицы N = ";
          cin>>n;
          if (n<=0)  cout<<Ощибка : N > 0 !\n";
     }while (n<=0); 
 
     double *A = new double [n*n];
     double *B = new double [n*n]; 
     double *C = new double [n*n]; 
 
     EnterMas(A,n); 
     PrintMas(A,n); 
     FillBMas(B,n); 
     PrintMas(B,n); 
     MnojMartixMas(A,B,C, n);  
     PrintMas(C,n); 
     delete(A);  delete(B);  delete(C);   
     cout<<"Хотите продолжить? (y/n)";
     cin >>a;
     
     }
     while (a=='y');
     system ("PAUSE");
     return 0;
}