调用maple
>> maple('dsolve','{diff(N(t),t)=n(t)*r*(1-N(t)/M),diff(n(t),t)=n(t)*s*(1-n(t)/K)}')
ans =
[{n(t) = K/(1+exp(-s*t)*_C2*K)},{N(t) = (exp(r/M*Int(n(t),t))*M+_C1)*exp(Int(-n(t)*r/M,t))}]
可以验证n(t)
>> n=dsolve('Dn=n*s*(1-n/K)')
n =
K/(1+exp(-s*t)*C1*K)
只要给出参数的具体数值和初始条件,有解.