While solving this post, I came across the Somos-5 sequence again (since it appeared in another post). Given the Heronian triangle with half-sides $(a,b,c)$, area $A$, and medians $(m_1,m_2,m_3)$,
$$A=\sqrt{(a + b + c) (b + c - a) (a + c - b) (a + b - c)}$$ $$m_1=\sqrt{2 (b^2 + c^2) - a^2}$$ $$m_2=\sqrt{2 (c^2 + a^2) - b^2}$$ $$m_3=\sqrt{2 (a^2 + b^2) - c^2}$$
Let $(a,b,c)$ be integers. By solving an elliptic curve, it is possible to make all $m_k$ as integers. However, if $A$ is to be a non-zero integer as well, that is still an open problem.
If we start with the area $A$, Schubert conjectured in 1905 that solving the first three equations in non-zero integers was impossible, but was proven wrong by Buchholz's 1989 thesis.
I. Three quadratics
We do the cyclic substitution $(a,b,c) = (p+q,\,q-r,\,r-p)$ to simplify the equations,
\begin{align} |16\, p\, q\, r\, (-p - q + r)| &=\, A^2\\[5pt] (p - q)^2 - 4 (p + q) r + 4 r^2 &= m_1^2\\[5pt] (2 p + q)^2 - 2 (2 p - q) r + r^2 &= m_2^2\\[5pt] (p + 2 q)^2 + 2 (p - 2 q) r + r^2 &= m_3^2 \end{align}
So four quadratics in $r$ to be made squares. For this system, three at a time is infinitely solvable.
II. Somos-5 Sequences
Define the two Somos-5 sequences (starting with $n=0$),
$S_n = 1,\, 1,\, 1,\, 2,\, 3,\, 5,\, 11,\, 37,\, 83,\, 274,\, 1217,\; 6161,\; 22833,\dots$ (A006721)
$T_n = 0, 1, -1, 1, 1, -7, 8, -1, -57, 391, -455, -2729, 22352,\dots$ (A360381)
Then in page 6 of Andrew Hone's 2022 paper, we find the formulas (modified),
$$A = |4\alpha\beta|$$
$$\alpha = S_{n}\,S_{n+1}\,S_{n+2}^2\, S_{n+3}\, S_{n+4}$$ $$\beta = T_{n}\,T_{n+1}\,T_{n+2}^2\, T_{n+3}\, T_{n+4}$$
and the $(p,q,r)$ as,
$$p_n = S_{n+1}\,S_{n+2}^3\, S_{n+3}\, T_{n+2}$$
$$q_n = \;S_{n}^2\, S_{n+1}\, T_{n+3}\, T_{n+4}^2$$
$$r_n = T_{n+1}\, T_{n+2}^3 \,T_{n+3}\, S_{n+2}$$
$$s_n = -\,T_{n}^2\, T_{n+1}\, S_{n+3}\, S_{n+4}^2$$
with the formula for $s_n$ added by yours truly. So these $(A,p,q,r,s)$ are highly factorable, being products of 12 and 6 terms, respectively (counting multiplicity). Specifically,
$$p_n = -2,\, 24,\, 270, -28875, 1969880, -46246189,\dots$$
$$q_n\; =\; 1,\; 49,\; -896,\; 96,\; -146205, -2396409675,\dots$$
$$r_n\; = -1,\, -2,\; -21, -13720,\; 39424,\; 16872,\dots\;\quad$$
$$\quad s_n \,=\; 0,\; 75, -605, -15059, 1784251, -2442672736,\dots$$
Example:
$(p_1,q_1,r_1) = (24,\,49,-2)$ solves the first three quadratics above while $s_n \to r_n$ $(p_1,q_1,s_1) = (24,\,49,\,75)\,$ solves the first two and the fourth. And so on for all $n$.
III. Question
Q: What are the recurrence relations of the $(p,q,r,s)$? Are they also Somos-5 sequences but with different initial values?
P.S. Part 2 in in this post.