R误差x和y长度不同

Question

我正在尝试在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

Answer 1

需要切换第6行和第7行。在创建变量之前，不能更改变量的长度。另外，我认为你的意思是将for循环输出分配给变量。

dx=0.01
x=seq(-5,5,by=dx)
A=0.5
B=1
C=-2
Y=x*0
length(Y)=length(x)-2
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")