doty=@(x,y) -2*y+2*x^2+2*x;  %定义匿名函数
[x,y]=ode45(doty,[0,0.5],1)
