Implementing while loop in Heun's Method

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

回答(1 个)

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.

类别

帮助中心File Exchange 中查找有关 Programming 的更多信息

产品

版本

R2019a

标签

Community Treasure Hunt

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

Start Hunting!

Translated by