Info

此问题已关闭。 请重新打开它进行编辑或回答。

What did I do wrong?

1 次查看（过去 30 天）
Zifeng Qiu
Zifeng Qiu 2020-7-6
关闭： Stephen23 2020-7-6
I am trying to solve the second question which is k in the given model, what did I do wrong?
function [tv, uv] = RK4(f, u0, T, n)
dt = T/n; % Differential of time.
tv = transpose(0:dt:T); % Evaluation times.
uv = zeros(n+1, 1); % Solution.
uv(1) = u0; % Initial value
for i=1:length(tv)-1
k1 = f(tv(i), uv(i))
k2 = f(tv(i) + (dt/2), uv(i) + dt*k1/2)
k3 = f(tv(i) + (dt/2), uv(i) + dt*k2/2)
k4 = f(tv(i) + dt, uv(i) + dt*k3)
uv(i+1) = uv(i) + dt/6*[k1+2*k2+2*k3+k4]
end
end
u0 = 5000;
lambda = 0.03;
pm = 9000;
f = @(t,k) lambda*500*(1-500/pm)-k;
[t,k] = RK4(f,u0,100,10)
Kf = k(end)

请先登录，再回答此问题。

回答（0 个）

标签

Community Treasure Hunt

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

Start Hunting!

Translated by