support for higher order boundary conditions

This should allow support for higher order boundary conditions by using the values of the higher order derivatives of the Chebyshev polynomials at +-1 and scaled by the domain.
Convenience routines to Evaluation: ldiffbc, rdiffbc, and diffbcs.