Julien Cubizolles <j.cubizolles at free.fr> writes:
> I just noticed that some old (3-4 years old) maxima code I had written
> doesn't work anymore. It's using the dynamics package
>
> Here it is
>
>
> /* [wxMaxima: input start ] */
> load(dynamics)$load(draw)$
> /* [wxMaxima: input end ] */
>
> /* [wxMaxima: input start ] */
> Epot(x):=.5*x^2$
> F(x):=-diff(Epot(x),x)$
> /* [wxMaxima: input end ] */
>
> /* [wxMaxima: input start ] */
> context :aperiodique$assume(Q>0,Q<1/2)$
> /* [wxMaxima: input end ] */
>
> /* [wxMaxima: input start ] */
> context : pseudoperiodique;assume(Q> 1/2);
> /* [wxMaxima: input end ] */
>
> /* [wxMaxima: input start ] */
> x0 : 1$ v0 : 0$
> /* [wxMaxima: input end ] */
>
> /* [wxMaxima: input start ] */
> kill(Q)$context : pseudoperiodique $
> sol : ic2(ode2('diff(x,t,2)+'diff(x,t)/Q+x=0,x,t),t=0,x=x0,dy=v0);
The problem is happening in ic2 here. Are you sure this ever worked? I
would expect something like 'diff(x,t) = v0 rather than dy=v0.
Rupert
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 315 bytes
Desc: not available
URL: <http://www.math.utexas.edu/pipermail/maxima/attachments/20121217/bea1bc5c/attachment.pgp>;