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(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 | + | 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
%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