Let us put the equation of the ellipse under its canonical form :
$$\frac{x^2}{(\sqrt{13})^2}+\frac{y^2}{2^2}=1\tag{1}$$
Therefore, the abscissas of its foci are $f=\pm \sqrt{a^2-b^2}=\pm 3$.
A recall : the locus of points "seing" a same segment (here $FF'$) under a same (non-oriented acute) angle $\theta$ is the curve $C_{\theta}$ which is the union of 2 circular arcs as represented here connected in points $F,F'$. These curves are "embedded" like russian dolls. The smaller the angle, the larger the curve. We are looking here for the larger value of $\theta$ for which curve $C_{\theta}$ is tangent to line $L$.
Have a look at the figure below featuring the two points on the line which are the two extremal solutions :

Fig.1 : The two points (green dots) which are extremal solutions. Strictly speaking, it is the upper one which exhibits the larger angular value.
Here is a description of the way these points have been obtained :
- Step 1: It is not difficult to see that the generic equation of the circles passing through both foci $F(-3,0)$ and $F'(3,0)$ is
$$x^2+y^2-2dx-9=0\tag{2}$$
where the center of the circle is in $(0,d)$
(such a set of circles is called an "elliptical pencil of circles")
- Step 2 : You probably know the classical rule giving the equation of the tangent of a circle with an equation like (2) in a point $(x_0,y_0)$ is by polarizing it (I will look for a hopefuly more convenient term in English) in the following way :
$$xx_0+yy_0-d(y+y_0)-9=0$$
(sanity check : suppressing indices $_0$, one finds back equation (2).)
Otherwise said :
$$xx_0+y(y_0-d) -(dy_0+9)=0$$
- Step 3 : the former equation is eqivalent to the given equation $2x+y-12=0$ iff their coefficients are proportional:
$$\frac{x_0}{2}=\frac{y_0-d}{1}=\frac{dy_0+9}{12}\tag{3}$$
As $(x_0,y_0)$ belongs to line $L$, i.e., we have as well :
$$2x_0+y_0-12=0\tag{4}$$
- Step 4: Equations (3)+(4) constitute a system of 3 equations with 3 unknowns ; it is not difficult to solve it ; here are the two solutions (with the coordinates $(x_0,y_0)$ of the points and the ordinates $d$ of the centers of the corresponding circles (where we have set $k=\sqrt{15}$).
$$\begin{cases}x_0& =& -(3/5)k + 6, \ \ \ \ &y_0 &=& +(6/5) k, \ \ \ \ &d &=& +(3/2) k - 3 \\ x_0 &=& +(3/5)k + 6, \ \ \ \ &y_0 &=& -(6/5) k, \ \ \ \ &d &=& -(3/2) k - 3
\end{cases}$$
As said in the legend of Fig. 1, only the first point should be finally considered.
- Step 5: Once the coordinates of the optimal point $P_0$ have been obtained, it is easy to compute the ratio $P_0F/P_0F'$ ; one finds $\sqrt{3}$.
For those who are interested in the way the figure has been drawn, here is its "Sage" program :
var('x y d')
s=solve([2*x+y==12,x/2==y-d,y-d==(d*y+9)/12],x,y,d)
x0=s[0][0].rhs();y0=s[0][1].rhs();
x1=s[1][0].rhs();y1=s[1][1].rhs();
print(s)
L=13
g=line([(-L,0),(L,0)],color='black')
g+=line([(0,-L),(0,L)],color='black')
d=s[0][2].rhs();g+=implicit_plot(x^2+y^2-2*d*y-9==0,(x,-L,L),(y,-L,L))
d=s[1][2].rhs();g+=implicit_plot(x^2+y^2-2*d*y-9==0,(x,-L,L),(y,-L,L))
g+=implicit_plot(2*x+y==12,(x,-L,L),(y,-L,L),color='green')
g+=implicit_plot(x^2/13+y^2/4==1,(x,-L,L),(y,-L,L),color='red')
g+=point(((-3,0),(3,0)),color='red',size=50)
g+=point(((x0,y0),(x1,y1)),color='green',size=50)
show(g)