function lorentz(x0,y0,z0)
    xyz=@(t,y) [-10*y(1)+10*y(2);-y(1)*y(3)+28*y(1)-y(2);y(1)*y(2)-8/3*y(3)];
    span=[0,50];
    t0=[x0;y0;z0];
    [t,y]=ode45(xyz,span,t0);
    plot3(y(:,1),y(:,2),y(:,3));
    view(3);
%    axis('off');
    grid('on');
    title(['X0=',num2str(x0),' Y0=',num2str(y0),' Z0=',num2str(z0)]);
    xlabel('x');ylabel('y');zlabel('z');
    axis([-30,30,-30,30,0,50]);