Weird mathematica problem

What don't you like about this representation? If you want to have the terms ordered by the total degree of the monomials, there seems not to be a built-in way to do it. (But it shouldn't be too hard to come up with a work-around oneself ...)

 

Try using the function Collect:

Collect[Series[Exp[x y], {x, 0, 8}, {y, 0, 8}] , {x, y}]

If that's not what you want, then I fail to understand. I'm sure there must be a way to do what you are asking, however.

 

Now I understand the problem. That's strange that Mathematica does not have much material about multivariate series, so we can add it ourselves.

Code (Text):

multiVarSeries[f_, x_List, a_List, k_Integer] :=  Block[{n, F},
  Evaluate[Fold[Sum,
     Product[(1/n[i]!) (x[[i]] - a[[i]])^(n[i]), {i, 1, Length[x]}]
 ((Fold[D, F@@x, Table[{x[[i]], n[i]}, {i, 1, Length[x]}]]) /.
        Table[x[[i]] -> a[[i]], {i, 1, Length[x]}]),
     Table[{n[i], 1, k}, {i, 1, Length[x]}]]] /. F -> Function[x, f]]

Copy and paste that function into a new cell and then execute it with shift + enter. After that you can invoke the function. Here is a simple example:

multiVarSeries[Exp[x y], {x, y}, {0,0}, 2]

This says to expand the function Exp[x y] with respect to the variables x and y around the point {0,0} up to order 2 (in both variables, I don't let you specify the order separately for each individual variable). The output is of course:

[tex]\frac{x^3 y^3}{6}+\frac{x^2 y^2}{2}+x y[/tex]

 

Yes, I understand your complaint with the function. One thing you can do is make n larger than you need and then use a filter to get only the terms you want. It is more work than it is worth for me to change the function to match that behavior.

There is no good way. If you are really wealthy, you can do workshops online with Wolfram Inc that will teach you how to program Mathematica. Otherwise you have to do what I did, which is to read the built-in help and practice for months.

 

What doesn't work and what function are you trying to do this with?