c++用前代法求解下三角方程组

#include<stdio.h>
#include<math.h>
#define n 3
main()
{
int i,j;
double m0=0;
double y[n];
double b[n];
double L[n][n];
for(i=0;i<n;i++)
{
printf("输入矩阵元素：\n");
for(j=0;j<n;j++)
{
scanf("%lf",&L[i][j]);
}
printf("b[%d]=",i);
scanf("%lf",&b[i]);
}
for(i=0;i<n;i++)
{
m0=0;
for(j=1;j<i-1;j++)
m0=m0+L[i][j]*y[j];
y[i]=(b[i]-m0)/L[i][i];
printf("结果为 y[%d]=%lf\n",i,y[i]);
}
｝

怎么修改？？