我正在学习如何使用wxMaxima来解决微分方程,但我已经遇到了非常简单的ODE问题。例如,如果我想解决逻辑微分方程并使用ode2:
diffeq: 'diff(S,t)=g*S*(1-S/K);
ode2(diffeq,S,t);
千里马回归
(log(S-K)-log(S))/g=t+%c
而且我不知道为什么Maxima不能解决S(t)或者我如何能够以S(t)= xxx的形式获得简单的解。我试图使用logexpand()等但无法摆脱log()。有没有办法获得更直接的解决方案,或者,有人可以推荐另一个免费软件来处理微分方程吗?
谢谢!
我可以说服Maxima生产一种更直接可用的解决方案,只需要一点点刺激。我会让你决定这是否有用。
(%i1) eq: (log(S-K)-log(S))/g=t+%c;
log(S - K) - log(S)
(%o1) ------------------- = t + %c
g
(%i2) logcontract (%);
S - K
log(-----)
S
(%o2) ---------- = t + %c
g
(%i3) map (exp, %);
S - K
log(-----)
S
----------
g t + %c
(%o3) %e = %e
(%i4) lhs(%)^g = rhs(%)^g;
S - K g (t + %c)
(%o4) ----- = %e
S
(%i5) solve (%, S);
K
(%o5) [S = - ----------------]
g t + %c g
%e - 1
至于其他软件,你可以试试Sympy(sympy.org)。我对它如何处理微分方程一无所知。