八度程序解決ODE

Question

我有一個ODE

dy=x dx,  y(0)=2

該方程的解為y =x^2/2 + K.

K become 2.

現在我要繪制圖形。 什么時候

 x=0, y=2
 x=1, y=2.5
 x=2, y=4
 x=3, y=6.5
 x=4, y=10

我必須編寫一個Octave program來生成這些值

我的代碼是test.m

function xdot = f (x,t)
 xdot=x;
 endfunction
 x=lsode("f",2,(t=linspace(0,4,5)));
 #plot(t,x)
 x

我在cmd中運行pgm但它給

2.0000
5.4366
14.7781
40.1711
109.1963

預期的結果是

2
2.5
4
6.5
10

請幫我..

Answer 1

終於得到答案了。

采用

xdot=t;

代替

xdot=x;