function dx = Lorenz(t, x)
a = 10; b = 28; c = 8/3;
dx = zeros(1,3);
dx(1) = a*(x(2) - x(1));
dx(2) = x(1)*(b - x(3)) - x(2);
dx(3) = x(1)*x(2) - c*x(3);
function dx = Lorenz(t, x)
a = 10; b = 28; c = 8/3;
dx = zeros(1,3);
dx(1) = a*(x(2) - x(1));
dx(2) = x(1)*(b - x(3)) - x(2);
dx(3) = x(1)*x(2) - c*x(3);