%Script file to drive the Lorenz solver:
y0=[3,3,6]; y0=y0';
tspan=0:0.005:50;

%options=odeset('OutputFcn','odeplot');
%options=odeset('OutputFcn','odephas3');

%[t,y]=ode23('lorenzeq',tspan,y0,options);

[t,y]=ode23s('lorenzeq',tspan,y0);

%Plot an individual time series (here, the second coordinate)
plot(t,y(:,2));