Added PDE model

IC_sane gives initial conditions, PDE_implementation is a finite
difference implementation of the system, main runs the simulation from
PDE_implementation and plots results