求解 ODE 的 Octave 程序

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;