function stated = MyLorenz(t, state)
x = state(1);
y = state(2);
z = state(3);
% these are our constants
sigma = 10.0; rho = 28.0; beta = 8/3;
% compute the state derivatives
xd = sigma * (y-x);
yd = (rho-z)*x - y;
zd = x*y - beta*z;
% return the state derivatives
 stated = [xd; yd; zd];
end
