diff(atan(b*s/(a*sqrt(s^2+a^2+b^2))),s)