How to re-index summation bounds in mupad

I am attempting to find the series solution to and ODE with mupad and I am trying to make mupad reindex the bounds on the solution when trying to substitute n=n-1 or n=n-2 several steps in and change the starting bound by + 1 in the derivative of the assumed solution

y = x^m*sum(a(n)*x^n, n=0..infinity) , for example when m = 0 the a(0) term should disappear in y' thus changing the bounds from n=0..infinity to n=1..infinity.

I can change the bounds using subs a couple of times but is there a way for mupad to change the bounds when taking the derivative of y and re-index the bounds and exponents when necessary

If anybody knows how to accomplish this your efforts would be appreciated, thank you.

Answer 1

[Walter Roberson]

The MuPAD operation subsop can be used to alter the bounds. I don't think there is any direct MATLAB interface to this, however.