Fix in pseudospectral.py: address case where evaluation point is one of the interpolation points.
Did it: Added on age constraint in ad hoc box in middle of domain.  Adding more should be easy.  Evaluating performance of the method may be hard ...