Add PolynomialSpace for expansion in general OPs

[dlfivefifty]

Essentially L is the conversion operator between new polynomial space and Chebyshev U. So we could plug in to the ApproxFun ODE solving. Will have to think of a good use