Обратная матрица

#include <stdio.h>
#include <stdlib.h>
 
void showinfo();
void readmatrix(double **, int);
void showmatrix(double **, int);
double determinant(double** ,int );
double addings(double ** a,int n,int x,int y);
void inversion(double ** a,double **inv,int n);
void transponation(double **a,int n);
void multiply(double **a,double *b,double *c,int n);
 
int main()
{
    int n;
    showinfo();
    printf("\n");
 
    printf("Input number of equations\n");
    scanf("%d",&n);
    printf("\n");
    double **m=new double*[n];
    double **m1=new double*[n];
    for (int i=0;i<n;i++)
    {
        m[i]=new double[n];
        m1[i]=new double[n];
    }
 
    readmatrix(m,n);
    printf("\n");
    showmatrix(m,n);
    printf("\n");
 
    double *w=new double[n];
    double *q=new double[n];
    
    printf("Input B matrix\n");
    for (int i=0;i<n;i++)
    {
        scanf("%lf",&w[i]);
    }
    printf("\n");
 
    inversion(m,m1,n);
    multiply(m1,w,q,n);
    
    printf("Result is\n");
    for (int i=0;i<n;i++)
    {
        printf("%6.2lf",q[i]);
    }
    system("pause");
    return 0;
}
 
void showinfo()
{
    printf("    | a11   ....   a1n |        |  x1  |       |  b1  |\n");
    printf("A = | ....  ....  .... | ,  x = | .... | , b = | .... |\n");
    printf("    | am1   ....   amn |        |  xn  |       |  bm  |\n");
}
 
void readmatrix(double **a, int l)
{
    printf("Input A-matrix\n");
    for (int i=0;i<l;i++)
    {
        for (int j=0;j<l;j++)
        {
            scanf("%lf",&a[i][j]);
        }
    }
}
 
void showmatrix(double **a,int l)
{
    printf("Your matrix\n");
    for (int i=0;i<l;i++)
    {
        for (int j=0;j<l;j++)
        {
            printf("%6.2lf",a[i][j]);
        }
        printf("\n");
    }
}
 
double determinant(double **a, int n)
{
    int k,l,j,k1,l1,sign,n1=n-1;
    double d1,d;
    if (n==1)
        d1=a[0][0];
    if (n==2)
        d1=a[0][0]*a[1][1]-a[0][1]*a[1][0];
    if (n>2)
    {
        d1=0;
        for (j=0;j<n;j++)
        {
            if (j%2==0) sign=1;
            else sign=-1;
            k=1;l=0;n1=n-1;
            double **temp=new double*[n1];
            for (int i=0;i<n1;i++)
                temp[i]=new double[n1];
            for (k1=0;k1<n1;k1++,k++)
            {
                l=0;
                for (l1=0;l1<n1;l1++,l++)
                {
                    if (l1==j) l++;
                    temp[k1][l1]=a[k][l];
                }
            }
            d=determinant(temp,n1);
            for (int i=0;i<n1;i++)
                delete[] temp[i];
            delete[] temp;
            d1+=a[0][j]*sign*d;
        }
    }
    return d1;
}
 
double addings(double **a,int n,int x,int y)
{
    double res,d;
    int sign,n1,k=0,l=0;
    if (n==1)
        res=a[0][0];
    if (n==2)
        res=a[0][0]*a[1][1]-a[0][1]*a[1][0];
    if (n>2)
    {
        res=0;
        n1=n-1;
        double **temp_2=new double*[n];
        for (int i=0;i<n1;i++)
            temp_2[i]=new double[n1];
        if ((x+y)%2==0) sign=1;
        else sign=-1;
        for (int i=0;i<n;i++)
        {
            if (i==x)
                continue;
            k=0;
            for (int j=0;j<n;j++)
            {
                if (j==x)
                    continue;
                temp_2[l][k]=a[i][j];
                k++;
            }
            l++;
        }
        d=determinant(temp_2,n1);
        res=sign*d;
        for (int i=0;i<n1;i++)
            delete[] temp_2[i];
        delete[] temp_2;
    }
    return res;
}
 
void inversion(double **a,double **inv,int n)
{
    double det=determinant(a,n);
    if (det==0)
    {
        printf("Matrix nevyrodgena\n");
        exit(0);
    }
    for (int i=0;i<n;i++)
    {
        for (int j=0;j<n;j++)
        {
            inv[i][j]=addings(a,n,i,j)/det;
        }
    }
    transponation(inv,n);
}
 
void transponation(double **a,int n)
{
    double vault;
    for (int i=0;i<n;i++)
    {
        for (int j=i;j<n;j++)
        {
            vault=a[i][j];
            a[i][j]=a[j][i];
            a[j][i]=vault;
        }
    }
}
 
void multiply(double **a,double *b,double *q,int n)
{
    for (int i=0;i<n;i++)
    {
        q[i]=0;
    for (int j=0;j<n;j++)
    {
        q[i]+=a[i][j]*b[j];
    }
    }
}