常微分方程的初值问题一般可以ode45（）函数命令求解，其计算精度比其他ode（）函数要高。
实现代码：
y0=[1,0];
[x,y]=ode45(@odefun,[0,20],y0)
自定义函数
function [dydx]= odefun(x,y)
dydx=zeros(2,1);
dydx(1)=y(2);
dydx(2)=1+3*y(2)-2*y(1);
end
运行结果
