Вычислить определитель матрицы СЛАУ:(10 0 5 20 2 1 25 1 5 ).Вектор правой части СЛАУ заполнен единицами

//Вычислить определитель матрицы СЛАУ:(10 0 5 20 2 1 25 1 5 )Вектор правой части СЛАУ заполнен единицами.
#include <iostream>
#include "matrio.h"
#include "solve.h"
#define n 3
#define m 3
 
using namespace std;
 
int main()
{
    int **a;
    float **b;
    a=new int *[n];
    for (int i=0;i<n;i++)
        a[i]=new int [n];
    b=new float *[n];
    for (int i=0;i<n;i++)
        b[i]=new float [n];
 
    InputFileMatr(a,n,m);
    cout<<endl<<"A: "<<endl;
    OutputDescMatr(a,n,m);
 
    triangMatr(a,b,n);
    cout<<"Triang A:"<<endl;
    OutputDescMatrF(b,n,m);
 
 
 
    for (int i=0;i<n;i++)
        delete[] a[i];
    delete[] a;
    a=NULL;
    for (int i=0;i<n;i++)
        delete[] b[i];
    delete[] b;
    b=NULL;
    return 0;
}
bool InputFileMatr (int** M,int n,int m)
{
    ifstream f("matrix.txt");
 
    if(!f)
    {
        cout<<"Error,there are no file to input matrix"<<endl;
        return false;
    }
 
    for (int i=0;i<n;i++)
        for (int j=0;j<m;j++)
        if(!(f>>M[i][j]))
        {
            cout<<"Error:not enough data or bed data in the file"<<endl;
            f.close();
            return false;
        }
    f.close();
    return true;
}
void OutputDescMatr (int** M,int n,int m)
{
    for (int i=0;i<n;i++)
    {
        for (int j=0;j<m;j++)
            {
                cout.setf(ios::fixed);
                cout.precision(2);
                cout<<setw(5)<<M[i][j];
            }
        cout<<endl;
    }
    cout<<endl;
}
 
void OutputDescMatrF (float** M,int n,int m)
{
    for (int i=0;i<n;i++)
    {
        for (int j=0;j<m;j++)
            {
                cout.setf(ios::fixed);
                cout.precision(2);
                cout<<setw(10)<<M[i][j];
            }
        cout<<endl;
    }
    cout<<endl;
}
void triangMatr (int **A,float **B,int N)
{
    int i,j,k;
    float koef;
    for (i=0;i<N;i++)
        for (j=0;j<N+1;j++)
            B[i][j]=A[i][j];
    for (k=0;k<N;k++)
        {
            for (i=k+1;i<N;i++)
            {
                koef=-1.*B[i][k];
                for (j=k;j<N+1;j++)
                    B[i][j]=B[i][j]+B[k][j]*koef;
            }
        }
}