If $x^3+y^3-3\alpha xy=1$, is there an expression for the integral $$\int_0^z \frac{\mathrm dx}{y^2-\alpha x}$$ in terms of more familiar functions?

A.C. Dixon introduced the elliptic functions $\operatorname{sm}(u,\alpha)$ and $\operatorname{cm}(u,\alpha)$ now named after him in this article. It can be shown (see e.g. my writeup here) that these functions can be expressed in terms of the more conventional Weierstrass elliptic functions, e.g.

$$\operatorname{sm}(u,\alpha)=-\frac{2\wp\left(u;g_2,g_3\right)+\frac{\alpha^2}{2}}{\wp^\prime\left(u;g_2,g_3\right)+\alpha\wp\left(u;g_2,g_3\right)-\frac{\alpha^3+4}{12}}$$

where

$$\begin{align*} g_2&=\frac{\alpha}{12}\left(\alpha^3-8\right)\\ g_3&=\frac{8-20\alpha^3-\alpha^6}{216} \end{align*}$$

The problem stated above is then equivalent to asking for a (hopefully simpler) explicit expression for the inverse Dixon elliptic function, $\operatorname{sm}^{(-1)}(z,\alpha)$.

Lagrangian inversion of the Maclaurin series for $\operatorname{sm}(u,\alpha)$ yields the series

$$z-\frac{\alpha z^2}{2}+\frac{5 \alpha^3+2}{12} z^4-\frac{\alpha\left(7 \alpha ^3+4\right)}{15} z^5+\frac{44 \alpha^6+40 \alpha^3+5}{63} z^7+c_8 z^8+\cdots$$

where the Maclaurin coefficients $c_n$ satisfy the recurrence

$$n(n+1)(n+2)c_n-(n+3)\left(18 \alpha ^3+4 \alpha ^3 n^2+n^2+18 \alpha ^3 n-3 n-24\right)c_{n+3}-(n+6)\left(180 \alpha ^3+\left(4 \alpha ^3+1\right) n^2+3 \left(18 \alpha ^3+7\right) n+84\right)c_{n+6}+(n+7)(n+8)(n+9) c_{n+9}=0$$

and with a little more work, one can derive the differential equation satisfied by $w=\operatorname{sm}^{(-1)}(z,\alpha)$:

$$(z+1)\left(z^2-z+1\right)\left(z^6-2\left(2 \alpha^3+1\right) z^3+1\right)w^{(3)}(z)+6 z^2\left(z^6+\left(1-\alpha^3\right) z^3-3\alpha^3-2\right)w^{\prime\prime}(z)+2z\left(3z^6+\left(\alpha^3+7\right)z^3-5\alpha^3-2\right) w^\prime(z)=0$$

Unfortunately, I have not succeeded in making further headway. I have reason to suspect that a (generalized) hypergeometric function (e.g. Appell's $F_1$) is involved, considering that the special value $\operatorname{sm}^{(-1)}(1,\alpha)$ is expressible in terms of the Gaussian hypergeometric function:

$$\operatorname{sm}^{(-1)}(1,\alpha)=\frac13 B\left(\frac13,\frac13\right) {}_2F_1\left({{\frac13,\frac13}\atop{\frac23}}\middle|-\alpha^3\right)+\frac{\alpha}{3} B\left(\frac23,\frac23\right) {}_2F_1\left({{\frac23,\frac23}\atop{\frac43}}\middle|-\alpha^3\right)$$

where $B(x,y)=\dfrac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}$ is the usual beta function.

This has been out for a while, but I have not yet gotten anything that is computationally useful (or even something I can numerically experiment with).

So, to spell things out *very explicitly*:

Given the standard Weierstrass reduction outlined above for $x^3+y^3-3\alpha xy=1$, can someone please write out the explicit representation of the inverse Dixon function $\operatorname{sm}^{(-1)}(z,\alpha)=\int_0^z \frac{\mathrm dx}{y^2-\alpha x}$, perhaps in terms of standard elliptic integrals or multivariate hypergeometric functions?

(I can worry about implementation details later; I just need a starting point.)

2more comments