
T=50;

[t,Y]=ode45(@(t,Y) [f(Y(1),Y(2),Y(3),10); g(Y(1),Y(2),Y(3),28); 
    h(Y(1),Y(2),Y(3),8/3)],[0 T],[0,1,0]);
figure(1);
plot3(Y(:,1),Y(:,2),Y(:,3),'b', 'LineWidth',1);
hold on;
figure(2);
plot(t, Y(:,1), 'b','LineWidth',2);
hold on;

[tt,Y]=ode45(@(tt,Y) [f(Y(1),Y(2),Y(3),10); g(Y(1),Y(2),Y(3),28); 
    h(Y(1),Y(2),Y(3),8/3)],[0 T],[0,1.001,0]);

figure(1);
plot3(Y(:,1),Y(:,2),Y(:,3),'r','LineWidth',1);
hold off

figure(2);
plot(tt, Y(:,1),'r','LineWidth',1);
hold off;



function z=f(x,y,z,s)
z=s*(y-x);
end

function z=g(x,y,z,r)
z=r*x-y-x*z;
end

function z=h(x,y,z,b)
z=-b*z+x*y;
end




