PDE interpolation approximation
next commit, adding faster approximation and linearity