Question

c语言 追赶法

#include<stdio.h>
#define N 3
main()
{
float a[N-1]={1.0,2.0};
float b[N]={1.0,2.0,3.0};
float c[N-1]={1.0,1.0};
float d[N]={2.0,4.0,5.0};
int i;
float x[N],y[N],r[N];
y[0]=c[0]/b[0];
r[0]=d[0]/b[0];
for(i=1;i<N-1;i++)
{
y[i]=(d[i]-a[i-1]*y[i-1])/(b[i]-a[i-1]*r[i-1]);
r[i]=c[i]/(b[i]-a[i-1]*r[i-1]);
} printf("%f",y[1]);
y[N-1]=(d[N-1]-a[N-2]*y[N-2])/(b[N-1]-a[N-2]*r[N-2]);
x[N-1]=y[N-1];
for(i=N-2;i>=0;i--)
{
x[i]=y[i]-r[i]*x[i+1];
}
printf("x=(");
for(i=0;i<N;i++)
{
printf(" %f ",x[i]);
}
printf(")\n");
}
这一段程序哪里错了

Answer 1

b[1]=2
a[0]=1
r[0]=d[0]/b[0]=2/1=2

i等于1时，b[i]-a[i-1]*r[i-1] = 0
分母为0了！