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

[Artemm]

Люди помогите пожалуйста, вот написал код , а он не пашет.

#include<stdio.h>
#include<math.h>
#include<conio.h>
float a[10][20],x[10][10],E[10][10],c,d[10][10];
int i,j,k,n,h;
main()
{
m2: printf("\n vvedit n \n");
	scanf("\n %d",&n);
	printf("\n n=%d",n);
	printf("\n vvedit matrizu A \n");
	for (i=0;i<n;i++)
    for (j=0;j<n;j++)
		scanf("\n %f",&a[i][j]);
	for(i=0;i<n;i++)
	{
		a[i][i+n]=1.0;
		printf("\n");
		for(j=0;j<n;j++)
			printf("%f",a[i][j]);
	}
	for(i=0;i<n;i++)
	for (j=0;j<n;j++)
	{
		d[i][j]=a[i][j];
	}
	for(j=0;j<n;j++)
	{
		for(i=j+1;i<n;i++)
		{
			c=a[i][j]/a[i][j];
			for(k=0;k<2*n;k++)
			{
				a[i][k]=(c*a[i][k])-a[j][k];
			}
		}
	}
	c=1;
	for(i=0;i<n;i++)c=c*a[i][j];
	if(c==0)
	{
		printf("\n Det=0");
	}
	for(j=0;j<n;j++)
	{
		for(i=n-1;i>=0;i--)
		{
			x[i][j]=a[i][j+n];
			for(k=n-1;k>i;k--)
			{
				x[i][j]=x[i][j]-(a[i][k]*x[k][j]);
			}
			x[i][j]=x[i][j]/a[i][j];
		}
	}
	printf("\n Obernena matriza= ");
	for(i=0;i<n;i++)
	{
		printf(" \n ");
		for(j=0;j<n;j++)
			printf(" %f", x[i][j]);
	}
	for(i=0;i<n;i++)
		for(j=0;j<n;j++)
		{
			E[i][j]=0;
			for(h=0;h<n;h++)
				E[i][j]=E[i][j]+(d[i][j]*x[i][j]);
		}
		printf("\n E=a*x/n");
		for(i=0;i<n;i++)
		{
			printf("\n");
			for(j=0;j<n;j++)
				printf("%f",E[i][j]);
		}
	goto m2;
	return 0;
	getch ();


    
}

[Levsha100]

Что конкретно за ошибка?. И еще дайте само задание

[Artemm]

Задание таково.
Для заданной матрицы найти оборотную и сделать проверку.

1.6 1.6 1.7 1.8
1.6 2.6 1.3 1.3
1.7 1.5 3.6 1.4
1.8 1.3 1.4 4.6


ошибка в том что оно не правильно считает .

[Sazary]

Переписал в соответствии с этим алгоритмом:
http://algolib.narod.ru/Math/ObrMatrix.gif

#include <stdio.h>
#include <conio.h>

int main()
{
float M[100][100],obr[100][100];
int i,j,k,N;
printf("N = ");
scanf("%d",&N);
for(i=0; i<N; i++)
 for(j=0; j<N; j++)
  {
   printf("N[%d][%d] = ",i,j);
   scanf("%f",&M[i][j]);
  }
printf("----------------\n");
printf("Matrix: \n");
for(i=0; i<N; i++)
 {
  for(j=0; j<N; j++)
   printf("%5.1f",M[i][j]);
  printf("\n");
 }

for(k=0; k<N; k++)
{
  for(i=0; i<N; i++)
   for(j=0; j<N; j++)
    {
      if (i==k&&j==k)
       obr[i][j] = 1/M[i][j];
      if(i==k&&j!=k)
       obr[i][j] = -M[i][j]/M[k][k];
      if (i!=k&&j==k)
       obr[i][j] = M[i][k]/M[k][k];
      if(i!=k&&j!=k)
       obr[i][j] = M[i][j] - M[k][j] * M[i][k]/M[k][k];
     }
 for(i=0;i<N;i++)
  for(j=0; j<N;j++) M[i][j]= obr[i][j];
}

printf("\nObr matrix: \n");
for(i=0; i<N; i++)
 {
  for(j=0; j<N; j++)
   printf("%5.1f",M[i][j]);
  printf("\n");
 }

getch ();
return 0;
}

[Artemm]

Супер Спасибо , а проверку можно еще сюда прилепить ?

[Sazary]

Сообщение от Artemm

а проверку можно еще сюда прилепить ?

Что за проверка?
Это вы уж как-нибудь сами.

[Artemm]

Ок , спасибо.