Implementing while loop in Heun's Method

2 vues (au cours des 30 derniers jours)
Nisrina Younes
Nisrina Younes le 12 Nov 2022
Modifié(e) : Alan Stevens le 12 Nov 2022
For dh/dt=-k√H where k is = 0.06
how can I implement while loop (while H>0) in Heun's method to find the time taken for the tank to drain ?. given the initial value of H is 4
Thanks for the help in advance

Réponses (1)

Alan Stevens
Alan Stevens le 12 Nov 2022
Modifié(e) : Alan Stevens le 12 Nov 2022
Here's a rough and ready way (I'll leave you to implement Heun's method):
f = @(H) -0.06*sqrt(H);
H = 4;
t = 0;
dt = 0.1; % choose your own value
while H>0
t = t+dt;
% Use Heun's method to get H at t+dt
end
disp(t) % time to drain
Note that, because of the square root, the last step is likely to result in a complex value for the final value of H.

Catégories

En savoir plus sur Programming dans Help Center et File Exchange

Produits


Version

R2019a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by