[Richard Fateman]
>>> Here's a suggestion.
>>> Look at
>>> http://www.cs.berkeley.edu/~fateman/papers/partition.pdf
>>>
>>> Write a program gammasimp(p) which does the following [...]
[Barton Willis]
>>
>> There is a function gatherargs that might be helpful. To use gatherargs,
>> first
>> do load("opsubst"). ?The quotes are important :(
>>
>> --Barton
[David Billinghurst]
> There are simplification functions for Bessel and Exponentional Integral
> functions that use gatherargs and then partition the expressions based on
> arguments in share/contrib/diffequations/contrib_ode.mac.
>
> ? ?David
Thanks everyone. You have been extremely helpful. I am a bit short
of time at the moment; but in a few weeks' time I hope I will be able
to revisit the problem. The idea sketched in Richard Fateman's e-mail
and paper seems far better than my naive pattern matching (I suppose I
might still need to find a solution for infinite loops of
substitutions).
Thanks again,
Pouya