I am dealing relatively often with vector functions of a variable vector. Specifically these are (in general nonlinear) mappings from R^6 to R^6, and often I have one or more extra parameter that controls the map.
Often I'd like to get a first-order expansion of these. The general scheme I use to do that is to map mtaylor over the components of the function vector. Like so:
map(mtaylor,F(<x1...x6>),[x1...x6,p1...pn],order);# F(<...>) is <F1(<..>)..F6(<...>)>
This works fine until I hit a situation where the Taylor expansion for one of the parameters p does not exist. Then it bombs. Actually, I often can use assumptions on the parameters to prevent bombing; however, in that case I often get the unchanged function F returned. This is Bad as the whole thing happens in a loop with many concatenated functions F, and now expression swell makes Maple lock up.
In at least a significant subset of these functions, the series command works where taylor fails. "mseries", however, does not exist in Maple. So my question is whether someone has written an "mseries" command for Maple. I know how to do it in principle (but caution is needed to make it work as intended), however, I wonder whether either someone has done it before and might share his/her code, or whether there is a reason why "mseries" actually cannot work.