Определитель матрицы

#include <cstdlib>
#include <iostream>
#include <conio.h>
#include <stdio.h>
#include <iomanip>
using namespace std;
const int str=10,stolb=10;
int k,l,p;
void input_matrix(double A[][stolb],int m, int n, char name) // Функция ввода
{ 
     cout<<"Input elements of "<<name<<endl;
     for (int i=0;i<m;i++)
     { 
         for(int j=0;j<n;j++)
         { 
                 cin>>A[i][j];
         }
     }
}
void output(double result[][stolb],int m,int n, char name) // функция вывода
{ cout<<"Matrix:"<<name<<endl;
     for(int i=0;i<m;i++)
         {for (int j=0;j<n;j++)
              { 
                   cout<<"  "<<result[i][j];   
              }
              cout<<endl;
         }
     cout<<"____________"<<endl;
}
void multiply(double A[][stolb],double B[][stolb],double C[][stolb]) // функция умножения
{     for(int i=0;i<k;i++)
     {
             for (int j=0;j<p;j++) //J должно быть меньше P, а не L
             {
                 C[i][j]=0;
                 for(int t=0;t<l;t++)
                 {
                         C[i][j] +=A[i][t]*B[t][j];
                 }
             }
     }
}
void trans(double C[][stolb],double D[][stolb]) //Транспонирование
   
{     for(int i=0;i<k;i++)
     {
             for (int j=0;j<p;j++)
             {
 
                         D[i][j]=C[j][i];
             }
     } 
}    
 
int main()
{
START: cout<<" Input number lines and column matrix's A"<<endl; // Вводим количество строк и столбцов матрицы А
cin>>k;
cin>>l;
cout<<"Input number column matrix's B"<<endl; // Вводим кол-во столбцов матрицы B
cin>>p;
cout<<endl;
 double A[str][stolb];
 double B[str][stolb];
 double C[str][stolb];
 double D[str][stolb];
 input_matrix(A,k,l,'A');
 input_matrix(B,l,p,'B');
 
 output(A,k,l,'A');
 output(B,l,p,'B'); 
 
 multiply(A,B,C);
 output(C,k,p,'C');
 
trans(C,D);
output(D,k,p,'D');
 
    
  system("PAUSE");
    goto START; // повторяем цикл
  return 0;
      }

#include <stdio.h>
#include <conio.h>
 
void printMatrix(float** B, int n)  //print the matrix
{
    printf("\n");
    int x ,y;
    for(y=0; y<n; y++)
    {
        for(x=0; x<n; x++)
            printf("%7.2f", B[x][y]);
 
        printf("\n");
    }
    printf("\n");
}
 
//determinant razlozheniem po 1-oj stroke
float det(float **p, int n)
{
    int i,h,x,y;
    if(n==1)
        return **p;
  
    float d=0;
    for(i=0; i<n; ++i)
    {
        float **add = new float*[n-1];
        for(h=0; h<n-1; ++h)
            add[h] = new float[n-1];
    
        for(y=1; y<n; y++)
            for(x=0; x<n; x++)
            {
                if(x==i) continue;
                if(x<i)
                    add[x][y-1] = p[x][y];
                else
                    add[x-1][y-1] = p[x][y];
            }
 
        if(i%2) p[i][0] = -p[i][0];
        d += p[i][0]*det(add, n-1);
 
        for(h=0; h<n-1; ++h)
            delete[] add[h];
        delete[] add;
    }
 
    return d;
}
 
int main()
{
    int n=4,
        x, y;
  
    float **array = new float*[n];
    for(x=0; x<n; x++) 
        array[x] = new float[n];
 
    for(y=0; y<n; ++y)
        for(x=0; x<n; ++x)
            scanf("%f",  &array[x][y]);
    
    printMatrix(array, n);
    printf("determinant = %f\n", det(array, n));
    getch();
 
    for(x=0; x<n; ++x)
        delete[] array[x];
    delete[] array;
 
    return 0;
}