improved euler integration using the butterfly effect
1 vue (au cours des 30 derniers jours)
Afficher commentaires plus anciens
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 commentaire
Torsten
le 26 Mar 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.
Réponses (0)
Voir également
Catégories
En savoir plus sur Function Creation dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!