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

[MafoR]

Здравствуйте. Возникла не большая проблемка с задачей по перемножению матриц.
ЗАДАЧА: произведение разреженный матриц (матрицы хранятся во входном файле сначала размер, потом данные).
Вот что у меня получилось:

#include "stdafx.h"
#include <iostream>
#include <fstream>
using namespace std;
int main(){
setlocale(0,"Rus");
int n,m,s,p,x;
ifstream fin;
ofstream fout;
fin.open("dan.txt");
if(fin.fail()){cout<<"Error open"; return 2;}
fin>>n;
fin>>m;
//cout<<n<<endl<<m<<endl;
int **Mas= new int*[n];
for(int i=0; i<n; ++i){
Mas[i] = new int[m];
for(int j=0; j<m; ++j){
fin>>x;
Mas[i][j]=x;
cout.width(3);
cout<<Mas[i][j];
}
cout<<endl;
}
cout<<"Умножается на"<<endl;
fin>>s;
fin>>p;
//cout<<s<<endl<<p<<endl;
int **Mas2= new int*[s];
for(int i=0; i<s; ++i){
Mas2[i] = new int[p];
for(int j=0; j<p; ++j){
fin>>x;
Mas2[i][j]=x;
cout.width(3);
cout<<Mas2[i][j];
}
cout<<endl;
}
cout<<endl;
if(m==s){
int **MasRez= new int*[n];
for(int i=0; i<n; ++i){
MasRez[i] = new int[p];
for(int j=0; j<p; ++j){
for(int r=0; r<s;++r ){
MasRez[i][j]+=Mas[i][r]*Mas2[r][j];
if(r==s-1){
cout.width(3);
cout<<MasRez[i][j];
}
}
cout<<endl;
}
}
}
return 0;
}

Проблема в том, что у меня даже умножение простых матриц не получается. Сделал все как нужно по алгоритму, по всем правилам, а он выводит всякую ахинею. Не могу понять, где допустил оплошность. Помогите, пожалуйста. Заранее при много благодарен.

[MafoR]

Нашел в чем ошибка. Забыл обнулить массив =)
MasRez[i] = new int[p];
for(int j=0; j<p; ++j){
MasRez[i][j]=0;
for(int r=0; r<s;++r ){
MasRez[i][j]+=Mas[i][r]*Mas2[r][j];
подскажите как можно перемножить разреженные матрицы?