我正在尝试在R中绘制ODE的解决方案,并输入:
dx=0.01
x=seq(-5,5,by=dx)
A=0.5
B=1
C=-2
length(Y)=length(x)-2
Y=x*0
Y[1]=A
Y[2]=dx*B+A
for(k in 1:(length(x)-3))(C*dx^2-1)*Y[k]+2*Y[k+1]
plot(x[1:(length(x)-2)],Y,type="l",col="blue")
它显示了错误:
Error in xy.coords(x, y, xlabel, ylabel, log) : 'x' and 'y' lengths differ