Calculating growth using ODE 45 from t=0 to t=100

4 次查看（过去 30 天）

Faris Amzar Mohamad Azrai 2020-6-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/541640-calculating-growth-using-ode-45-from-t-0-to-t-100

评论： James Tursa 2020-6-3

Hello, I need help with my code using ODE 45. Can someone help me check my code?

t = 0:0.01:100;
r = [0.05 0.1 0.5 1 10]; 
k = 15; 
B0 = 1; 
B = @(x,y) 30-x.^2-5*(cos((x.*y)/5)).^2;
colourmap = [228 26 28; 55 126 184; 77 175 74; ...
    152 78 163; 255 127 0]/255;
z = zeros(length(t),length(r));
T = zeros(length(t),length(r));
% 
for i = 1:length(r)
    dBdt = @(B) r.*B*(1-B./k);
    [z,T] = ode45(dBdt,tspan, B0);
    
    plotting
    figure(6) %figure 6
    loglog(t,T,'-','color',colourmap(i,:))
    hold on
    grid on
end