function plot_lorenz(t0,t1,y0,b,s,r)

tspan = [t0, t1];
yzero = y0; 
options = odeset('AbsTol',1e-7,'RelTol',1e-4);

[t,y] = ode45(@lorenzode,tspan,yzero,options,b,s,r);
figure(1)
plot3(y(:,1),y(:,2),y(:,3))

figure(2)
subplot(2,2,1)
plot(y(:,1),y(:,2)) , xlabel('x'), ylabel('y')
subplot(2,2,2)
plot(y(:,2),y(:,3)) , xlabel('y'), ylabel('z');
subplot(2,2,3)
plot(y(:,1),y(:,3)) , xlabel('x'), ylabel('z');
subplot(2,2,4)
plot(t,y(:,1)) , xlabel('t'), ylabel('x');