Riga 4: Riga 4:
 
   % dy/dt = x - y - (sqrt(27)+x)z
 
   % dy/dt = x - y - (sqrt(27)+x)z
 
   % dz/dt = sqrt(27)*(x+y) - z + xy
 
   % dz/dt = sqrt(27)*(x+y) - z + xy
 
 
   function dy = lorenz(t,x)
 
   function dy = lorenz(t,x)
 
     dy = zeros(3,1);
 
     dy = zeros(3,1);

Revisione 16:56, 6 Feb 2006

 %Equazione differenziale Ordinaria dell'attrattore di lorenz
 %-----------------------------------------------------------
 % dx/dt = -10 x - 10 y
 % dy/dt = x - y - (sqrt(27)+x)z
 % dz/dt = sqrt(27)*(x+y) - z + xy
 function dy = lorenz(t,x)
    dy = zeros(3,1);
    dy(1) = -10*x(1) + 10 * x(2);
    dy(2) = x(1) - x(2) -(sqrt(27)+x(1))*x(3);
    dy(3) = sqrt(27)*(x(1)+x(2))-x(3)+x(1)*x(2);
 end
 %Equazione differenziale Ordinaria dell'attrattore di lorenz
 %-----------------------------------------------------------
 % dx/dt = -10 x - 10 y
 % dy/dt = x - y - (sqrt(27)+x)z
 % dz/dt = sqrt(27)*(x+y) - z + xy
 function dy = lorenz(t,x)
    dy = zeros(3,1);
    dy(1) = -10*x(1) + 10 * x(2);
    dy(2) = x(1) - x(2) -(sqrt(27)+x(1))*x(3);
    dy(3) = sqrt(27)*(x(1)+x(2))-x(3)+x(1)*x(2);
 end