と思ってまた問題を見ると、求めるべき角度が違った。
θ=acos((b^2+r^2-a^2)/(4*b*r))