%yp = [ ]; y0 = [.2 .3 .4]; options = odeset('RelTol',.0000001,'AbsTol',.000000001) r = 1.5; tspan = [0:.01:50]; [tt,yy] = ode45('lor',tspan,y0,options,r); plot3(yy(:,1),yy(:,2),yy(:,3));