symmetric skew Macdonald polynomials

edit

Is there a correct way to compute skew Macdonald polynomials in sage? The routine skew_by does not seem to give the right answer (far from it). Probably due to the fact that "zee" is not defined properly somewhere.