Switched to doing a non-polynomial solution on a non-reference element.  This converges to an incorrect solution.
Added dz() implementations for several Function subclasses.