Code is running and seems stable. Need a relatively large number of cells per axis dimension for an accurate result.