integrate[exp(-t(t a^2-u b^2+t b^2)/(a^2 b^2)), t]