D[Q0*exp(-r/h)*(cosh[z/z0])