       Help with project needed

• To: mathgroup at smc.vnet.net
• Subject: [mg96619] Help with project needed
• From: Aaron Fude <aaronfude at gmail.com>
• Date: Tue, 17 Feb 2009 06:27:26 -0500 (EST)

```Hi,

I'm about to attempt a project in which so many Mathematica related
things are unclear to me that I will describe the project, rather than

Suppose that C is a function of alpha, but it also depends on a
parameter epsilon, so I think of it as C[epsilon][alpha] -- but I'm OK
with it if the implementation treats C as a function of two variables.
This C is given. It's complicated, but one could easily obtain a
series for it in epsilon.

Now, I have the following operator that acts on functions of "r", and
"alpha" (d is partial):

L = C''(1-C')/(1+C^2) * d/dr + C'(1-C'') /(1-C')* d/dalpha

For a given function f[r, alpha], I need to take successive orders in
epsilon of the operator and apply it to f.

I realize that for a given function  f, I can form L[f] and the do
Series[L[f], {epsilon, 0, 5}]. However, I'm very interested in seeing
the operator decomposed into a series. I want to be able to say,
here's the epsilon term of the operator, here's the epsilon squared
term, and so forth. Can this be done. Can I decompose L into a series:

L = L0 + L1 + L2 + ....

and then be able to apply Ln to f?