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

[Macmeprag]

#include <cstdio>
#include <iostream>

typedef long double** matr;
int p;

matr mult(matr, matr);

int main(){
	freopen("input.txt","r",stdin);
	freopen("output.txt","w",stdout);
	int n; scanf("%d%d", &n, &p);
	matr a = new long double *[p];
	for (int i = 0; i < p; i++) a[i] = new long double [p];
	for (int i = 0; i < p; i++)
		for (int j = 0; j < p; j++) scanf("%lf", &a[i][j]);
	a = mult(a, a);
	return 0;
}

matr mult(matr a, matr b){
	matr subMatr = new long double *[p];
	for (int i = 0; i < p; i++) subMatr[i] = new long double [p];
	for (int i = 0; i < p; i++)
		for (int j = 0; j < p; j++) subMatr[i][j] = 0;	
	for (int i = 0; i < p; i++)
		for (int j = 0; j < p; j++)
			for (int k = 0; k < p; k++)
				subMatr[i][j] += (a[i][k] * b[k][j]);					
	return subMatr;
}

вот код, алгоритм вроде верный, но где-то в умножении матрицы на матрицу все теряется..a и b передаются нормально. Помогите пожалуйста!!

[onewho]

у меня ваш код работает (вроде даже верно перемножил)
кроме того что n считывается и не используется.

	for (int i = 0; i < p; i++) {
		for (int j = 0; j < p; j++) {
			printf("%.3lf ", a[i][j]);
		}
		puts("");
	}

тестовый вывод.

[Macmeprag]

n нужно будет для возведения матрицы в степень, где умножение - подзадача..
у меня что-то не так с компьютером? я уже и на винде и на убунте запускал и там и там в subMatr чушь какая-то, причем разная(

[onewho]

попробуй long double на просто double замени