实对称正定矩阵求逆矩阵

根据线性代数理论,实对称正定矩阵显然有逆矩阵,而且逆矩阵也是对称矩阵,即aij=aji(i不等于j).以下为书上原程序,但运行后结果不对(结果不是对称矩阵),正确结果为:

68-41-1710

-412510-6

-17105-3

10-6-32

请将程序稍加修改,运行结果正确后再答复给我.

/*试验未通过*/

/*正定矩阵求逆*/

#include

#include

/*#include*/

/*新系数的计算公式:*/

intGJ(int,double**);/*a11'=1/a11*/

double**TwoArrayAlloc(int,int);/*a1j'=-a1j/a11j=2,3,...,n*/

voidTwoArrayFree(double**);/*ai1'=ai1/a11i=2,3,...,n*/

/*aij'=aij-ai1*a1j/a11i,j=2,3,...,n*/

voidmain()

{

inti,j,n;

double**a;

n=4;

a=TwoArrayAlloc(n,n);

a[0][0]=5;a[0][1]=7;a[0][2]=6;a[0][3]=5;

a[1][0]=7;a[1][1]=10;a[1][2]=8;a[1][3]=7;

a[2][0]=6;a[2][1]=8;a[2][2]=10;a[2][3]=9;

a[3][0]=5;a[3][1]=7;a[3][2]=9;a[3][3]=10;

if(!GJ(n,a))

{

printf("矩阵求逆失败n");

exit(1);

}

printf("该矩阵的逆为：n");

for(i=0;i