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); |
%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