-1
我試圖繪製解決在R A ODE,我打字:半徑誤差x和y的長度不同
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
是,'長度(X [1:(長度(X)-2)])''給出和999''長度(Y)''給出1001'。 –
for循環的要點是什麼?它不分配任何東西。你的意思是'Y [k] < - (C * dx^2 -....'? –
這個情節適用於我,但是你的調用順序是錯誤的,'length(Y)'在改變之前被改變 –