Gave current working version of dirirchlet boundary script a more descriptive name and started working on the neumann boundary conditions using a copy of the previous script.