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