求下三角单位矩阵 ,为什么运行乱码
#include "stdio.h"void main()
{
float A[15][15],L[15][15],U[15][15];
int i,j,k;
float s;
A[1][1]=2;
A[1][2]=1;
A[1][3]=1;
A[2][1]=1;
A[2][2]=3;
A[2][3]=2;
A[3][1]=1;
A[3][2]=2;
A[3][3]=2;
for(j=1;j<=3;j++)
U[1][j]=A[1][j];
for(i=2;i<=3;i++)
L[i][1]=A[i][1]/U[1][1];
for(i=2;i<=3;i++)
{
for(j=i;j<=3;j++)
{
s=0;
for(k=1;k<=i-1;k++)
s=s+L[i][k]*U[k][j];
U[i][j]=A[i][j]-s;
}
for(j=i+1;j<=3;j++)
{
s=0;
for(k=1;k<=i-1;k++)
s=s+L[j][k]*U[k][i];
L[j][i]=(A[j][i]-s)/U[i][i];
}
}
for(i=1;i<=3;i++)
{
for(j=1;j<i;j++)
printf("%4.1f",L[i][j]);
printf(" 1\n");
}
for(i=1;i<=3;i++)
{
for(j=1;j<=3;j++)
printf("%4.1f",U[i][j]);
printf("1\n");
}
}