實對稱正定矩陣求逆矩陣
根據(jù)線性代數(shù)理論,實對稱正定矩陣顯然有逆矩陣,而且逆矩陣也是對稱矩陣,即aij=aji(i不等于j).以下為書上原程序,但運行后結(jié)果不對(結(jié)果不是對稱矩陣),正確結(jié)果為:
68 -41 -17 10
-41 25 10 -6
-17 10 5 -3
10 -6 -3 2
請將程序稍加修改,運行結(jié)果正確后再答復給我.
/* 試驗未通過 */
/* 正定矩陣求逆 */
#include
#include
/* #include*/
/* 新系數(shù)的計算公式:*/
int GJ(int,double **); /* a11'=1/a11 */
double **TwoArrayAlloc(int,int); /* a1j'=-a1j/a11 j=2,3,...,n */
void TwoArrayFree(double **); /* ai1'=ai1/a11 i=2,3,...,n */
/* aij'=aij-ai1*a1j/a11 i,j=2,3,...,n */
void main()
{
int i,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