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

[BDA]

Добавил ввод. Кое-что чуть поменял. (100% работоспособность не гарантирую)

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

#define E 0.000001

void
mover(double *a,int to,int from,int n)
{
    double tmp[2*n];
    int i;
    memcpy(tmp,(a+2*from*n),2*sizeof(double)*n);
    memcpy((a+2*from*n),(a+2*to*n),2*sizeof(double)*n);
    memcpy((a+2*to*n),tmp,2*sizeof(double)*n);
}

void
adder(double *a,int to,int from,int n)
{
    double d=*(a+2*to*n+from)/(*(a+2*from*n+from));
    int i;
    for (i=0; i<2*n; i++)
        *(a+2*to*n+i)-=(*(a+2*from*n+i))*d;
}

void
printmatrix(double *a, int n)
{
    int i,j;
    for(i=0; i<n; i++)
    {
        for (j=0; j<n; j++)
            printf("%3.3lf  ",*(a+2*n*i+j));
        putchar('\n');
    }
}

void
printinvmatrix(double *a, int n)
{
    int i,j;
    for(i=0; i<n; i++)
    {
        for (j=0; j<n; j++)
            printf("%3.3lf  ",*(a+2*n*i+j+n));
        putchar('\n');
    }
}

int
main(void)
{
    int n,i,j,q;
    printf("Input size of matrix: ");
    scanf("%d",&n);

    double a[n][2*n];
    memset(a,0,2*n*n*sizeof(double));

    printf("Input elements of matrix:\n");
    for(i=0; i<n; i++)
    {
        for (j=0; j<n; j++)
            scanf("%lf",&(a[i][j]));
    }

    for (i=0; i<n; i++)
        a[i][n+i]=1;

    printf("Source matrix:\n");
    printmatrix(a,n);

    for (i=0; i<n; i++)
    {
        j=i;
        while((fabs(a[j][i])<E)&&(j<n))
            j++;
        if (j==n)
        {
            printf("Matrix does not have the inverse matrix");
            return 0;
        }
        if (j!=i)
            mover(a,i,j,n);
        for (q=i+1; q<n; q++)
            adder(a,q,i,n);
    }

    for (i=n-1; i>=0; i--)
    {
        for (q=i-1; q>=0; q--)
            adder(a,q,i,n);
        for (q=2*n-1; q>=i; q--)
            a[i][q]/=a[i][i];
    }

    for (i=0; i<n; i++)
        for (j=n; j<2*n; j++)
            if (fabs(a[i][j])<E)
                a[i][j]=0;

    printf("\nThe inverse matrix:\n");
    printinvmatrix(a,n);
    return 0;
}

УПД Еще чуток оптимизирована.

[nat@]

Спасибо .. сама бы я этого не написала.))