function [x,y,z] = lorenz_att1(x0,y0,z0,h,steps)
% x` = sigma(y - x)
% y` = x(ro - z) - y
% z` = xy - beta*z

        beta = 8/3;
        sigma = 10;
        ro = 28;

[x y z] = rk4(x0,y0,z0,beta,sigma,ro,h,steps);

