“…The algorithm in Section 2 finds the root β λ of the function F λ (β) in (12), which together with (x (0, λ, β λ ), y(0, λ, β λ )) give the initial condition of a level curve of u which solves (6). To find β λ numerically, we use the boundary points which are joined by the level curve of u corresponding to λ, say (x 1 , y 1 ) and (x 2 , y 2 ).…”