大连理工大学矩阵与数值分析课又叫计算机科学计算的上机实验报告,内容详尽,有程序代码和预算结果

}

for(k=1;k<=n-1;k++) {

ark=max(k);

if(ark==0)

{

printf("\n\nIt’s wrong!");message();

}

else if(flag!=k)

exchange(flag,k);

for(i=k+1;i<=n;i++)

for(j=k+1;j<=n+1;j++)

A[i][j]=A[i][j]-A[k][j]*A[i][k]/A[k][k];

}

x[n]=A[n][n+1]/A[n][n];

for( k=n-1;k>=1;k--)

{

float me=0;

for(j=k+1;j<=n;j++)

{

me=me+A[k][j]*x[j];

}

x[k]=(A[k][n+1]-me)/A[k][k];

}

for(i=1;i<=n;i++)

{

printf(" \n\nx%d=%f",i,x[i]);

}

message();

return 0;

}

void exchange(int r,int k)

{

int i;

for(i=1;i<=n+1;i++)

A[0][i]=A[r][i];

for(i=1;i<=n+1;i++)

A[r][i]=A[k][i];

for(i=1;i<=n+1;i++)

A[k][i]=A[0][i];

}

float max(int k)