(b^2/a^2+a^2/b^2-2)*((a+b)/(b-a)+(b-a)/(a+b))*([1/(a^2)+1/(b^2)]/[1/(b^2)-1/(a^2)]-[1/(b^2)-1/(a^2)]/[1/(a^2)+1/(b^2)])