Maxima没有找到差异解决方案

Question

在我的wxmaxima 11.08.0（Ubuntu的12.04，千里马版本：5.24.0）
的练习我跟着一个example from P.Lutus和他的第二个例子不适合我。

eq: y(t) = -r*c*'diff(y(t),t)+m*sin(%omega*t); 
sol:desolve(eq, y(t)); 
Is %omega zero or nonzero? nonzero 

然后Maxima没有反应，直到我重新启动它。
是否有最大的变化，我需要激活或先定义才能得到结果？

预期输出应为：

有我的问题的第二部分的情况下，我定义用手等式：

sol: y(t) = (m * sin(%omega*t))   /(%omega^2*c^2*r^2 + 1) - 
      (%omega*c*m*r*cos(%omega*t)) /(%omega^2*c^2*r^2 + 1) + 
      (%omega*c*m*r*%e^-((1*t)/(c*r)))/(%omega^2*c^2*r^2 + 1); 

为一个连续过程的初始条件：

init_val:-(c*m*r*(%e^-(t/r*c))*%omega)/(c^2*r^2*%omega^2+1); 
atvalue(y(t),t=0, init_val); 
try2 : desolve(sol,y(t)); 
"Is "%omega" zero or nonzero?" nonzero; 

这里最后一个词还在那里。这些问题是基于使用三角函数的吗？

最好的问候， 马库斯

我通过PPA更新，以wxMaxima 13.04.0 &千里马5.29.1。现在解散了，但最后一个词似乎很复杂。

使用negativ last term和desolve命令执行init_val仍然会在方程中留下％e ^（..）* ...。

Answer 1

在Maxima邮件列表中，您可能会对此问题感兴趣。参见：http://maxima.sourceforge.net/maximalist.html

对于Lutus实施例2的第一个版本，我得到：

y(t) = m*sin(%omega*t)/(%omega^2*c^2*r^2+1) 
-%omega*c*m*r*cos(%omega*t)/(%omega^2*c^2*r^2+1) 
+(y(0)*%omega^2*c^3*r^3+%omega*c^2*m*r^2+y(0)*c*r)*%e^-(t/(c*r)) 
    /(c*r*(%omega^2*c^2*r^2+1))$ 

这是一样的预期结果，如果y（0）= 0。然而，我没有看到在那里假设。

atvalue(y(t),t=0,init_val)后，我得到了相同的结果Lutus，即：

y(t) = m*sin(%omega*t)/(%omega^2*c^2*r^2+1) 
-%omega*c*m*r*cos(%omega*t)/(%omega^2*c^2*r^2+1)$ 

我与千里马5.31.1，与CLISP建工作，在Linux上。