% by Dr. Ute Feldmann (14.11.11)
s=10;b=8/3;r=28;
y0=[18,18,27];
options = odeset('OutputFcn',@odephas3)
%options = odeset('OutputFcn',@odeplot)
[t,l]=ode23(@lorenzD,[0 100],y0,options,s,b,r);
plot3(l(:,1),l(:,2),l(:,3));
% use view-->figure toolbar-->rotate3d and have fun turning the butterfly