Add third density derivatives to sat_Z_accuracy.cpp (remarkably, even better than the first derivative)

Evidently not returning a value from a function that should causes some weird behaviors with (ierr != 0) not working properly