arcsin((r-13)/r)+arcsin((r-16)/r)=pi/2-arcsin((r-20)/r)