Implementing while loop in Heun's Method

2 visualizzazioni (ultimi 30 giorni)
Nisrina Younes
Nisrina Younes il 12 Nov 2022
Modificato: Alan Stevens il 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

Risposte (1)

Alan Stevens
Alan Stevens il 12 Nov 2022
Modificato: Alan Stevens il 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.

Categorie

Scopri di più su Programming in Help Center e File Exchange

Prodotti


Release

R2019a

Community Treasure Hunt

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

Start Hunting!

Translated by