function dydt = lorenzField(t, Y, options, flag, sigma, rho, beta)

x = Y(1);
y = Y(2);
z = Y(3);

dydt = [sigma*(y-x); 
        x*(rho - z) - y;
        x*y - beta*z];