Added code to test how discretising works with a simple PDE