nijso beishuizen <nijso at hotmail.com> writes:
> The two methods come up with different constants of integration,
> without expand it is y^2/2 and with expand it is 0.
> I use integrate twice to compute the integration factor for
> an ode somewhere and the first time the integration is over x
> and the second time over y. If I don't use expand then this goes
> wrong.
> so in the end I do this:
> return(ratsimp(integrate(expand(mu),y) - integrate(expand(mu*phi),x)=%c))
I'm not sure exactly what you're doing, but it seems that if you care
about the constant of integration, you probably shouldn't be using an
indefinite integral. Why not do something like
integrate( subst( u, x, expression ), u, 0, x );
This forces the integration to have what I suspect you're calling a zero
constant of integration.
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/20091112/cc64d471/attachment.pgp