Esta es una respuesta incompleta, sólo he cubierto el caso de que un óvalo es "grande" en el sentido de que contiene tanto de sus focos en razonable detalles.
Para simplificar el análisis, vamos a empezar con $(x_1,y_1) = (-1,0), (x_2, y_2) = (1,0)$ y definir
$r, s, u, v$ tal que
$$ \begin{cases}r^2 = (x+1)^2 + y^2\\s^2 = (x-1)^2 + y^2\end{casos}
\quad\text{ y }\quad\begin{cases}u = \frac12 (r-s)\\v = \frac12(r+s)\end{casos}
$$
Vamos a ver en la zona de la mitad de un óvalo en la mitad superior del plano dado por
$$\mathscr{O} = \big\{\; (x,y) \in \mathbb{R}^2 : y \ge 0,\, r + a s \le k \;\big\}$$
WOLOG, podemos suponer $a \ge 1$. Geométricamente, hay 3 posibilidades:
- $k > 2 a$ - El óvalo es "grande" en el sentido que encierra, tanto de sus focos.
- $k = 2 a$ - La elipse se convierte en una forma de lágrima con una cúspide en el foco $(-1,0)$.
- $2a > k > 2$ - El óvalo que encierra uno y sólo uno de los focus $(1,0)$.
En el siguiente análisis, nos concentraremos en el $1^{st}$ caso de que $k > 2 a$.
Aviso de $u^2 + v^2 = \frac12 (r^2 + s^2) = x^2 + 1 + y^2$. En la mitad superior del plano, tenemos
$$
\begin{cases}
x = \frac14 (r^2 - s^2) = uv\\
\\
y = \sqrt{u^2 + v^2 - u^2 v^2 -1 } = \sqrt{(1-u^2)(v^2-1)}
\end{casos}$$
La zona se convierte en elemento
$$dx \wedge dy = \frac{dx \wedge d y^2}{2y} = \frac{duv \wedge d(u^2 + v^2)}{2y}
= \frac{(u dv + v du)\wedge(u du + v dv)}{\sqrt{(1-u^2)(v^2-1)}}\\
= \frac {v^2 - u^2) du \wedge dv}{\sqrt{(1-u^2)(v^2-1)}}
= \frac{( (2v^2-1)+(1-2u^2)) du\wedge dv}{2\sqrt{(1-u^2)(v^2-1)}}
$$
Aviso
$$d u\sqrt{1-u^2} = \frac{1-2u^2}{\sqrt{1-u^2}}\quad\text{ y }\quad
d v\sqrt{v^2-1} = \frac{2v^2-1}{\sqrt{v^2-1}}
$$
Tenemos
$$dx \wedge dy = \frac12 d \left[\frac{u(1-u^2)dv - v(v^2-1)du}{\sqrt{(1-u^2)(v^2-1)}}\right]$$
Al $k > 2a$, la mitad oval $\mathscr{O}$ asignado a un cuadrilátero $\mathscr{Q}$ $uv$- plano.
$$\mathscr{Q} = \big\{\; (u,v) \in \mathbb{R}^2 : -1 \le u \le 1; 1 \le v \le \alpha + \beta u \;\}$$
donde $ \displaystyle \alpha = \frac{k}{a+1},\;\beta = \frac{a-1}{a+1}$.
Esto nos da
$$2 \text{Área}(\mathscr{S}) = 2\int_{\mathscr{O}} dx \wedge dy
= \int_{\parcial\mathscr{Q}} \frac{u(1-u^2)dv - v(v^2-1)du}{\sqrt{(1-u^2)(v^2-1)}}
$$
Es fácil de comprobar en $\partial\mathscr{Q}$, el horizontal y los dos bordes verticales contribuyen en nada.
Como resultado, el área de la plena oval puede escribirse como
$$2 \text{Área}(\mathscr{S}) = \int_{-1}^1 \frac{(\alpha+\beta u)((\alpha+\beta u)^2-1) - \beta u(1-u^2)}{\sqrt{(1-u^2)((\alpha+\beta u)^2-1)}} du\etiqueta{*1}
$$
Abusando de la notación, utilizaremos $v$ como una abreviación de $\alpha + \beta u$ a partir de ahora.
Considere la ecuación cuadrática $t = (\alpha+\beta t)(\alpha t + \beta)$. Al $k > 2a$, se puede comprobar por encima de la ecuación tiene dos raíces reales $\eta$$\eta^{-1}$$t$. Vamos a dejar que $\eta$ ser la única raíz con $|\eta| < 1$.
Definir $\lambda, \mu$ e introducir una nueva variable $z$ tal que
$$\begin{cases}
\lambda &= \alpha + \beta \eta\\
\mu &= \alpha \eta + \beta
\end{casos}
\quad\text{ y }\quad
u = \frac{z + \eta}{1 + \eta z}\;\;\ffi\;\;
z = \frac{u \eta}{1 - \eta u}
$$
Uno puede mostrar que $|\mu| < 1$ y hemos
$$
du = \frac{(1-\eta^2)}{(1+\eta z)^2}dz,\quad
\begin{cases}
u & = \frac{z+\eta}{1+\eta z }\\
1 + u & = \frac{(1+\eta)(1+z)}{1+\eta z}\\
1 - u & = \frac{(1-\eta)(1-z)}{1+\eta z}
\end{casos}
\quad\text{ y }\quad
\begin{cases}
v &= \frac{\lambda + \mu z}{1 + \eta z}\\
v + 1 &= \frac{(\lambda+1)(1 + \mu z)}{1+\eta z}\\
v - 1 &= \frac{(\lambda-1)(1 - \mu z)}{1+\eta z}
\end{casos}$$
Esto implica
$$\begin{align}
v(v^2 - 1) & = \frac{\lambda^2 - 1}{(1+\eta z)^3}(\lambda + \mu z)(1 - \mu^2z^2)\\
u(1-u^2) &= \frac{1 - \eta^2}{(1+\eta z)^3}( z + \eta )(1 - z^2)\\
\frac{du}{\sqrt{(1-u^2)(v^2-1)}} &= \sqrt{\frac{1-\eta^2}{\lambda^2-1}}\frac{dz}{\sqrt{(1-z^2)(1-\mu^2 z^2)}}
\end{align}$$
Aviso de $\beta( 1 - \eta^2 ) = \mu - \lambda \eta = \mu ( 1 - \lambda^2)$, podemos reescribir $(*1)$
$$2 \text{Área}(\mathscr{S}) =
\sqrt{(1-\eta^2)(\lambda^2-1)} \int_{-1}^1 \frac{P(z) dz}{
(1+\eta z)^3 \sqrt{(1-z^2)(1-\mu^2 z^2)}}\etiqueta{*2}$$
donde $P(z) = (\lambda + \mu z)(1 -\mu^2 z^2) + \mu (z + \eta)(1-z^2)$.
Por la fuerza bruta, uno puede mostrar que
$$\frac{P(z)}{(1+\eta z)^3} = A(\eta,\mu) + \left[ B(\eta,\mu) + C(\eta,\mu)\frac{\partial}{\partial \eta} + D(\eta,\mu)\frac{\partial^2}{\partial \eta^2}\right]\frac{1}{1+\eta z}$$
donde
$$\begin{cases}
A(\eta,\mu) &= -\frac{\mu^3+\mu}{\eta^3},\\
B(\eta,\mu) &= \frac{\mu^4 + (\eta^4+1)\mu^2+\eta^4}{\eta^3\mu},\\
C(\eta,\mu) &= -\frac{{\mu}^{4}+\left( -2\,{\eta}^{4}+2\,{\eta}^{2}+1\right) \,{\mu}^{2}-2\,{\eta}^{4}}{{\eta}^{2}\,\mu}\\
D(\eta,\mu) &= \frac{{\mu}^{4}+\left( {\eta}^{4}-4\,{\eta}^{2}+1\right) \,{\mu}^{2}+{\eta}^{4}}{2\,\eta\,\mu}
\end{casos}$$
Sustituir esto en $(*2)$, en términos de los siguientes completa de las integrales elípticas de $1^{st}, 2^{nd}$ $3^{rd}$ tipo:
$$\begin{align}
K(\mu) &= \int_0^1 \frac{dz}{\sqrt{(1-z^2)(1-\mu^2 z^2)}}\\
E(\mu) &= \int_0^1 \sqrt{\frac{1-\mu^2 z^2}{1-z^2}} dz\\
\Pi(\eta^2,\mu) &= \int_0^1 \frac{dz}{(1 - \eta^2 z^2)\sqrt{(1-z^2)(1-\mu^2 z^2)}}
\end{align}$$
Tenemos
$$2\text{Área}(\mathscr{S}) = 2\sqrt{(1-\eta^2)(\lambda^2-1)}
\left(
Un K(\mu) + \left[ B + C \frac{\partial}{\partial \eta} + D + \frac{\partial^2}{\parcial \eta^2}\right]\Pi(\eta^2,\mu)
\right)$$
Por favor, tenga en cuenta que las derivadas parciales de $\Pi(\eta^2,\mu)$ puede ser expresada en términos de la completa de las integrales elípticas a sí mismo.
$$\frac{\partial\Pi(n,\mu)}{\partial n} = \frac{1}{2n(\mu^2 - n)(n-1)}\left(n E(\mu) + (\mu^2 - n)K(\mu) + (n^2 - \mu^2)\Pi(n,\mu)\right)$$
Por lo que en principio, podemos deshacernos de las derivadas parciales de arriba y expresar el área de la elipse únicamente en términos de completar las integrales elípticas directamente.
Nos deja cambiar a la $3^{rd}$ caso de que $2 < k < 2a$. El óvalo ahora cruzan el
segmento de línea que une los dos focos y encierra uno de los focus $(1,0)$.
Deje $\displaystyle \sigma = \frac{1-\alpha}{\beta}$.
En el $uv$-plano, la línea de $v = \alpha + \beta u$ cruzan la línea de $v = 1$ a
$( \sigma, 1 )$ antes de la línea de $u = -1$. La mitad oval $\mathscr{O}$ se correlacionan con un triángulo $\mathscr{T}$:
$$\mathscr{T} = \big\{\; (u,v) \in \mathbb{R}^2 : u \le 1; 1 \le v \le \alpha+\beta u\;\big\}$$
Una vez más, podemos expresar el área de medio óvalo como una integral de línea sobre $\partial \mathscr{T}$ y los bordes horizontales y verticales no contribuye en nada. Como resultado,
el área de la plena oval está dada por
$$2\text{Área}(\mathscr{S}) = \int_\sigma^1
\frac{(\alpha+\beta u)((\alpha + \beta u)^2 - 1) - \beta u (1-u^2)
}{\sqrt{(1-u^2)((\alpha + \beta u)^2-1)}}\etiqueta{*1'}
$$
Desde el límite inferior $\sigma$ ahora depende de $\alpha, \beta$. Necesitamos una diferente
Mobius transformar el mapa que entre el $u$ y una nueva variable $\tilde{z}$ tal que
$$\begin{array}{rcr}
u& &\tilde{z}\\
\frac{-1-\alpha}{\beta} & \longleftrightarrow & -\frac{1}{\tilde{\mu}}\\
-1 & \longleftrightarrow & \frac{1}{\tilde{\mu}}\\
\sigma = \frac{1-\alpha}{\beta} & \longleftrightarrow & -1\\
1 & \longleftrightarrow & 1
\end{array}$$
para llevar el integrando de nuevo a la forma canónica. Si uno realizar esta sustitución,
uno va a obtener algo muy similar a $(*2)$ por encima. La derivación va a ser muy complicado y voy a parar aquí. Sin embargo, no es claro si uno hacer el trabajo sucio a lo largo de esta dirección, en última instancia, vamos a ser capaces de expresar el área de oval de nuevo en condiciones de completar las integrales elípticas.
Actualización
De vuelta al problema original. Observe que la distancia entre el$(2,3)$$(3,1)$$\sqrt{5}$.
Si tenemos la escala de la elipse por un factor de $\frac{2}{\sqrt{5}}$, se reducirá para el caso 3.
hemos cubierto antes con los siguientes parámetros:
$$( k, a ) = (\frac{8}{\sqrt{5}}, 2 )\quad\implica\quad(\alpha,\beta\sigma) =
( \frac{8}{3\sqrt{5}}, \frac13, 3 -\frac{8}{\sqrt{5}})$$
El área de la original oval está dado por la integral en $(*1')$ escalado por un factor de
$\left(\frac{\sqrt{5}}{2}\right)^2$:
$$\text{Area} = \frac{5}{4} \int_{3 -\frac{8}{\sqrt{5}}}^1 \frac{\left( {\left( \frac{u}{3}+\frac{8}{3\,\sqrt{5}}\right) }^{2}-1\right) \,\left( \frac{u}{3}+\frac{8}{3\,\sqrt{5}}\right) -\frac{u\,\left( 1-{u}^{2}\right) }{3}}{\sqrt{\left( {\left( \frac{u}{3}+\frac{8}{3\,\sqrt{5}}\right) }^{2}-1\right) \,\left( 1-{u}^{2}\right) }} du$$
Utilizando esta expresión para WA nos da $$3.2040508691540993294644483966405124836692674716427904\cdots$$