Octave程序解决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;