微分方程，求matlab 代码，不胜感激！求t=[0,20],对应的速度V

微分方程，求matlab 代码，不胜感激！求t=[0,20],对应的速度V

用ode23（）函数可以较好求得微分方程的数值解。
v0=[2.3,0];
h=0.5;
a=0;
b=20;
[t,v]=ode23(@fun,a:h:b,v0);
求解结果：
         t       V(t)       dV(t)
         0    2.3000         0
0.5000    5.2562    3.9856
1.0000    4.4588   -7.0559
1.5000    2.8150    6.5266
2.0000    5.5829    0.7497
2.5000    3.3187   -7.8804
3.0000    3.9351    7.8954
3.5000    5.4400   -2.5105
4.0000    2.4192   -3.2536
4.5000    4.9425    5.3320
5.0000    4.7890   -5.9263
5.5000    2.5304    4.4832
6.0000    5.4779    1.8393
6.5000    3.6658   -7.9648
7.0000    3.5824    7.9280
7.5000    5.4794   -1.6476
8.0000    2.5728   -4.7430
8.5000    4.7258    5.9901
9.0000    4.9100   -5.2930
9.5000    2.4462    3.2745
10.0000    5.3931    2.3321
10.5000    3.7920   -7.8113
11.0000    3.4505    7.7593
11.5000    5.4600   -1.3779
12.0000    2.6260   -5.0057
12.5000    4.6576    6.0813
13.0000    4.8927   -5.2044
13.5000    2.4588    3.1977
14.0000    5.3626    2.2714
14.5000    3.7309   -7.7293
15.0000    3.5055    7.6847
15.5000    5.4053   -1.6398
16.0000    2.5687   -4.3342
16.5000    4.7281    5.6945
17.0000    4.7518   -5.6085
17.5000    2.5561    4.1328
18.0000    5.3754    1.6801
18.5000    3.4960   -7.5912
19.0000    3.7339    7.5863
19.5000    5.2974   -2.4283
20.0000    2.4512   -2.6040
追问您好，初始条件只有速度V=2.3，没有加速度，刚刚搞错了。
能否麻烦你把定义函数的代码也一并说下，所有代码，谢谢追答用ode23（）函数可以较好求得微分方程的数值解。执行代码如下：
v0=0
h=0.5;a=0;b=20;
[t,v]=ode23(@fun,a:h:b,v0);
求解结果：
         t       V(t)
   0                         0
0.5          3.63716291516561
  1          3.64158930086147
1.5           3.6428487523781
  2          3.63885309188913
2.5           3.6339010775317
  3           3.6342955704192
3.5          3.63569824585008
  4          3.63869545971492
4.5          3.64207939462941
  5          3.64279555495064
5.5          3.64076897493905
   6           3.6375611660217
 6.5          3.63456757976042
   7          3.63417276743521
 7.5          3.63640542725326
  8          3.63970151119111
 8.5          3.64250253471956
  9          3.64247879596381
 9.5          3.64003118285705
10          3.63670450855706
10.5          3.63413885787979
11          3.63455245488448
11.5          3.63719178528953
。。。。。。。
18.5          3.63397281939756
  19          3.63569830040656
19.5          3.63887628526499
20          3.64241702048224

追问请问这个里面定义函数那个代码是怎么样的，谢谢

ode45可以追问能否给个代码，谢谢