dY=inline('[y(2);(y(1)*y(4)^2+9.8*(cos(y(3)-1)))/2;y(4);-(9.8*sin(y(3))+2*y(2)*y(4))/y(1)]','t','y');
ode45(dY,1,[.5 0 0 .5])
