Riga 1: Riga 1:
%Equazione differenziale Ordinaria dell'attrattore di lorenz
+
  %Equazione differenziale Ordinaria dell'attrattore di lorenz
%-----------------------------------------------------------
+
  %-----------------------------------------------------------
% dx/dt = -10 x - 10 y
+
  % dx/dt = -10 x - 10 y
% 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);
    dy(1) = -10*x(1) + 10 * x(2);
+
    dy(1) = -10*x(1) + 10 * x(2);
    dy(2) = x(1) - x(2) -(sqrt(27)+x(1))*x(3);
+
    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);
+
    dy(3) = sqrt(27)*(x(1)+x(2))-x(3)+x(1)*x(2);
end
+
  end

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