function angle = phi_i(r, focus)
    angle = atan(r / focus);
end