Partial differential equation solvers written in Haskell