Question

请用matlab解答

Answer 1

该微分方程可以dsolve函数，解出其解析解。

clc

syms x(t)

Dx=diff(x);

x=dsolve(Dx==0.02*(1-x/600)*x,x(1800)==10)

t=1800:2000;

x=eval(x);

plot(t,x,'r-')

xlabel('t'),ylabel('x(t)')

title('微分方程的解析解') 

运行结果

x =600/(exp(log(59) - t/50 + 36) + 1)