f=inline('[x(2); -sin(x(1))]','t','x');
options=odeset('OutputFcn','odephas2');
ode45(@(t,x) f(t,x),[0 24],[0 2],options);
hold on;
ode45(@(t,x) f(t,x),[0 7.2],[0 1],options);
ode45(@(t,x) f(t,x),[0 .1],[0 0],options);