function dxdt=lorenz63_dxdt(t,x)
%% 
%    UMD AOSC 615   SPRING 2009
%   dxdt for Lorenz '63 systems 

%% parameters
global ode_prmtr
% ode_prmtr=[a b r];    

%% dxdt
dxdt=zeros(size(x));
dxdt(1)=                   - ode_prmtr(1)*(x(1)-x(2));
dxdt(2)=-x(1)*x(3) + ode_prmtr(3)*x(1)-x(2); 
dxdt(3)=  x(1)*x(2) - ode_prmtr(2)*x(3); 

%% end
return