[MWorks Syslab科学计算] 微分方程数值解ode45

using TyBase
using TySymbolicMath

using TyMath
using TyPlot

 function odefun(t,y)
     dydt = [y[2]; -t*y[1]+exp(t)*y[2]+3*sin(2*t)]
     return dydt
 end
# function odefun(t,y)
#     f=[2*t];
#     return f
# end

tspan = [1 2];
y0 = [0; 0];
t,y, = ode45(odefun, tspan, y0);

plot(t,y[:,1],"-o",t,y[:,2],"-o")