Defined a demand function:
The Slutsky Matrix is calculated as Dpx(p,w) + Dwx(p,w) x(p,w)^T with Dpx(p,w) = Jacobian for p1, p2, p3 and Dwx(p,w) is a column vector of the derivate of x1, x2, x3 in respect to w.
Step by step:
Dxp := Jacobian(x(p,w),[p,p,p])
yields the correct Jacobian.
For Dwx I've tried several different methods, i.e.
Dwx := Jacobian(x(p,w),[w]) or diff(x(p,w),w) or D(x)(p,w)
which all appear to yield the correct column vector.
The problem is that Maple tells me 'Error, (in rtable/Product) invalid arguments), when I try the product:
If I instead define Dwx through:
Dwx := Vector(3,(i)->diff(x(p,w)[i],w))
Maple correctly calculates Dwx x(p,w)^+ as a 3x3 matrix.