We search a given line for a point that yields a summation of angles (for both zero and pole angles) equal to an odd multiple of \$180^\circ\$. Take a look at the following picture
The line in which a point that yields a summation of angles equal to an odd multiple of 180 deg is shown as the Green line in the above picture. The angle of this line is \$ \beta = \cos^{-1} (\zeta) \$ and the radius that satisfies the aforementioned condition is searched by a software that we code. In this example, this is Matlab code that searches for the target point.
r=.75;
while r < .76
zeta=.5;
a=acos(zeta);
x=r*cos(a);
y=r*sin(a);
%poles' angles
p1=0;
Th1=pi-atan( y/x );
p2=-1;
Th2=atan(y/(abs(p2)-x));
p3=-3;
Th3=atan(y/(abs(p3)-x));
angle=-(Th1+Th2+Th3);
fprintf('Th1=%.3f : Th2=%.3f : Th3=%.3f : angle=%.3f : r=%.3f\n', rad2deg(Th1),rad2deg(Th2), rad2deg(Th3),rad2deg(angle),r);
r = r + .00001;
end
If you run the code, you notice the following result
Th1=120.000 : Th2=46.102 : Th3=13.898 : angle=-180.000 : r=0.750
The aforementioned condition is satisfied if \$\theta_1 = 120^\circ, \theta_2=46.102^\circ \$ and \$\theta_3 = 13.898^\circ\$ which yields a radius \$r=0.75\$. The point \$P\$ is then \$ (r\cos(\beta), r\sin(\beta))\$, or \$P=-0.3750+j0.6495\$. Now, to compute the gain K, we substitute P in the following formula:
$$ \begin{align} K &= \frac{1}{|G(s)||H(s)|} \\ &= s(s+1)(s+3) \Big|_{s=P} \\ &= 1.8281. \end{align} $$$$ \begin{align} K &= \frac{1}{|G(s)||H(s)|} \\ &= |s(s+1)(s+3)| \Big|_{s=P} \\ &= 1.8281. \end{align} $$