function yp=sistema(t,y);
yp=zeros(2,1);
yp(1)=y(2);
yp(2)=(32+2*t^3-y(1)*y(2))/8;