%yp = [ ];
clear y1max y2max;
load lorin;
y0 = lorin;
options = odeset('RelTol',.0000001,'AbsTol',.000000001)
%y0 = y0 -.09;
r = 25.8;
tspan = [0:.01:50];
[tt,yy] = ode45('dlor',tspan,y0,options,r);
plot3(yy(:,1),yy(:,2),yy(:,3));
sy = size(yy);
xy = yy(sy(1),:)';
xy = yy(sy(1),:)';
save xyout xy -ascii