derive(derive(1/(sqr(pi)*a^(3/2))*exp(-r/a),r),r)