How to solve 2nd order DE with Explicit Euler method?
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Onur Metin Mertaslan
le 25 Mar 2022
Commenté : Onur Metin Mertaslan
le 25 Mar 2022
Hello, I have a second order DE equation with 2 initial conditions. I have to solve it with Eurler's method but I don't kow how to do that in second order DE.
Is there anyone who can share a script or at least show me a way to do second order DE in euler method?
thanks a lot!
0 commentaires
Réponse acceptée
Torsten
le 25 Mar 2022
Modifié(e) : Torsten
le 25 Mar 2022
g = 9.81;
L = 1.0;
T = 1.0;
dt = 0.01;
y_0 = pi/2;
v_0 = 0;
f = @(t,y)[y(2),-g/L*sin(y(1))];
t = (0:dt:T).';
nt = numel(t);
node = 2;
y0 = [y_0 v_0];
y = zeros(nt,node)
y(1,:) = y0;
for it = 1:nt-1
y(it+1,:) = y(it,:) + dt*f(t(it),y(it,:));
end
y_linear = v_0/sqrt(g/L)*sin(sqrt(g/L)*t) + y_0*cos(sqrt(g/L)*t);
plot(t,[y(:,1),y_linear])
But try to start earlier with your assignment next time. Or is it a challenge for you ?
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Mathematics and Optimization 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!