Question: Differentiating output of dsolve

Good Day!

I numerically solved a system of ODEs with dsolve for the functions p(x) and q(x) over the domain x=[-1,+1].  Boundary conditions are specified at x=-1 and x=+1.  The output is a listprocedure.  

I want to determine the n-th derivatives of p(x) and q(x) at x=0.  Alternatively, I want the power series representation of p(x) and q(x) about x=0.

I would appreciate any advice on how to use Maple 10 to accomplish this.

Thanks.

Glenn

Please Wait...