T=100.0;
x0=0.1;
y0=-2.05;
[t,z]=ode45('pend',[0 T],[x0 y0 0 0]);
plot(z(:,1),z(:,2));