moved integral code to the last iteration of each Newton solve, rather than all iterations, also changed timing code to print # iterations and time