Question

C语言顺序高斯消元法求上三角，不知哪里错了，请各位大哥大姐帮帮我吧！

#include<stdio.h>
#include<math.h>
void main(void)
{
long i,j,k,a[3][3];
printf(" -^-^-^-*-*------please input the a[3][3].------*-*-^-^-^-\n");
for(i=0;i<3;i++)
for(j=0;j<3;j++)
scanf("%d",&a[i][j]);
for(k=0;k<2;k++)
for(i=k+1;i<3;i++)
{
a[i][k]=-a[i][k]/a[k][k];
for(j=k+1;j<3;j++)
a[i][j]+=a[i][k]*a[k][j];
}
for(i=0;i<3;i++){
for(j=0;j<3;j++)
printf("%5d ",a[i][j]);
printf("\n");}
printf("shun xu gao si xiao yuan fa jie guo:\n");
for(i=0;i<3;i++){
for(j=0;j<3;j++)
{if(i>j)
{a[i][j]=0;}
else{a[i][j]=a[i][j];}
printf("%5d ",a[i][j]);}
printf("\n");}
getchar();
getchar();
getchar();
}

Answer 1

for(k=0;k<2;k++)
for(i=k+1;i<3;i++)
{if (a[k][k]!=0)
for(j=k;j<3;j++) a[i][j]=(a[i][j]*a[k][k]-a[k][j]*a[i][k])/a[k][k];
}