+ results printing to csv file
+ to do: find boundary conditions that give limited high order partial
derivatives values