Матрицы не хотят умножаться (внутри практически рабочий код!)

#include <iostream>
#include <math.h>
#include <fstream>
using namespace std;
 
double determinant(double **matrica_a, int n)
{ int i,j,k,r; double c,M,max,det=1, **a;
a=new double *[n];
for(i=0;i<n;i++) a[i]=new double[n];
for(i=0;i<n;i++)
 for(j=0;j<n;j++)
a[i][j]=matrica_a[i][j];
for(k=0;k<n;k++)
{
max=fabs(a[k][k]); r=k;
for(i=k+1;i<n;i++)
if (fabs(a[i][k])>max)
{ max=fabs(a[i][k]); r=i; }
if (r!=k) det=-det;
for(j=0;j<n;j++)
{c=a[k][j]; a[k][j]=a[r][j]; a[r][j]=c; }
for(i=k+1;i<n;i++)
for(M=a[i][k]/a[k][k],j=k;j<n;j++)
a[i][j]-=M*a[k][j];
}
for(i=0;i<n;i++)
det*=a[i][i];
return det;
}
 
int main()
{
    int n,i,j,max=0;
    cout<<"Ввод количества элементов массива \n";
    cin>>n;
    cout<<endl;
    double q;
    double **a;
    double **b;
    int *c=new int [n];
    a=new double *[n];
    for (i=0;i<n;i++)
    a[i]=new double (n);
    fstream fi;
    fi.open("/home/vitalya/abcd.txt", ios::in);
    if (!fi)
    cout<< "ошибка!\n"; 
    for(i=0;i<n;i++)
    fi>>c[i];
    cout<<"Чудо-массив \n";
    for(i=0;i<n;i++)
    cout<<c[i]<<' ';
    cout<<endl;
    fi.close();
    cout<<"Транспонированный массив \n";
    for (i=0;i<n;i++)
    cout<<c[i]<<endl;
    cout<<endl;
    for (i=0;i<n;i++)
    for (j=0;j<n;j++)
    a[i][j]=c[i]*c[j];
    cout<<"Произведение (Матрица A[i][j]) \n";
    for (i=0;i<n;cout<<endl,i++)
    for (j=0;j<n;cout<<"\t",j++)
    cout<<a[i][j];
    cout<<endl;
    max=a[0][0];
    for (i=0;i<n;i++)
    for (j=0;j<n;j++)
    if (max<a[i][j])
    max=a[i][j];
    b=new double *[n];
    for (i=0;i<n;i++)
    b[i]=new double [n];
    for (i=0;i<n;i++)
    for (j=0;j<n;j++)
    b[i][j]=max/a[i][j];
    cout<<"Матрица, состоящая из A(i,j)/max|A| (Матрица B[i][j]) \n";
    for (i=0;i<n;cout<<endl,i++)
    for (j=0;j<n;j++)
    cout<<b[i][j]<<"\t";
    cout<<endl;
    for (i=0;i<n;i++)
    for (j=0;j<n;j++)
    if (i==j)
    b[i][j]=3*b[i][j]-1;
    else 
    b[i][j]=3*b[i][j];
    cout<<"3*В-Е \n";
    for (i=0;i<n;cout<<endl,i++)
    for (j=0;j<n;cout<<"\t",j++)
    cout<<b[i][j];
    cout<<endl;
    q=determinant(a,n);
    cout<<"Определитель равен "<<q<<", то есть нулю";
}