improved euler integration using the butterfly effect
3 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
i have created the following code but it appears as error. can someone please help me? thank you
%euler's Method
% Initial conditions and setup
sigma=10;
beta=8/3;
rho=28;
t1=100;
t0=0;
x0=0;%initial x value
y0=1;%initial y value
z0=0; %initial z value
Delta=0.0001;
t=t0:Delta:t1;
i=1;
x=1*ones(1,length(t));
while (i<length(t))
x(i+1)=sigma*(y-(1+Delta));
i=i+1;
end
y=1*ones(1,length(t));
while (i<length(t))
y(i+1)=x*(rho-z)-(1+Delta);
i=i+1;
end
z=1*ones(1,length(t));
while(i<length(t))
z(i+1)=x*y-beta*(1+Delta);
i=i+1;
end
plot(x,y,z);
1 comentario
Torsten
el 26 de Mzo. de 2020
In the first while loop, y is undefined.
In the second while loop, z is undefined.
Before entering the second and third while loops, you forgot to reset i to 1.
I further think that you will have to use x(i-1), y(i-1) and z(i-1) instead of x,y and z on the right-hand sides of the expressions within the while loops.
If errors persist, please include the error messages.
Respuestas (0)
Ver también
Categorías
Más información sobre Loops and Conditional Statements en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!