\section{Introduction}\label{sec:intro}
\begin{figure*}
\centering
\includegraphics[scale=1]{./figures/Fig1.eps}
\caption[]{\footnotesize \textbf{Bending rules of swimming fish.} (A) Schematic of flexion parameters in fish with flexion distance $l$, ratio $l/L$ and maximum angle defined as in \cite{Lucas2014}. (B) Flexion ratio $l/L$ and maximum flexion angle $A$  observed in different fish species; data taken from \cite{Lucas2014}. The observed flexion ratio and angles are fairly consistent among different fish species, despite large variations in length scales.
}
\label{fig:biodata}
\end{figure*}
Millions of years of natural selection endowed fish with remarkable abilities to swim efficiently compared to underwater man-made propulsors~\citep{Sfakiotakis1999,Lauder2005}.
Several fish use body and caudal fin (BCF) deformations for propulsion. Details of BCF deformations have been used to classify fish swimming modes \citep{Sfakiotakis1999,Shadwick2005,Low2010,Smits2019}.
Invariably, in most BCF swimming modes, anterior to posterior bending of the fish body seems to play an important role in swimming efficiency, and it is often linked to increased flexibility towards the fish tail or caudal peduncle \citep{Combes2003,Tytell2004,Tytell2010,Tytell2016,Gemmell2016}.

In an effort to document the bending rules of fluid-based propulsors, both aerial and aquatic, \cite{Lucas2014} collected morphometric data of the flexion parameters across length scales and animal taxa. They identified two parameters to characterize the bending behavior: flexion ratio $l/L$ between the flexion distance $l$ and total length $L$ of the propulsor and maximum flexion angle $A$. They found that flexion ratio and maximum flexion angle of all surveyed animals, including fish, clustered in a limited design space: bending occurs at about $70 \%$ of body length at maximum flexion angle of about $30 ^\circ$, as shown in Figure~\ref{fig:biodata} for swimming fish based on the data collected in \cite{Lucas2014}.
It is unclear the extent to which this anterior-to-posterior bending is active or whether it follows passively due to the interaction of a flexible posterior with the fluid motion.
Either way, these findings raise the question of whether hydrodynamics could have provided a selective force for driving this convergent bending design.

To address these questions, we analyze the influence of bending on the swimming speed and efficiency of a simplified fish model that consists of anterior and posterior sections connected via a rotational joint at the flexion point (see Figure~\ref{fig:biodata}A). The fish anterior undergoes  periodic planar pitching while the posterior either (i) moves in synchrony with the anterior as if the two parts were a single rigid body, (ii) bends actively at distinct amplitude and phase relative to the fish anterior, or (iii) bends passively due to interactions with the flow generated by the fish anterior. We find that swimming with passive bending could be more efficient than rigid flapping but at the cost of diminished swimming speed. Active bending provides more possibilities to alter the swimming performance through not only the flexion ratio and maximum flexion angle as reported in~\cite{Lucas2014} but also the phase difference $\phi$ between the flapping motions of the anterior and posterior parts. Importantly, we find that antiphase anterior-to-posterior flexion can simultaneously enhance swimming speed and efficiency in a region of the design parameter space $(l/L, A)$. Flexion ratios and angles that lead to significant improvements in speed and efficiency overlap with the observations of real fish reported in \cite{Lucas2014}. We analyze in depth the hydrodynamic mechanisms underlying these improvements in swimming performance.

Details of the flow field around swimming fish have received a great deal of attention.
Several studies used particle image velocimetry to measure the flow field around live fish and analyze the interplay between body deformations and thrust production; see, e.g.,~\cite{Muller1997, Muller2002, Liao2003, Tytell2004, Gemmell2016}.
In-silico models of various degrees of fidelity to fish morphology and kinematics have also been used to examine the offsets of body deformations on swimming speed and efficiency~\citep{Eldredge2006, Kern2006, Tytell2010, Eloy2013}.
Importantly, several experimental and numerical studies have shown that plates and foils undergoing pitching or heaving motions provide good approximations of the fluid-structure interactions in swimming fish \citep{Blondeaux2005,Wen2013,Lauder2011,Menon2019}; including the reverse von K{\'a}rm{\'a}n vortex wake left behind swimming fish and flapping foils \citep{Taneda1965,Triantafyllou1993}.

A variety of fluid-structure interaction models have been proposed to analyze the effect of body flexibility on bending in flows; see, e.g.,~\cite{Shelley2011} and references therein. Here, we present a brief  literature review  focused on this topic.
\cite{Heathcote2007} conducted experiments on a flapper with a rigid leading edge and flexible tail fixed in a water channel and found that flexibility can enhance both efficiency and thrust production. \cite{Eldredge2008} and \cite{Eldredge2010} numerically simulated the flapping motion of articulated rigid links and found that joint flexibility can reduce the power required for flapping.
\cite{Alben2008} used a filament of uniform flexibility to model the tail of swimming fish in the context of the vortex sheet method and predicted enhancement in efficiency rather than thrust when choosing parameters (dimensionless rigidity and reduced pitching frequency) that are consistent with biological data.
\cite{Quinn2015} conducted a large set of experiments on two-dimensional pitching and heaving flexible plates at various stiffness values, kinematic parameters, and incoming flow speeds. By combining grid search and gradient-based optimization methods, they found that optimizing the pitching angle with heaving can almost double the propulsor efficiency compared to heave-only motions.
\cite{Hoover2018} conducted simulations combining three-dimensional Navier-Stokes equation with one-dimensional Euler–Bernoulli beam theory to analyze the motion of heaving flexible plates, and identified local peaks in swimming speed over a parameter space consisting of the beam material property and heaving frequency.
Also using three-dimensional simulations of pitching plates of uniform flexibility, \cite{Dai2012} found that the phase delay $\phi$ between the leading and trailing edge of the plate decreases with increasing stiffness $\kappa$. For large stiffness, the plate moves in no-neck mode (in-phase in our notation), in which thrust production is close to that of a rigid pitching plate with similar trailing edge displacement.

The effect of uniform flexibility on swimming speed, thrust generation, swimming energetics, and stability have been analyzed in numerous other experimental and computational studies; see, e.g., \cite{Shoele2016,Feilich2015,Michelin2009,Miao2006,Combes2001,Hua2013,Wang2020,Ryu2019}.
Specifically,~\cite{Liu1997,Heathcote2008,Tangorra2010} have indicated that flexibility could lessen or prevent thrust production.
Flexible propulsors have also been widely used in man-made biomimetic underwater autonomous vehicles (UAV); see, e.g., \cite{Fujiwara2017,Katzschmann2018,Gibouin2018,Zhu2019,White2021}. Most notable is the Tunabot design of
\cite{Zhu2019} and \cite{White2021} which mimics the shape and bending kinematics of yellowfin tuna (\textit{Thunnus albacares}) and Atlantic mackerel (\textit{Scomber scombrus}).

While most studies have focused on uniformly flexible bodies, \cite{Combes2003} and \cite{Lucas2014} noted that the stiffness along the fish body is not uniform, but decreases towards the tail, and that the propulsor becomes highly flexible at the flexion point of the body; see Figure~\ref{fig:biodata}(A).
To explore the effects of non-uniform flexibility on efficiency and thrust production, \cite{Lucas2015} considered flexible plates of inhomogeneous stiffness undergoing heaving and pitching motions in a water tunnel and found that non-uniform stiffness can improve thrust production, and that in order to achieve optimal propulsion, the morphologic factor (flexion ratio) and kinematic factor (motion type and motion parameters) should be considered simultaneously.
\cite{Vincent2020a,Vincent2020b} analyzed the effect of non-uniform flexibility on flight performance in the context of a tumbling wing model, and found that wing tip flexibility that follows the empirical rules reported in \cite{Lucas2014} leads to improved flight performance.

In this paper, we use a simplified two-link fish model to analyze the influence of active and passive bending on the swimming speed and efficiency. Two link models are commonly used to study the effect of flexibility on the performance of flapping bodies \citep{Eldredge2010, Wan2012,Li2015}. We solve for fluid-structure interactions in the context of the vortex sheet method as described in Section~\ref{sec:formulation}.
The vortex sheet model has been used extensively to analyze problems of fluid-structure interactions, including ring formation at the edge of a circular tube~\citep{Nitsche1994} and wakes of oscillating plates~\citep{Jones2003,Sheng2012}, falling cards~\citep{Jones2005}, flapping flexible flags~\citep{Alben2008,Alben2009}, swimming plates~\citep{Wu1971}, hovering flyers~\citep{Huang2016, Huang2018} and schooling of swimming plates~\citep{Heydari2020}.
Here, we use the implementation of \cite{Nitsche1994}.
In Section~\ref{sec:result1}, we report the effects of both active and passive bending on the swimming performance of sinusoidally pitching swimmers compared to rigid flapping.
In Section~\ref{sec:con}, we discuss these findings in the context of existing work and highlight the implications of our results on the design of underwater autonomous vehicles.

\section{Problem formulation}\label{sec:formulation}
\begin{figure*}
\centering
\includegraphics[scale=1]{./figures/Fig2.eps}
\caption[]{\footnotesize \textbf{Wake and flow velocity of free swimmers.} (A) Rigid swimmer undergoing periodic pitching (inset) of period $T=1$ and amplitude $a=15^\circ$.  (B,C) Active bending with both anterior and posterior sections undergoing periodic pitching albeit at different amplitude and phase. The anterior follows the same pitching motion (blue line in inset) as the rigid swimmer while the relative rotation of the posterior follows a prescribed Jacobi elliptic sine function (red line in inset) with flexion amplitude $ A = 30 ^ \circ$, flexion ratio $l/L=0.7$, elliptic modulus $M=0.9$, and phase  $\phi= 0$ (in-phase)  and $\phi=\pi$ (antiphase) in B and C, respectively. (D) Passive bending of posterior section while anterior section follows the same prescribed pitching as the rigid swimmer. Joint parameters are set to $\kappa =0$ and $c=1$.
The dissipation time is set to be $T_{\rm{diss}}=1.65T$ in (A-C), $\sqrt{2.09}T$ in (D).}
\label{fig:flowfield}
\end{figure*}

We model the flexible swimmer as a planar two-link body of total length $L$, negligible thickness $e\ll L$, and total mass per unit depth $m=\rho eL$, where $\rho$ is both the swimmer and fluid density, assuming a neutrally-buoyant fish.
The flexion point indicates where the anterior link (of length $l$) is joined to the posterior link; see Figure~\ref{fig:biodata}(A).
The anterior link undergoes sinusoidal pitching motion $\theta_a(t) = a\sin(2\pi f t)$, where $\theta_a$ is the angle relative to the swimming direction, taken to be parallel to the $x$-axis. Here, $a$ is the flapping amplitude, $f={2 \pi}/T$ the flapping frequency, and $T$ the flapping period.
When the posterior is connected rigidly to the anterior link at zero flexion, the posterior motion $\theta_p(t)$ is equal to $\theta_a(t)$ and the flexion angle $\alpha = \theta_p - \theta_a$ is identically zero for all time. The two links form a single rigid plate (Figure~\ref{fig:flowfield}A) whose swimming motion due to sinusoidal pitching has been extensively analyzed~\citep{Jones2003,Sheng2012,Moored2019,Heydari2020,Labasse2020}. To explore the effects of body bending on swimming, we consider two cases: (i)
active bending where the flexion angle $\alpha(t)$ is controlled by the swimmer, and (ii) passive bending where the flexion angle $\alpha(t)$ is dictated by the physics of fluid-structure interactions.

When the two-link swimmer bends actively, we allow the anterior link to have a phase advantage of magnitude $\phi$ relative to the flapping motion of the posterior link. At $\phi = 0$, both anterior and posterior links flap in phase and the swimmer bends in the direction of flapping (Figure~\ref{fig:flowfield}B); for $\phi = \pi$, they flap antiphase resulting in bending in the opposite direction to the anterior pitching motion (Figure~\ref{fig:flowfield}C). The  flexion angle $\alpha(t)$ follows a Jacobi elliptic sine function $\alpha(t)=A \operatorname{sn}(4 K f t, M)$, where $A$ is the maximum flexion angle and $M$ is the elliptic modulus that controls the shape of the elliptic sine function. As $M \to 0$, the elliptic sine function tends to a sinusoidal function and as $M \to 1$, it approaches a square wave shape. The parameter $K$ is introduced to ensure that the flapping frequency of the posterior link is the same as that of the anterior link ($K$ is related to the elliptic modulus $M$ via the elliptic integral and it is given by $K = \mathrm{ellipticF}(\pi / 2, M)$ in Matlab). In this paper, without other specification, we fix $M=0.9$ and explore the effects of the anterior-to-posterior phase difference $\phi$, flexion ratio $l/L$, and maximum flexion angle $A$ on the swimming performance.

We write the equations governing the self-propelled motion of the two-link swimmer in non-dimensional form. To this end, we  scale all parameter values using $L/2$ as the characteristic length scale, $1/f$ as the characteristic time scale, and $\rho (L/2)^2$ as the characteristic mass per unit depth. Accordingly, velocities are scaled by $Lf/2$, forces by $\rho f^2 (L/2)^3$, moments by $\rho f^2 (L/2)^4$, and power by $\rho f^3 (L/2)^4$.
The equation of motion governing the free swimming $x(t)$ is given by Newton's second law
\begin{equation}
m \ddot{x}=  F_x - D_x .
\label{eq:eom}
\end{equation}
Here, $F_x$ and $D_x$ denote the components of the hydrodynamic  pressure and drag forces $F$ and $D$ in the swimming direction. Specifically, $F$ is the total force  normal to the swimmer due to hydrodynamic pressure and $D$ is the force tangential to the swimmer due to skin drag. The hydrodynamic pressure force is calculated in the context of the inviscid vortex sheet model~\citep{Nitsche1994, Huang2016, Huang2018, Heydari2020}, and the drag force $D$ is introduced to emulate the effect of fluid viscosity~\citep{Fang2016,Ramananarivo2016}. A brief overview of the vortex sheet method and its numerical implementation is given in Appendix~\ref{app:VSM} and \ref{app:numerics}. Detailed expressions of the fluid forces acting on the swimmer are given in Appendix~\ref{app:forces}.

When the swimmer bends passively, the relative rotation $\alpha(t)$ of the posterior end is not prescribed a priori and follows from the physics of fluid-structure interactions. Considering that the rotational joint at the flexion point is equipped with a torsional spring of stiffness $\kappa$ and damping coefficient $c$, we write the equation governing the rotational motion of the posterior link
\begin{equation}
I_{p}(\ddot{\theta}_a + \ddot \alpha) +  c \dot{\alpha} + \kappa \alpha =  M_{p} + M_{\rm inertia},
\label{eq:eomalpha}
\end{equation}
where $I_{p}$ and $M_p$ are the moment of inertia and hydrodynamic moment acting on the posterior link about the flexion point, $M_{\rm{inertia}}$ is an inertial moment that arises because the flexion point about which the moments are balanced is moving;  see details in Appendix~\ref{app:forces}.

To assess the swimming performance of the two-link swimmer, we introduce five metrics: the period-average swimming speed $U=\int_t^{t+T}{\dot x dt}$ at steady state, the thrust force $F_x$, the period-average input power $P=\int_t^{t+T}{ P(t) dt}$ required to maintain the prescribed flapping motions, and the propulsion efficiency $mU^2/2PT$ defined as the kinetic energy of the swimmer divided by the input work over one flapping period; see Appendix~\ref{app:forces} for more details.

\begin{figure*}
\centering
\includegraphics[scale=1]{./figures/Fig4.eps}
\caption[]{\footnotesize \textbf{Active anterior-to-posterior bending of free swimmers.} Time-dependent speed, thrust, input power and circulation for (A) rigid swimmer undergoing pitching at $a=15 ^ \circ$, $T=1$, active flexion (B) at phase difference $ \phi=0$, (C) at phase difference $ \phi=\pi$. In (B) and (C), flexion ratio $l/L=0.7$ and flexion angle $A=30 ^ \circ$.
Solid lines represent the instantaneous values and dashed lines represent time-period averages. The results are shown after the swimmers have reached steady state $t_{\rm{s}}=15T$. The dissipation time is set to be $T_{\rm{diss}}=1.65T$.
}
\label{fig:compare}
\end{figure*}
\begin{figure*}
\centering
\includegraphics[scale=1]{./figures/Fig5.eps}
\caption[]{\footnotesize \textbf{Active anterior-to-posterior bending can minimize lateral  forces and negative thrust.} Force hodograph of (A) rigid flapping, (B) in-phase flexion, (C) antiphase flexion for the same cases shown in Figure~\ref{fig:compare}. The arrow indicates the direction of time. The white, blue, red, yellow points are: $t/T= 0$, $0.25$, $0.5$, $0.75$. Bottom row shows force distribution at the instants indicated in red stars.
Blue arrows represent pitching direction of the anterior link motion while red arrows  of the posterior link.
}
\label{fig:hodograph}
\end{figure*}

\section{Results} \label{sec:result1}
We compare the free swimming that results from flapping while undergoing active and  passive bending to that of rigidly flapping. All swimmers have the same total length $L$ and undergo the same sinusoidal pitching motion about their leading edge $\theta_a = a\sin(2\pi f t)$ with $a= 15^o$ and $f=1$, albeit exhibiting distinct bending patterns. Figure~\ref{fig:flowfield} shows snapshots of the wake represented by the free vortex sheet and velocity field generated by a swimmer undergoing (A) rigid flapping, (B) in-phase active bending with flexion amplitude $A = 30^o$ and flexion ratio $l/L = 0.7$, (C) antiphase active bending at the same flexion amplitude and ratio, and (D) passive bending. All snapshots are taken at the same instant in the flapping cycle (at $0.25T$ after steady state has been reached).
Compared to the rigid swimmer, in-phase flexion produces wider wakes and larger leading edge circulation and instantaneous flow speeds, while antiphase flexion is characterized by a leaner, longer wake with weaker leading edge circulation and lower flow speeds.
The main features of the instantaneous flow during antiphase flexion, namely, the leaner wake and weaker leading edge circulation are also observed in passive bending of the swimmer body. These flow features are important for the hydrodynamic forces exerted on the swimmer as discussed later.

\begin{figure*}
\centering
\includegraphics[scale=1]{./figures/Fig6.eps}
\caption[]{\footnotesize \textbf{Active anterior-to-posterior bending of swimmers fixed in oncoming flow.} Time-dependent drag force, input power, and wake circulation for (A) rigid swimmer,  (B) in-phase active flexion ($\phi=0$), (C) antiphase active flexion ($\phi=\pi$). In all cases, the swimmer is fixed in a uniform oncoming flow at $U=9$.
Results are shown after the swimmers have reached steady state $t_{\rm{s}}=11T$. From top to bottom, drag force, input power and circulation of wake are shown. Solid lines represent the instantaneous values and dashed lines represent time-period averages. The dissipation time is set to be $T_{\rm diss}=1.65T$.
}
\label{fig:compare2}
\end{figure*}

We quantitatively evaluate the steady state motion of the rigid and actively-bending swimmers in Figure~\ref{fig:flowfield}(A-C). In Figure~\ref{fig:compare}, from top to bottom, we report the instantaneous (solid lines) and period-average (dashed lines) values of the swimming speed, thrust force, input power, and circulation.
On average, rigid flapping produces the lowest swimming speed while antiphase flexion the highest. Fluctuations around the average swimming speed are smallest for the swimmer undergoing antiphase flexion.
The discrepancy in average swimming speeds between the three flapping modes is surprising at first sight given that the average values of the thrust force are comparable.  However, a closer look at the instantaneous thrust shows that the swimmer undergoing antiphase flexion hardly experiences negative thrust over its flapping cycle. In-phase flexion leads to negative thrust of high magnitudes over larger subintervals of the flapping cycle, as highlighted further in Figure~\ref{fig:hodograph}.
Consequently, the required input power for in-phase flexion is largest compared to both rigid flapping and antiphase flexion.
This is also true of the overall wake circulation.  It is worth noting that, by Kelvin's circulation theorem, circulation around the leading edge must be equal to the overall circulation in the wake. Therefore, compared to rigid flapping, in-phase flexion increases the circulation around the leading edge of the swimmer while antiphase flexion decreases leading edge circulation as noted qualitatively in Figure~\ref{fig:flowfield}.

\begin{figure*}
\centering
\includegraphics[scale=1]{./figures/Fig7_v2.eps}
\caption[]{\footnotesize \textbf{Performance of active anterior-to-posterior bending as a function of phase.}
(A) Free swimmer: swimming speed (black),  power requirement (red), and efficiency (blue) are scaled by the corresponding values of a rigid swimmer.  (B) Swimmer fixed in oncoming flow of uniform speed $U$: drag force (black), power requirement (red) and efficiency (blue) are scaled the corresponding values of a pitching rigid plate fixed oncoming uniform flow. Parameters values are set to $a=15^\circ$, $l/L=0.7$, $A=20^\circ$, $M=0.9$, $U=9$. Flexion agreement parameter $Z$ between the relative velocity of an actively bending posterior and the fluid velocity generated by a passively bending swimmer at zero stiffness $\kappa = 0$ (green) as a function of phase $\phi$ during (A) free swimming  and (B) holding station in oncoming flow $U=5$.
}
\label{fig:mix}
\end{figure*}

The results in Figure~\ref{fig:compare} indicate that the swimmer undergoing antiphase flexion achieves higher swimming speed at lower power requirement and energetic cost. To elucidate the hydrodynamic forces at play, we report in Figure~\ref{fig:hodograph} the force hodograph defined as a plot of the lateral pressure force $F_y$ versus thrust $F_x$ acting on each swimmer.
Snapshots of the distribution of hydrodynamic pressure forces along the swimmers are depicted in bottom row of Figure~\ref{fig:hodograph}, and indicate that the $x$-component of the forces on the anterior and posterior sections during antipase flexion act opposite to each other as in a tug-of-war, leading to overall reduction in thrust values.
Importantly, antiphase flexion also reduces the lateral force and negative thrust, with negative thrust experienced only over a  small subinterval of the flapping period, as noted earlier.
In contrast, in-phase flexion significantly increases the lateral force and negative thrust.

To explain the effect of flexion on the lateral force experienced by the swimmer, it is instructive to re-examine the flow field around the rigid and actively-bending swimmers in Figure~\ref{fig:flowfield}(A-C).
A large leading edge vortex (LEV) is known to generate large lift in flapping flight; see, e.g.,~\cite{Ellington1984,Dickinson1999}. In swimming, larger leading edge circulation creates larger lateral force, which explains why, compared to rigid flapping, in-phase flexion increases the lateral force acting on the swimmer while antiphase flexion decreases it.
Lift is beneficial for flight but large lateral forces are detrimental to swimming speed, as noted in~\cite{Drucker2000} for fish and recapitulated here in the context of our swimmer model.

To further analyze the difference in the swimming performance between rigid flapping and flapping with inphase and antiphase active bending, we fix the swimmer in an oncoming uniform flow of speed $U$  and we compute the hydrodynamic drag forces in each case. In Figure~\ref{fig:compare2}, we report the drag force, input power, and circulation in the wake of the fixed swimmer.
Compared to rigid flapping, antiphase flexion reduces instantaneous drag, power, and circulation, while in-phase flexion increases all three quantities. Reduced drag implies lower thrust requirement during steady state swimming, which provides another perspective for understanding the improved performance of antiphase flexion.

We next examine the period-average performance of actively bending swimmers as a function of anterior-to-posterior phase difference $\phi$. In Figure~\ref{fig:mix}(A), we consider the case of free swimming, we fix the flexion ratio $l/L=0.7$ and flexion amplitude $A=20^\circ$, and we plot the swimming speed $U$, input power $P$ and efficiency $\eta$ versus $\phi$, all scaled by the corresponding values of a rigidly flapping swimmer $U_{\rm rigid}$, $P_{\rm rigid}$, $\eta_{\rm rigid}$, respectively.
We find that active bending is always beneficial in terms of enhanced speed relative to rigid flapping, albeit at an increased power requirement. Importantly, as the anterior-to-posterior bending changes from in-phase flexion to flexion at a phase lag, the scaled speed increases and the scaled power requirement decreases. Optimal performance occurs at $\phi=0.9$ and $\phi=0.8$ in terms of maximum swimming speed and minimum input power and maximum efficiency, respectively. In Figure~\ref{fig:mix}(B), we fix the swimmer in oncoming flow of uniform speed $U$ and compute the scaled drag force $D$, input power $P$ and efficiency $\eta$ as a function of $\phi$ scaled by the corresponding values of a fixed rigid flapper.
Analogous to the free swimmer, drag and input power are minimal at $\phi=0.8$. Taken together, the results imply that antiphase active flexion is near optimal for enhancing speed and efficiency and reducing drag force and power requirement.

To complete this analysis, we also explored the effect of the flapping parameter $M$ on the swimming performance. We found that for a fixed phase, the swimming speed and efficiency change monotonically with $M$ with maximum speed and minimum efficiency as $M\to 1$; see Figure~\ref{fig:appendix_M} in the Appendix. That is, reversing the bending direction with a quick flicker improves speed at the cost of decreasing efficiency.

\begin{figure*}
\centering
\includegraphics[scale=1]{./figures/Fig3.eps}
\caption[]{\footnotesize \textbf{Passive anterior-to-posterior bending of free swimmers.} (A) Top to bottom: passive flexion angle $\alpha$ of posterior end at zero stiffness ($\kappa=0$),  moderate stiffness ($\kappa=20$), large stiffness ($\kappa=40$). The anterior link is pitching at $a=15^\circ$. The damping ratio is set to $c=1$, and flexion ratio to $l/L=0.7$. (B) The bending parameters (phase $ \phi$ and maximum flexion angle $A$). (C) Swimming speed $U$ and propulsion efficiency $\eta$ as a function of stiffness $\kappa$ for three flexion ratios $l/L=0.6$, $0.7$, $0.8$ reported in orange, green and black, respectively. The results are shown after the swimmers have reached steady state $t_{\rm{s}}=25T$. The dissipation time is set to be $T_{\rm{diss}}=\sqrt{2.09}T$.
}
\label{fig:passive}
\end{figure*}

Is active bending necessary for obtaining this enhancement in swimming speed and efficiency over rigid flapping? To address this question, we examine the free swimming of a passively bending swimmer, where the posterior end flaps passively under the effect of hydrodynamic forces and moments generated by the pitching motion of the anterior section. Elastic forces due to a spring of stiffness $\kappa$ located at the flexion point are also at play. In Figure~\ref{fig:passive}(A), we keep all parameter values the same as those used for the actively bending swimmer, and, from top to bottom, we report the flapping motions of the anterior and posterior ends for stiffness values $\kappa=0$, $20$, and $40$. At zero stiffness, flexion introduces no restoring forces and moments. The posterior part rotates antiphase relative to the flapping motion of the anterior part, at an amplitude comparable to the anterior pitching amplitude. The associated wake, shown in Figure~\ref{fig:flowfield}(D) and represented by the free vortex sheet in the inset of Figure~\ref{fig:passive}(A), shares similar features to the wake obtained during antiphase active flexion. At moderate spring stiffness $\kappa=20$, the flexion amplitude increases ($\alpha_{\max} \approx 80^\circ$), and the wake also exhibits larger lateral dispersion.
At large stiffness $\kappa=40$, the  posterior part rotates in-phase with the anterior part at the same flapping amplitude in a way reminiscent to rigid flapping; as reported in \cite{Dai2012} for flexibble pitching plates.

The relative motion of the posterior part is close to a sinusoidal function for all stiffness values $\kappa$. Therefore, for each $\kappa$ value, we fit $\alpha(t)$  by a sine function $\alpha(t)=A \sin (\omega t- \phi)$ using a
standard algorithm~\citep{More1983}. For all fitting, we have at least $95\%$ confidence and frequency $\omega \approx 2\pi$, thus ensuring  convergence of the fitting. In Figure~\ref{fig:passive}(B), we report the fitted flexion amplitude $A$ and phase difference $\phi$ as a function of spring stiffness $\kappa$ for three distinct flexion ratios $l/L =0.6$, $0.7$ and $0.8$. We find that for all $l/L$, as stiffness $\kappa$ increases, the phase difference ${\phi}$ decreases monotonously from $\pi$ to $0$ implying that the posterior flapping motion changes from antiphase to in-phase. The maximum flexion angle ${A}$ first increases with increasing $\kappa$, then decreases to nearly zero at large stiffness implying rigid flapping of both anterior and posterior ends.

\begin{figure*}
\centering
\includegraphics[scale=1]{./figures/Figure6_vXX.eps}
\caption[]{\footnotesize \textbf{Active bending in agreement with passive hydrodynamics  improves swimming performance.}
 Illustration of the interaction with the flow field for (A) in-phase and (B) anti-phase active flexion. Maximum improvement in swimming performance occurs when the tail benefits from flows created by the anterior portion of the fish body.
Light grey arrows represent the flow direction. Red arrows and blue arrows represent the flapping direction of the posterior and anterior sections of the fish body, respectively.
}
\label{fig:cartoon}
\end{figure*}

We compute the associated swimming speed and efficiency for each stiffness value $\kappa$ and we scale the results by those of a rigid swimmer; see Figure~\ref{fig:passive}(C). Clearly, the swimmer with passive flexion never surpasses the swimming speed of a rigid swimmer. At very low stiffness ($\kappa \approx 0$), passive flexion results in a swimming speed close to that of rigid flapping while doubling the swimming efficiency. The increase in swimming efficiency at small stiffness comes purely from a decrease in power requirement compared to rigid flapping. This is in contrast to active flexion where the enhancement in swimming speed and efficiency noted in Figure~\ref{fig:mix} comes at an increase in power requirement relative to rigid flapping.
As $\kappa$ increases, the scaled swimming speed  and propulsion efficiency decrease, indicating that moderate flexibility is detrimental to both speed and efficiency. For large $\kappa$, the speed and efficiency converge to the same speed and efficiency as the rigid swimmer, consistent with the results of~\cite{Dai2012}. We repeat this analysis for three flexion ratios $l/L=0.6$, $0.7$, $0.8$. The scaled speed seems to increase monotonically with increasing $l/L$, but the scaled efficiency seems to peak at $l/L = 0.7$ but only for a range of small $\kappa$ values.
These findings imply that, unlike flapping insect wings~\citep{Ellington1985,Huang2015a}, restoring elastic forces seem to be detrimental to swimming performance. Swimming efficiency peaks at low stiffness values when the restoring spring forces are weak and the posterior end is driven passively by the fluid forces.
\begin{figure*}
\centering
\includegraphics[scale=1]{./figures/Fig8.eps}
\caption[]{\footnotesize \textbf{Swimming performance of actively flexion swimmer scaled by performance of rigid swimmer.} Average speed and efficiency versus flexion ratio $l/L$ and flexion angle $A$ for (A) in-phase flexion ($\phi=0$) and (B) antiphase flexion ($\phi=\pi$). The amplitude for the proximal part and the elliptic modulus are $a=15^\circ$ and $M=0.9$. Dark grey areas indicate regions of diminished performance while light gray areas indicate improved performance over a rigid swimmer. The dissipation time is set to be $T_{\rm{diss}}=1.65T$.
}
\label{fig:colormap}
\end{figure*}

Could the swimmer learn from passive flexion to improve its performance by bending actively in a way that exploits the hydrodynamic forces generated naturally during passive flexion? The results in Figures~\ref{fig:flowfield}-\ref{fig:mix} for actively bending swimmers suggest that maximum benefit occurs for near anti-phase flexion, whereas the results in Figure~\ref{fig:passive} show that maximum efficiency for passively bending swimmers occurs for zero stiffness $(\kappa = 0)$ for which the posterior bends anti-phase. Importantly, the main features of the instantaneous flow during antiphase active flexion (Figure~\ref{fig:flowfield}C) are also observed in passive bending at zero stiffness (Figure~\ref{fig:flowfield}D). We thus posit that active bending is most beneficial when the swimmer actively beats its tail in a direction that takes advantage of the natural flows that arise during passive bending.
To test this hypothesis, we define a flexion agreement parameter $Z$ that aims to relate passive and active bending. Starting from a swimmer bending passively at zero spring stiffness $\kappa=0$ (Figures~\ref{fig:flowfield}D and~\ref{fig:passive}A), we assume a hypothetical posterior that is actively flapping about the flexion point of the swimmer at a relative angle  $\alpha = A \sin(2\pi t - \phi)$. We compute the fluid velocity $\mathbf{u}(s,t)$ induced by the passively-bending swimmer along the hypothetical posterior section that is bending actively (here, $s$ is the distance from the flexion point). The flexion agreement parameter $Z$ is given by
\begin{equation}
Z=\frac 1 T \frac 1 {L-l}\int_0^T{\int_0^{L-l}{\mathbf{u}(s,t) \cdot \mathbf{v}(s,t)ds}dt},
\label{eq:flexion_agreement}
\end{equation}
where $\mathbf{v}(s,t) = s \dot{\alpha} \,\mathbf{n}$ is the relative velocity of the hypothetical actively-flapping posterior. Positive values of the flexion agreement parameter imply a beneficial interaction between  the flow generated during passive flexion and the velocity of the hypothesized posterior during active flexion, whereas negative values indicate a detrimental one.

In Figure~\ref{fig:mix}(A), we set $A$ to be equal to the maximum flexion angle of the passively-bending swimmer, and we vary $\phi$ from $0$ to $\pi$. We find that the agreement parameter $Z$, normalized by its maximum value, is largest for antiphase active flexion and smallest near in-phase active flexion.  This result indicates that antiphase active flexion matches the local flow created during passive flexion best, whereas in-phase active flexion acts opposite to these flows. That is, the swimmer during antiphase flexion can utilize better the flow field generated by the pitching motion of its anterior section, and thus it can achieve higher swimming speed and efficiency compared to in-phase flexion.

To emphasize the effect of the interaction between the flow field and kinematics of active flexion on swimming performance,
we schematically summarize the two cases of in-phase and antiphase flexion in Figure~\ref{fig:cartoon}.
During in-phase active flexion, the flow field is characterized by a strong leading edge vortex around the anterior section of the fish and large lateral forces. Further, the posterior part bends opposite to the local flow field of a passively bending swimmer.
During antiphase active flexion, the leading edge vortex is smaller and it is followed by a counter-rotating vortex around the fish mid-section such that the posterior part is moving in synchrony with the downwash flow induced by this counter-rotating vortex. The flow field is helping the motion of the posterior end. These results indicate that active body deformation that are in agreement with local flows produced during passive deformations are more advantageous for enhancing swimming speeds and efficiencies.

\begin{figure*}
\centering
\includegraphics[scale=1]{./figures/Fig9_v5.eps}
\caption[]{\footnotesize \textbf{Relation to fish swimming behavior. }
(A) Comparison of fish flexion parameters (black dots, cite from \cite{Lucas2014}) and regions of optimal performance predicted by our model (anti-phase active flexion swimmer): pink region corresponds to 200\% enhancement in swimming speed, and green region corresponds to 300 \% enhancement in swimming efficiency, both compared to a swimmer of the same total length rigidly flapping with no flexion. Overlap of the two regions is indicated in beige. The contour grey line encloses a region of 600\% enhancement in efficiency.
}
\label{fig:compareFish}
\end{figure*}

Lastly, we explore the effect of maximum flexion angle $A$ and flexion ratio $l/L$ on the period-average values of the swimming speed and efficiency for both in-phase and antiphase active flexion. Specifically, we examine the range $l/L \in [0.4, 1.0]$ and $A \in [0^\circ, 60^\circ]$ for $\phi = 0$ and $\phi = \pi$.  In Figure~\ref{fig:colormap}, we report the period-average values normalized by the corresponding values for a rigid swimmer with pitching amplitude equal to the anterior part amplitude.
We highlight in light and dark grey respectively the regions in the parameter space where the flexible swimmer outperforms and underperforms the rigid swimmer. The swimmer with in-phase flexion swims slower than the rigid swimmer for small flexion ratios ($l/L < 0.7$) and high flexion amplitudes ($A >40^\circ$), and swims faster than rigid swimmer otherwise. This swimmer, however, is always less efficient than the rigid swimmer for reasons explained previously. In Figure~\ref{fig:colormap}(B), for most parameter values, the antiphase swimmer outperforms the rigid swimmer in terms of swimming speed and efficiency. Note that the region with the highest swimming speed advantage lies in $l/L \in [0.6, 0.8]$ and $ A \in [30^\circ,  60^\circ]$, and the region with the highest efficiency advantage lies in $l/L \in [0.6, 0.8]$ and $ A \in [25^\circ, 50^\circ]$.

We compare the regions of highest swimming speed and efficiency obtained during antiphase flexion in Figure~\ref{fig:colormap}(B) to the design parameters of biological fish reported in Figure~\ref{fig:biodata}. In Figure~\ref{fig:compareFish}, we plot the flexion angle $A$ as a function of the flexion ratio $l/L$ for the fish data in Figure~\ref{fig:biodata} and we superimpose on this design space the regions of 200\% enhancement in speed and 300\% enhancement in efficiency from Figure~\ref{fig:colormap}(B).
As shown in Figure~\ref{fig:compareFish}, there is significant overlap between these regions of improved performance and the biological data. Indeed, all biological data lie within the region of improved efficiency.

Many of these fish are known to exhibit migratory behavior that requires efficient swimming. Even baby clownfish are reported to migrate over long distances~\citep{Simpson2014}. Tuna can cover $7600$ km  in one traveling phase~\citep{Itoh2003} and tiger shark are capable of traveling long distances in short time~\citep{Simpfendorfer2009}.  Our results in Figure~\ref{fig:compareFish} are consistent with these facts: tiger shark and clownfish lie in the intersection region of improved speed and efficiency, and tuna lies in the region with highest increase in efficiency. On the hand, butterflyfish, who are only known to migrate over short distances during spawning \citep{Yabuta1997}, and Koi, for which there is no evidence of migration, lie in the regions characterized by smaller increase in speed and efficiency.

\section{Conclusion}\label{sec:con}

We analyzed the swimming performance of flapping swimmers undergoing active and passive deformations. Whereas fish exhibit a variety of swimming modes~\cite{Sfakiotakis1999}, we simplified body deformations to account for only anterior-to-posterior bending, with one degree of freedom describing the relative rotation between the two sections. We explored the effects of morphological and kinematic parameters on the swimming speed and efficiency.

We found that passive body bending, at negligible body stiffness and minor elastic forces, caused anterior-to-posterior antiphase flexion. This antiphase flexion is dictated by the flow physics and causes the swimmer's morphology to get more streamlined compared to rigid flapping with no flexion, thus creating leaner wakes that reduce drag and power requirement and increase efficiency. While drag reduction is desirable for improved efficiency, passive bending also reduced thrust production, thus diminishing swimming speed.  Interestingly, restoring elastic forces at moderate stiffness values seemed detrimental to both swimming speed and efficiency. These findings are consistent with the hypothesis that for maximum efficiency, the fish tail and posterior body should flex like water, exhibiting little or no resistance to the flows generated by the flapping motion of the anterior portion of the body. This hypothesis could explain how anesthetized fish, with no muscle activity, placed in periodic wakes generate oscillatory body deformations that allow the fish to swim upstream~\citep{Beal2006}.

We also found that a swimmer that actively creates antiphase anterior-to-posterior bending enjoyed the same benefits of leaner wake and reduced drag as a passively bending swimmer while mitigating the reduction in thrust and swimming speed.
To quantify the hydrodynamic mechanisms leading to this improved swimming performance, we introduced a flexion agreement parameter that compares the active flexion velocity of the swimmer's posterior to the local flow velocity during passive flexion. We found that during active in-phase flexion, the posterior beats opposite to the local flow that would naturally arise during passive flexion, leading to a negative flexion agreement parameter, and thus lower swimming speed and efficiency.
During active antiphase flexion, the posterior flaps consistently with the local flow, leading to a positive flexion agreement parameter, and improved speed and efficiency.

These findings suggest tremendous versatility in swimming performance, even when accounting only for coarse anterior-to-posterior bending motions. They indicate that fish can readily and fluidly transition from efficient (passive bending) to fast (active bending) by actively beating their tail in agreement with the local flow generated during passive bending.

To explore the role of flow physics on the convergent bending rules of \citep{Lucas2014}, we examined the effect of flexion ratio and maximum flexion angle on the performance of swimmers undergoing active antiphase anterior-to-posterior bending.  We found an optimal region in this design space that simultaneously enhances swimming speed and efficiency. Importantly, we found that this region has a significant overlap with the fish bending parameters reported in \citep{Lucas2014}; see Figure~\ref{fig:cartoon}(A).  In fact, all biological data fell in the region of enhanced efficiency predicted by our model.

Taken together, our results have two major implications on understanding the role of body bending in fish swimming. They are consistent with the hypothesis that fish that actively bend their bodies in a fashion that exploits the local hydrodynamics can at once improve speed and efficiency. They also support the hypothesis that flow physics could have provided a selective force for driving the evolution of fish bending patterns.

Beyond the trade-offs between speed and efficiency explored here, our model could be generalized in future studies to examine transient swimming maneuvers such as turning. \cite{Pollard2019} proposed that a passive posterior not only improves efficiency, but also improves fish maneuverability compared to a rigid body. \cite{Drucker2000} pointed out that although large lateral forces are detrimental to fish swimming speed, they can improve fish maneuverability. Our model predicts large lateral forces during in-phase active flexion, suggesting that this bending pattern, while not optimal for forward swimming, could be beneficial for turning motions. These considerations, as well as models of higher fidelity to the fish biomechanics and fluid environment, will be explored in future research.

Finally, our finding that active body bending and tailbeat patterns that match local flow velocities that would be produced naturally by anterior sections of the fish body could lead to improved performance and energy savings might have important implications on understanding the mechanisms driving body and caudal fin deformations in swimming~\citep{Ramakrishnan2011,Bozkurttas2007}, schooling~\citep{Heydari2020,Becker2015,Li2020}, and navigating ambient unsteady flows~\citep{Liao2003}.


\section*{Introduction}

In the study of linear ordinary differential equations (LODEs), one of the central questions is to find closed form solutions. Here, we give a solution to this question that can be algorithmically implemented for algebraic LODEs, i.e. LODEs whose full system of solutions can be expressed in terms of algebraic functions. We base our approach on the Picard-Vessiot theory which is a galoisian perspective developed by E. Picard \cite{PICARD1991III}, E. Vessiot \cite{VESSIOT1892} and E. R. Kolchin \cite{KOLCHIN1946}.

In the literature there are three different methods for describing the solutions of algebraic LODEs.

In \cite{SINGER1993}, M. Singer and F. Ulmer showed how to compute the minimal polynomial of the solutions of algebraic LODEs. The drawback of this method is that the degree of the minimal polynomial is relatively big, it can reach $120$ for second order equations and $216$ for third order.

A second approach has been given by J. Kovacic in \cite{KOVACIC1986}. J. Kovacic showed that the solutions of second order algebraic LODEs can be expressed as exponential of an integral of an algebraic function, where the degree of the minimal polynomial of the algebraic function is at most $12$. M. van Hoeij, J. F. Ragot, M. Singer, F. Ulmer and J.-A. Weil \cite{VANHOEIJ1999} extended this method to higher orders and showed that, for equations of third order, the minimal polynomial of the algebraic function involved in the integral of the exponential is of degree up to $45$.

A third and the most effective method for expressing the solutions of algebraic LODEs is based on a classical result of F. Klein \cite{KLEIN187711,KLEIN187712}. This result states that a full system of solutions to an algebraic second order LODE with rational coefficients can be given in closed form as solutions to hypergeometric equations precomposed with a rational function, multiplied by the solution to a first order LODE. Furthermore, for each finite group $G\subseteq SL_2(\mathbb{C})$, there exist a standard hypergeometric equation which is minimal in the sense that any other LODE with differential Galois group $G$ is projectively equivalent to a pullback by a rational map of this hypergeometric equation. B. Dwork and F. Baldassari \cite{BALDASSARRI1980,BALDASSARRI1979} presented Klein's result using modern tools and generalized it to second order algebraic LODEs defined over algebraic curves. Moreover, J.-A. Weil and M. van Hoeij \cite{VANHOEIJ2005} identified a standard equation for each of the primitive finite group in $SL_2(\mathbb{C})$ and implemented Klein's result in an algorithm. By introducing the concept of standard equation, M. Berkenbosch \cite{BERKENBOSCH2006} generalized Klein's result to order three. C. Sanabria  \cite{SANABRIA2017} extended Klein-Dwork-Baldassari's result to any order showing that any algebraic LODE is projectively equivalent to the pullback by a rational map of a standard equation. A drawback of this generalization is that the family of standard equations involved is infinite and continuous, making it unsuitable for an algorithmic implementation.

In this paper, we improve the result of C. Sanabria \cite{SANABRIA2017} by proving that for each finite primitive group $G\subseteq SL_n(\mathbb{C})$, there exist a positive integer $d$ and a standard equation such that the solutions to any LODE with differential Galois group isomorphic to $G$ can be expressed in terms of the solutions of this standard equation, together with its derivatives, functions in a field extension of degree $d$, and a solution to a first order LODE defined over this extension. With this result, we reduce the infinite family of standard equations to a singleton. In the last section, we show how to obtain closed form solutions to irreducible algebraic LODEs of order three in terms of generalized hypergeometric functions. The computations were carried out using MAPLE and can be obtained from my webpage \cite{MapleComputations2021}.

\section{Preliminary}

\subsection{Space of orbits}\label{orbitspace}

Let $G\subseteq GL_n(\mathbb{C})$ be a finite linear algebraic group. We define a left group action of $G$ on the coordinate ring $R=\mathbb{C}[X_1,\ldots,X_n]$ of $\mathbb{C}^n$ by $\mathbb{C}$-automorphisms
\[
X_j\mapsto \sum_{i=1}^n X_ig_{ij},
\]
for $(g_{ij})_{i,j=1}^n\in G$.
This $G$-action on $R$ induces a right $G$-action on $\mathbb{C}^n$ defined by
\[
(x_1,\ldots,x_n)\mapsto \Big(\sum_{i=1}^nx_ig_{i1},\ldots,\sum_{i=1}^nx_ig_{in}\Big).
\]
Since $G$ is finite, it is reductive and thus the $G$-invariant polynomials in $R$ separate the $G$-orbits in $\mathbb{C}^n$. Therefore, the coordinate ring of the orbit space $\mathbb{C}^n/G$ is the finitely generated $G$-invariant subring $R^G$.

Let $F_1,\ldots,F_N$ be homogeneous generators of $R^G$. We can embed $\mathbb{C}^n/G$ into $\mathbb{C}^N$ through the algebraic map
\[
\mathbf{x}\cdot G \mapsto \Big(F_1(\mathbf{x}),\ldots,F_N(\mathbf{x})\Big),
\]
where $\mathbf{x}=(x_1,\ldots,x_n)$. We will identify the orbit space $\mathbb{C}^n/G$ with its embedding in $\mathbb{C}^N$. Since $G$ is finite, the orbit space has dimension $n$ and therefore, there exits a dense subset of $\mathbb{C}^n$ where the derivative of the quotient map
\[
\mathbf{x}\mapsto \Big(F_1(\mathbf{x}),\ldots,F_N(\mathbf{x})\Big)
\]
is non-singular.

Similarly, these $G$-actions define $G$-actions on $\mathbb{P}^{n-1}(\mathbb{C})$ and on its homogeneous coordinate ring $R$. The orbit space of the action on $\mathbb{P}^{n-1}(\mathbb{C})$ is $\mathbb{P}^{n-1}(\mathbb{C})/G$. Let $\Lambda\in\mathbb{Z}_{>0}$ be such that the homogeneous elements of $R^G$ of degree $\Lambda$, $R^G_\Lambda$, form a homogeneous coordinate system for $\mathbb{P}^{n-1}(\mathbb{C})/G$, i.e. $$\mathbb{P}(R^G_\Lambda)\simeq\mathbb{P}^{n-1}(\mathbb{C})/G.$$
Let $\Phi_1,\ldots,\Phi_M$ be a basis of $R^G_\Lambda$. We can embed $\mathbb{P}^{n-1}(\mathbb{C})/G$ into $\mathbb{P}^{M-1}(\mathbb{C})$ with the map
\[
[\mathbf{x}]\cdot G\mapsto \Big[\Phi_1(\mathbf{x}):\ldots:\Phi_M(\mathbf{x})\Big],
\]
where $[\mathbf{x}]=[x_1:\ldots:x_n]$. As before, we identify the orbit space $\mathbb{P}^{n-1}(\mathbb{C})/G$ with its embedding in $\mathbb{P}^{M-1}(\mathbb{C})$. Again, there exits a dense subset of $\mathbb{P}^{n-1}(\mathbb{C})$ where the derivative of the quotient map
\[
[\mathbf{x}]\mapsto \Big[\Phi_1(\mathbf{x}):\ldots:\Phi_M(\mathbf{x})\Big]
\]
is non-singular.

\subsection{Schwarz maps}

Let $C_0$ be a compact Riemann surface and let $K=\mathbb{C}(C_0)$ be the field of meromorphic functions over $C_0$. Let $\delta:K\rightarrow K$ be any non-trivial derivation. Note that $\delta$ can be uniquely extended to the sheaf of meromorphic functions over any open set in $C_0$.

Let $L(y)=0$ be an LODE of order $n$ with rational coefficients
\[
L(y)=\delta^n(y)+a_{n-1}\delta^{n-1}(y)+\ldots+a_1\delta(y)+a_0y,
\]
where $a_0,a_1,\ldots,a_{n-1}\in K$. Let $S\subset C_0$ be the collection of singularities of $L(y)$. Given a non-singular point $p\in C_0$ together with a fundamental system of solutions $\mathbf{y}=(y_1,\ldots,y_n)$ over a neighborhood $U\subseteq C_0$ of $p$, we define the Schwarz map as the analytic extension of
\begin{eqnarray*}
[\mathbf{y}]: U & \longrightarrow & \mathbb{P}^{n-1}(\mathbb{C})\\
z & \longmapsto & \Big[y_1(z):\ldots:y_n(z)\Big].
\end{eqnarray*}

The monodromy of the Schwarz map is the projection of the monodromy group of the fundamental system of solutions $(y_1,\ldots,y_n)$ on $PGL_n(\mathbb{C})$. We will denote it by $G_0$ and call it the projective monodromy of $L(y)=0$. Let us assume $G_0$ is finite. Then, post-composing the Schwarz map with the quotient by the action of $G_0$ on  $\mathbb{P}^{n-1}(\mathbb{C})$, we obtain a single-valued map $$\psi:C_0\setminus S\rightarrow\mathbb{P}^{n-1}(\mathbb{C})/G_0,$$ which we will call the quotient Schwarz map.

\subsection{Algebraic Picard-Vessiot theory}

We will briefly recall the results from Picard-Vessiot theory that we will use in this paper. An extensive exposition of this subject is given in \cite{VANDERPUT2003}.

We will assume that $L(y)=0$ is irreducible and that all its solutions are algebraic over $K$. Under these assumptions, if $\mathbf{y}=(y_1,\ldots,y_n)$ is a fundamental system of solutions, then a Picard-Vessiot extension for $L(y)=0$ is the field
$$E=K[y_1,\ldots,y_n].$$
Let us denote by $I$ the kernel of the $K$-morphism
\begin{eqnarray*}
\Phi: K[Y_1,\ldots,Y_n] & \longrightarrow & K\\
Y_j & \longmapsto & y_j.
\end{eqnarray*}
An isomorphic representation of the differential Galois group of $L(y)=0$ is the finite group $G\subset GL_n(\mathbb{C})$ composed by the elements $(g_{ij})_{i,j=1}^n$ that send the ideal $I$ into itself under the $K$-automorphism
\begin{eqnarray*}
K[Y_1,\ldots,Y_n]  & \longrightarrow & K[Y_1,\ldots,Y_n]\\
Y_j & \longmapsto & \sum_{i=1}^nY_ig_{ij}.
\end{eqnarray*}
We will call $G$ the representation induced by $\mathbf{y}$. By the Galois correspondence, if $P\in K[Y_1,\ldots,Y_n]$ is $G$-invariant, then $\Phi(P)\in K$. Since $L(x)=0$ is irreducible, then $G$ is reductive and we have the following theorem.

\begin{thm}\label{compointthm}[Compoint's theorem for the algebraic case \cite{SANABRIA2017,VANDERPUT2020}] Let $L(y)=0$ be an algebraic irreducible LODE. If the differential Galois group of $L(y)=0$ is finite, then the ideal $I$ is generated by the $G$-invariants contained in it. In particular, if $P_1,\ldots,P_N\in\mathbb{C}[Y_1,\ldots,Y_n]$ is a set of generators of the $G$-invariant subring $\mathbb{C}[Y_1,\ldots,Y_n]^G$, then
\begin{equation*}\label{Igen}
I=\langle P_1-f_1,\ldots, P_N-f_N\rangle,
\end{equation*}
where $f_i=\Phi(P_i)$, $i=1,\ldots,N$.
\end{thm}

\begin{rem}
Given $L(y)=0$, as in the previous theorem, and homogeneous generators $P_1,\ldots, P_N$ of $\mathbb{C}[Y_1,\ldots,Y_n]^G$, the values $f_1,\ldots,f_n$ can obtained using the algorithm in \cite{VANHOEIJ1997}.
\end{rem}

\subsection{Gauge equivalence and invariants}

\begin{defn}
Let $L(y)=0$ and $L'(x)=0$ be LODEs of order $n$ with coefficients in $K$, let $K\subseteq K_1$ be a field extension, and let $\delta:K\rightarrow K$ be a non-trivial derivation. We say that $L'(x)=0$ is gauge equivalent to $L(y)=0$ over $K_1$ if there exist $f_0,f_1,\ldots,f_{n-1}\in K_1$ such that $x=f_0y+f_1\delta(y)+\ldots+f_{n-1}\delta^{n-1}(y)$ is a solution of $L'(x)=0$ whenever $y$ is a solution of $L(y)=0$.
\end{defn}

Let $V$ be a  $K$-vector space and let $\delta:K\rightarrow K$ be a non-trivial derivation. A connection $\nabla: V\rightarrow V$ over $(K,\delta)$ is an additive map that satisfies the Leibniz rule $\nabla fv=\delta(f)v+f\nabla v$. A differential module over $(K,\delta)$ is a $K$-vector space endowed with a connection. If $TV$ is a tensorial construction of $V$, then $TV$ inherits a connection $T\nabla$ and therefore, a differential module structure.

The LODE $L(y)=0$ defines the connection $\nabla_L: K^n \rightarrow K^n$ such that
$$\nabla_L(e_i) =  e_{i+1},\quad i=1,\ldots,n-1$$
and
$$\nabla_L(e_n) = - a_{n-1}e_n-\ldots-a_1e_2-a_0e_1,$$
for the canonical basis $(e_1,\ldots,e_n)$.
There is a well know correspondence between the homogeneous invariants of degree $d$ in $K[Y_1,\ldots,Y_n]^G_d$ and the solutions to $S^d\nabla_L v=0$, where $S^d$ is the $d$-th symmetric power.

Let $K\subseteq K_1$ be a field extension. Let $f_0,f_1,\ldots,f_{n-1}\in K_1$ and let
$$L'(x)=\delta^n(x)+b_{n-1}\delta^{n-1}(x)+\ldots+b_1\delta(x)+a_0x=0$$
be the LODE such that $x=f_0y+f_1\delta(y)+\ldots+f_{n-1}\delta^{n-1}(y)$ is a solution of $L'(x)=0$ whenever $y$ is a solution of $L(y)=0$. Then, if $(e_1,\ldots,e_n)$ is the canonical basis of $K_1^n$, $\nabla_L: K_1^n \rightarrow K_1^n$ is the connection defined by $L(y)=0$, and if we define
$$e'_1 = f_0e_1+f_1e_2+\ldots+f_{n-1}e_n$$
and
$$e'_i=\nabla_Le'_{i-1},\quad i=2,\ldots,n,$$
then
$$\nabla_L(e'_n)=- b_{n-1}e'_n-\ldots-b_1e'_2-b_0e'_1.$$
In particular, if $x_i=f_0y_i+f_1\delta(y_i)+\ldots+f_{n-1}\delta^{n-1}(y_i)$, $i=1,\ldots,n$, then the correspondence between homogeneous invariants and solutions to the symmetric power \cite{SANABRIA2014} implies that for every $P\in K[Y_1,\ldots,Y_n]^G_d$, $P(x_1,\ldots,x_n)$ can be written as a linear combination of elements in $\mathbb{C}[y_1,\ldots,y_n]^G_d$ with coefficients in $K_1$. We obtain the following result.

\begin{prop}\label{propgau}
Let $L(y)=0$ and $L'(x)=0$ be LODEs of order $n$ with coefficients in $K$ that are gauge equivalent over $K_1\supset K$. Let $\delta:K\rightarrow K$ be a non-trivial derivation. Then we have the following.
\begin{itemize}
\item[i)] The differential Galois groups of $L(y)=0$ and $L'(x)=0$ over $K_1$ are isomorphic.
\item[ii)] If $\mathbf{y}=(y_1,\ldots,y_n)$ and $\mathbf{x}=(x_1,\ldots,x_n)$ are fundamental systems of solutions of $L(y)=0$ and $L'(x)=0$, respectively, such that they induce the same representation $G\subset GL_n(\mathbb{C})$ of their differential Galois groups, then, given $P\in K[Y_1,\ldots,Y_n]^G_d$ and a basis $Q_1,\ldots,Q_m$ of $\mathbb{C}[Y_1,\ldots,Y_n]^G_d$, there exist $g_1,\ldots,g_m\in K_1$ such that
$$P(\mathbf{x})=g_1Q_1(\mathbf{y})+\ldots+g_mQ_m(\mathbf{y}).$$
\item[iii)] If $f_0,f_1,\ldots,f_{n-1}\in K_1$ are such that
$$x_i=f_0y_i+f_1\delta(y_i)+\ldots+f_{n-1}\delta^{n-1}(y_i),\ i=1,\ldots,n,$$
then $g_i\in \mathbb{C}[f_0,\ldots,f_{n-1}]^G_d$.
\end{itemize}
\end{prop}

\subsection{Schwarz maps for invariant curves}

Let $G$ be a finite primitive algebraic subgroup of $GL_n(\mathbb{C})$ and let $C\subseteq\mathbb{P}^{n-1}(\mathbb{C})$ be an algebraic $G$-invariant curve not contained in a projective hyperplane. Then, by the following theorem \cite{SANABRIA2017arxiv}, there exist an LODE $L_C(y)=0$ such that the image of its Schwarz is $C$. Furthermore, the field of definition of this equation is an abelian extension of the field of meromorphic functions of $C$.

\begin{thm}\label{theoalg}
Let $G$ be a finite primitive algebraic subgroup of $GL_n(\mathbb{C})$ and let $C\subseteq\mathbb{P}^{n-1}(\mathbb{C})$ be an algebraic $G$-invariant curve not contained in a projective hyperplane. Let $C_0$ be a compact Riemann surface  such that there is a non-constant morphism of curves
$$\psi: C_0\setminus S\rightarrow C/G\subseteq\mathbb{P}^{n-1}(\mathbb{C})/G,$$
with $S\subseteq C_0$ finite, and such that $C\rightarrow C/G$ is unramified over $\psi(C_0\setminus S)$ (i.e. $C\rightarrow C/G$ is a local isomorphism of complex curves over every point  of $\psi(C_0\setminus S)$). Let $K_0=\mathbb{C}(C_0)$ be the field of meromorphic functions over $C_0$. Then there exists a branched cover $\pi:C_1\rightarrow C_0$, where $C_1$ is a compact Riemann surface and $K_0\subseteq K_1=\mathbb{C}(C_1)$ is an abelian extension, and there exists a unique monic linear differential equation $L_C(y)=0$ of order $n$ with coefficients in $K_1$ admitting a fundamental system of solutions $(y_1,\ldots,y_n)$ over an open set $U\subseteq C_1$, such that the closure of the image of its Schwarz map is $C$ (see diagram).
Moreover, the restriction of $\psi\circ\pi$ to $C_1\setminus\pi^{-1}(S)$ is the quotient Schwarz map of $L(y)=0$ associated to the system of solutions $(y_1,\ldots,y_n)$.
\end{thm}

{\center
\hfill
\xymatrix{
 & & & C\ar[ddd]^\Pi\\
 & & & \\
C_1\ar@{..>}[rrruu]^{[\mathbf{y}]}\ar[d]_\pi &  &  & \\
C_0\ar@{..>}[rrr]_\psi^\simeq & & &  C/G
}
\hfill
}

\section{Invariant curves, pullbacks and gauge equivalence}

\subsection{Invariant curves and pullbacks}

Let $G$ be a finite primitive algebraic subgroup of $GL_n(\mathbb{C})$ and let $C\subseteq\mathbb{P}^{n-1}(\mathbb{C})$ be an algebraic $G$-invariant curve not contained in a projective hyperplane. Let $L_C(y)=0$ be as in Theorem \ref{theoalg}. Suppose now that we have another $n$-th order LODE $L(x)=0$ admitting a fundamental system of solutions $(x_1,\ldots,x_n)$ such that the image of the associated Schwarz map is $C$. We will prove that $L(x)=0$ is projectively equivalent to the pullback of $L_C(y)=0$ by a map in an abelian extension of the field of definition $L$.

\begin{thm}\label{theopullback}
Let $G$ be a finite primitive algebraic subgroup of $GL_n(\mathbb{C})$ and let $C\subseteq\mathbb{P}^{n-1}(\mathbb{C})$ be an algebraic $G$-invariant curve not contained in a projective hyperplane. Let $C_0$, $C_1$, and $L_C(y)=0$ be as in Theorem \ref{theoalg}. Let $X_0$ be a compact Riemann surface, let $F_0=\mathbb{C}(X_0)$ be the field of meromorphic functions over $X_0$, and let $L(x)=0$ be an LODE with coefficients in $F_0$ such that, for a fundamental system of solutions, the closure of the image of its associated Schwarz map is $C$. Then there exists a branched cover $p:X_1\rightarrow X_0$, where $X_1$ is a compact Riemann surface and $F_0\subseteq F_1=\mathbb{C}(X_1)$ is an abelian extension, and there exist two functions $f,h\in F_1$ and $\Lambda\in\mathbb{Z}_{>0}$ such that $(x_1,\ldots,x_n)$, with $$x_i=f^{1/\Lambda}\cdot y_i\circ h,\quad i=1,\ldots,n,$$ is a fundamental system of solutions of $L(x)=0$ over an open set $V\subseteq C_1$.
\end{thm}

\begin{proof}
Let $\mathbf{z}=(z_1,\ldots,z_n)$ be a fundamental system of solutions of $L(x)=0$ over an open set $U_0\subseteq X_0$ such that the image of its Schwarz map is $C$. We may assume that the domain $U$ of $\mathbf{y}=(y_1:\ldots:y_n)$ and $U_0$ are such that the image of $[\mathbf{z}]=[z_1:\ldots:z_n]$ is included in the image of $[\mathbf{y}]=[y_1:\ldots:y_n]$. Let $\Pi: \mathbb{P}^{(n-1)}(\mathbb{C})\rightarrow \mathbb{P}^{(n-1)}(\mathbb{C})/G$ be the projection onto the orbit space so that $\Pi\circ[z_1:\ldots:z_n]$ maps $U_0$ into $C/G$.

Let $\psi$ and $\pi$ be as in Theorem \ref{theoalg} and let $\phi$ be an inverse of $\psi$. Then $\phi\circ\Pi\circ[\mathbf{z}]$ defines an algebraic map $h_0:X_0\rightarrow C_0$. We will identify $K_0=\mathbb{C}(C_0)$ with $h_0^*(K_0)\subseteq F_0$. Let $F_1=F_0[K_1]$, where $K_1=\mathbb{C}(C_1)$. Then $F_0\subseteq F_1$ is an abelian extension. Let $X_1$ be the Riemann surface such that $\mathbb{C}(X_1)=F_1$ and let $p:X_1\rightarrow X_0$ and $h:X_1\rightarrow C_1$ be the morphisms induced, respectively, by the inclusions $F_0\subseteq F_1$ and $K_1\subseteq F_1$. We may assume that $U_0=p(V)$ and $U=h(V)$, for some open set $V\subseteq C_1$.

By definition, we have $h_0\circ p=\pi \circ h$. Furthermore, $h_0=\phi\circ\Pi\circ[\mathbf{z}]$ over $U_0$. On the other hand, we have $\Pi\circ [\mathbf{y}] = \psi\circ\pi$ over $U$. Therefore,
\begin{align*}
\Pi\circ  [\mathbf{y}]\circ h &=\psi\circ\pi\circ h=\psi\circ h_0\circ p\\
 &=\psi\circ \phi \circ\Pi\circ[\mathbf{z}]\circ p\\
 & =\Pi\circ[\mathbf{z}]\circ p
\end{align*}
over $V$. Then there exist $(g_{ij})_{i,j=1}^n\in G$ such that $(g_{ij}):  [\mathbf{y}]\circ h\mapsto [\mathbf{z}]\circ p$. Since the image of $[\mathbf{z}]$ is included in the image of $[\mathbf{y}]$, we have that $(g_{ij})$ is the identity and $[\mathbf{y}]\circ h=[z]\circ p$.

As in Section \ref{orbitspace}, let $\Lambda\in\mathbb{Z}_{>0}$ be such that the homogeneous elements of $\mathbb{C}[X_1,\ldots,X_n]^G$ of degree $\Lambda$ form a homogeneous coordinate system for $\mathbb{P}^{n-1}(\mathbb{C})/G$ and let $\Phi_1,\ldots,\Phi_M$ be a basis of $\mathbb{C}[X_1,\ldots,X_n]^G_\Lambda$. Then, since $\Pi\circ  [\mathbf{y}]\circ h=\Pi\circ[\mathbf{z}]\circ p$, we have
\[
\Big[\Phi_1(\mathbf{z}\circ p):\ldots:\Phi_M(\mathbf{z}\circ p)\Big]=\Big[\Phi_1(\mathbf{y}\circ h):\ldots:\Phi_M(\mathbf{y}\circ h)\Big].
\]
Hence, there exists $f\in F_1$ such that $\Phi_i(\mathbf{z}\circ p)=f\cdot\Phi_i(\mathbf{y}\circ h)$, $i=1\ldots,M$. Since $\Pi$ is locally a biholomorphism, we have $\mathbf{z}\circ p=f^{1/\Lambda}\cdot\mathbf{y}\circ h$. The theorem follows by taking $x_i=z_i\circ p$, $i=1,\ldots,n$.

{\center
\hfill
\xymatrix{
 & & & & C\ar[dddd]^\Pi\\
X_1\ar@{..>}[rrrru]^{[\mathbf{x}]}\ar[dd]_p\ar[rrd]^h & & & & \\
 & & C_1\ar@{..>}[rruu]_{[\mathbf{y}]}\ar[dd]_\pi &  & \\
X_0\ar@{..>}[rrrruuu]^{[\mathbf{z}]}\ar[rrd]^{h_0} & & & &\\
& & C_0\ar@{..>}@<0.7ex>[rr]^\psi & &  C/G\ar@{..>}@<0.7ex>[ll]^{\phi}
}
\hfill
}

\end{proof}

\subsection{Invariant curves and gauge equivalence}

Let $G$ be a finite primitive algebraic subgroup of $GL_n(\mathbb{C})$ and let $C$ and $X$ be algebraic $G$-invariant curves in $\mathbb{P}^{n-1}(\mathbb{C})$, none contained in a projective hyperplane. Let $L_C(y)=0$ and $L_{X}(x)=0$ be LODEs such that for two fundamental systems of solutions $\mathbf{y}=(y_1,\ldots,y_n)$ and $\mathbf{x}=(x_1,\ldots,x_n)$  the images of their Schwarz maps are $C$ and $X$, respectively, as in Theorem \ref{theoalg}. We will prove that if $d\in\mathbb{Z}_{>0}$ is the degree of $C$, then there exists a field extension of the field of definition of $L_{X}$, of degree at most $d$, over which $L_{X}(x)=0$ is gauge equivalent to a pullback of $L_C(y)=0$.

\begin{thm}\label{thetheo}
Let $G$ be a finite primitive algebraic subgroup of $GL_n(\mathbb{C})$ and let $C\subseteq\mathbb{P}^{n-1}(\mathbb{C})$ be an algebraic $G$-invariant curve of degree $d$ not contained in a projective hyperplane. Let $L_C(y)=0$ be an LODE such that for a fundamental system of solutions $(y_1,\ldots,y_n)$ the image of its Schwarz map is $C$, as in Theorem \ref{theoalg}. Let $K$ be the field of definition of $L_C$.

Let $L(x)=0$ be an LODE such that for a fundamental system of solutions $(x_1,\ldots,x_n)$ the image of its Schwarz map is an algebraic $G$-invariant curve $X\subseteq\mathbb{P}^{n-1}(\mathbb{C})$ not contained in a projective hyperplane. Let $F_0$ be the field of definition of $L_{X}$. Then there exist a field extension $F_0\subseteq E$, of degree at most $d$, $n-2$ functions $f_1,\ldots,f_{n-2}\in E$, an abelian extension $E\subseteq E_1$, two functions $f,h\in E_1$, and $\Lambda\in\mathbb{Z}_{>0}$ such that
$$x_i+f_1\delta_0(x_i)+\ldots+f_{n-2}\delta_0^{n-2}(x_i)=f^{1/\Lambda}\cdot y_i\circ h,$$
where $\delta_0:F_0\rightarrow F_0$ is a non-trivial derivation.
\end{thm}

\begin{proof}
We may assume that $y_1,\ldots,y_n$ are functions defined over an open set $U_0\subseteq C_0$, where $C_0$ is a Riemann surface such that $\mathbb{C}(C_0)=K$. Similarly, we may assume that $x_1,\ldots,x_n$ are functions defined over an open set $V_0\subseteq X_0$, where $X_0$ is a Riemann surface such that $\mathbb{C}(X_0)=F_0$.

Since $X$ is not contained in a proyective hyperplane, for a generic $t\in V_0$, we have that the vectors
\begin{align*}
\mathbf{x}(t) &=\left(x_1(t),\ldots,x_n(t)\right),\\
\delta_0\left(\mathbf{x}\right)(t) &=\left(\delta_0\left(x_1\right)(t),\ldots,\delta_0\left(x_n\right)(t)\right),\\
\vdots & \\
\delta_0^{n-2}\left(\mathbf{x}\right)(t) & =\left(\delta_0^{n-2}\left(x_1\right)(t),\ldots,\delta_0^{n-2}\left(x_n\right)(t)\right)
\end{align*}
are linearly independent in $\mathbb{C}^n$ and therefore, there exist $f_1(t),\ldots,f_{n-2}(t)\in \mathbb{C}$ such that
$$[\mathbf{x}+f_1\delta_0\left(\mathbf{x}\right)+\ldots+f_{n-2}\delta_0^{n-2}\left(\mathbf{x}\right)](t)\in C\subseteq\mathbb{P}^{n-1}(\mathbb{C}).$$

As in Section \ref{orbitspace}, let $\Lambda\in\mathbb{Z}_{>0}$ be such that the homogeneous elements of $\mathbb{C}[X_1,\ldots,X_n]^G$ of degree $\Lambda$ form a homogeneous coordinate system for $\mathbb{P}^{n-1}(\mathbb{C})/G$ and let $\Phi_1,\ldots,\Phi_M$ be a basis of $\mathbb{C}[X_1,\ldots,X_n]^G_\Lambda$. Let $\Psi_1,\ldots,\Psi_m\in\mathbb{C}[\Phi_1,\ldots,\Phi_M]\subseteq\mathbb{C}[X_1,\ldots,X_n]$ be homogeneous polynomials such that
$$C=\Big\{[X_1,\ldots,X_n]\in\mathbb{P}^{n-1}(\mathbb{C})\Big|\ \Psi_i([X_1,\ldots,X_n])=0,\ i=1,\ldots,m\Big\}.$$
Therefore, $f_1,\ldots,f_{n-2}$ satisfy the system of equations
$$\Psi_i\left(\mathbf{x}+f_1\delta_0\left(\mathbf{x}\right)+\ldots+f_{n-2}\delta_0^{n-2}\left(\mathbf{x}\right)\right)=0,\ i=1,\ldots,m.$$
Since $C$ is a curve of degree $d$, the system can be solved in a field extension $E$ of degree $d$ over the field of definition of the system. Now, we will prove that field of definition of the system is $F_0$. Let $d_i$ be the degree of $\Psi_i$. By Proposition \ref{propgau}, if $Q_1,\ldots, Q_{m_i}$ is a basis of $\mathbb{C}[Y_1,\ldots,Y_n]^G_{d_i}$, then
$$\Psi_i\left(\mathbf{x}+f_1\delta_0\left(\mathbf{x}\right)+\ldots+f_{n-2}\delta_0^{n-2}\left(\mathbf{x}\right)\right)=g_1Q_1(\mathbf{x})+\ldots+g_{m_i}Q_{m-i}(\mathbf{x}),$$
where $g_i\in \mathbb{C}[f_0,\ldots,f_{n-1}]^G_{d_i}$ with $f_0=1$ and $f_{n-1}=0$. By Galois correspondence, $Q_1(\mathbf{x}),\ldots,Q_{m-i}(\mathbf{x})\in K_1$ and therefore, the field of definition of the system is $K_1$.

Let $L'(x)=0$ be the LODE with fundamental system of solution $\mathbf{x}+f_1\delta_0\left(\mathbf{x}\right)+\ldots+f_{n-2}\delta_0^{n-2}\left(\mathbf{x}\right)$. Since the image of the Schwarz map of this fundamental system of solutions is $C$, the theorem follows from applying Theorem \ref{theopullback} to $L'(x)=0$.
\end{proof}

\begin{cor}\label{coro}
Under the hypothesis and notation of Theorem \ref{thetheo}, there exist $g_0,\ldots,g_n\in E$ such that
$$x_i=g_0f^{1/\Lambda}\cdot y_i\circ h+g_1\delta(f^{1/\Lambda}\cdot y_i\circ h)\ldots+g_n\delta^n(f^{1/\Lambda}\cdot y_i\circ h).$$
In particular, $L(x)=0$ is gauge equivalent over $E$ to an LODE projectively equivalent to a pullback by a map of $L_C(y)=0$ in $E_1$.
\end{cor}

\begin{proof}
Let
\[
L(x)=\delta^n(x)+a_{n-1}\delta^{n-1}(x)+\ldots+a_1\delta(x)+a_0x.
\]
By Theorem \ref{thetheo}, we have
$$x_i+f_1\delta_0(x_i)+\ldots+f_{n-2}\delta_0^{n-2}(x_i)=f^{1/\Lambda}\cdot y_i\circ h.$$
Differentiating with $\delta$ this equation $(n-1)$-times and using the relation
$$\delta^n(x)=-a_{n-1}\delta^{n-1}(x)-\ldots-a_1\delta(x)-a_0x,$$
we obtain
{\small $$
\left[\begin{array}{ccccc}
1 & f_{1} & \hdots &  f_{n-2} & 0\\
0 & 1+f'_{1} & \hdots & f_{n-3}+f'_{n-2} & f_{n-2}\\
\vdots & \vdots & \vdots & \vdots & \vdots\\
* & * & \hdots &  * & *
\end{array}\right]
\left[\begin{array}{c}
x_i\\
\delta(x_i)\\
\vdots\\
\delta^{(n-1)}(x_i)
\end{array}\right]=
\left[\begin{array}{c}
f^{1/\Lambda}\cdot y_i\circ h \\
\delta\big(f^{1/\Lambda}\cdot y_i\circ h\big) \\
\vdots\\
\delta^{n-1}\big(f^{1/\Lambda}\cdot y_i\circ h\big)
\end{array}\right].
$$}
Since $x_i$, $i=1,\ldots,n$ and $f^{1/\Lambda}\cdot y_i\circ h$, $i=1,\ldots,n$ form systems of solutions to irreducible LODEs, then the vectors
$$\left\{\left[\begin{array}{c}
x_i\\
\delta(x_i)\\
\vdots\\
\delta^{(n-1)}(x_i)
\end{array}\right],\quad i=1,\ldots,n\right\}$$
and
$$\left\{\left[\begin{array}{c}
f^{1/\Lambda}\cdot y_i\circ h \\
\delta\big(f^{1/\Lambda}\cdot y_i\circ h\big) \\
\vdots\\
\delta^{n-1}\big(f^{1/\Lambda}\cdot y_i\circ h\big)
\end{array}\right],\quad i=1,\ldots,n\right\}$$
are linearly independent and therefore, the matrix
$$\left[\begin{array}{ccccc}
1 & f_{1} & \hdots &  f_{n-2} & 0\\
0 & 1+f'_{1} & \hdots & f_{n-3}+f'_{n-2} & f_{n-2}\\
\vdots & \vdots & \vdots & \vdots & \vdots\\
* & * & \hdots &  * & *
\end{array}\right]$$
is invertible. The corollary follows from the first entry of the vector equality
{\small $$
\left[\begin{array}{c}
x_i\\
\delta(x_i)\\
\vdots\\
\delta^{(n-1)}(x_i)
\end{array}\right]=
\left[\begin{array}{ccccc}
1 & f_{1} & \hdots &  f_{n-2} & 0\\
0 & 1+f'_{1} & \hdots & f_{n-3}+f'_{n-2} & f_{n-2}\\
\vdots & \vdots & \vdots & \vdots & \vdots\\
* & * & \hdots &  * & *
\end{array}\right]^{-1}
\left[\begin{array}{c}
f^{1/\Lambda}\cdot y_i\circ h \\
\delta\big(f^{1/\Lambda}\cdot y_i\circ h\big) \\
\vdots\\
\delta^{n-1}\big(f^{1/\Lambda}\cdot y_i\circ h\big)
\end{array}\right].
$$}
\end{proof}

\begin{rem}
Note that in the case $n=2$ we have
$$ \left[\begin{array}{cc}
1 & 0\\
* & *
\end{array}\right]\left[\begin{array}{c}
x_i\\ \delta(x_i)
\end{array}\right]=\left[\begin{array}{c}
f^{1/\Lambda}\cdot y_i\circ h \\
\delta\big(f^{1/\Lambda}\cdot y_i\circ h\big) \\
\end{array}\right].
$$
By multiplying both sides of the equation by the inverse of the matrix, we recover the closed form from Klein's theorem
$$x_i=f^{1/\Lambda}\cdot y_i\circ h.$$
\end{rem}

\section{Third order algebraic LODEs}

In this section we will present generalized hypergeometric equations that can be used as standard equations for third order algebraic LODEs. We will also give the formula to obtain the pullback function in Theorem \ref{thetheo} for these equations. Moreover, we will illustrate Corollary \ref{coro}  with an example.

There are, up to isomorphism, $8$ primitive subgroups of $SL_3(\mathbb{C})$. As in \cite{SINGER1993}, we define the following matrices

$$
E_1=\left[\begin{array}{ccc}
1 & 0 & 0\\
0 & \xi^4 & 0\\
0 & 0 & \xi
\end{array}\right],\quad
E_2=\left[\begin{array}{ccc}
-1 & 0 & 0\\
0 & 0 & -1\\
0 & -1 & 0
\end{array}\right],\quad
E_3=\dfrac{1}{\sqrt{5}}\left[\begin{array}{ccc}
1 & 2 & 2\\
1 & s & t\\
1 & t & s
\end{array}\right],
$$
$$
E_4=\dfrac{1}{\sqrt{5}}\left[\begin{array}{ccc}
1 & 2\lambda_2 & 2\lambda_2\\
\lambda_1 & s & t\\
\lambda_1 & t & s
\end{array}\right], \quad
S=\left[\begin{array}{ccc}
\beta & 0 & 0\\
0 & \beta^2 & 0\\
0 & 0 & \beta^4
\end{array}\right], \quad
R=\dfrac{1}{\sqrt{7}i}\left[\begin{array}{ccc}
a  & b & c\\
b & c & a\\
c & a  & b
\end{array}\right],
$$
$$
S_1=\left[\begin{array}{ccc}
1 & 0 & 0\\
0 & \omega & 0\\
0 & 0 & \omega^2
\end{array}\right], \quad
T=\left[\begin{array}{ccc}
0 & 1 & 0\\
0 & 0 & 1\\
1 & 0  & 0
\end{array}\right], \quad
U=\left[\begin{array}{ccc}
\varepsilon & 0 & 0\\
0 & \varepsilon & 0\\
0 & 0  & \varepsilon\omega
\end{array}\right],
$$
$$
V=\rho\left[\begin{array}{ccc}
1 & 1 & 1\\
1 & \omega & \omega^2\\
1 & \omega^2  & \omega
\end{array}\right],\text{ and }
Z=\left[\begin{array}{ccc}
\omega & 0 & 0\\
0 & \omega & 0\\
0 & 0  & \omega
\end{array}\right].
$$
Here, $\xi$ is a primitive $5$th root of unity, $s=\xi^3+\xi^2$, and $t=\xi^4+\xi$. Note that $t-s=\sqrt{5}$. Furthermore, $\varepsilon$ is a primitive $9$th root of unity. Thus, $\varepsilon^6+\varepsilon^3+1=0$ and $\omega=-1-\varepsilon^3$ is a primitive $3$rd root of unity. We have that $\beta$ is a primitive $7$-th root of unity, $a=\beta^4-\beta^3$, $b=\beta^2-\beta^5$, $c=\beta-\beta^6$, $\lambda_1=\frac{-1+\sqrt{15}i}{4}$, $\lambda_2=\frac{-1-\sqrt{15}i}{4}$, and $\rho=\frac{1}{\omega-\omega^2}$. Note that $\beta^6+\beta^5+\beta^3-\beta^4-\beta^2-\beta=\sqrt{7}i$. Then, the primitive subgroups of $SL_3(\mathbb{C})$ are the following.

\begin{itemize}
\item The Klein group $G_{168}=\langle R,S,T\rangle$ and its direct product with the cyclic group of order three $G_{168}\times C_{3}=\langle R,S,T,Z\rangle$. Their projection into $PGL_3(\mathbb{C})$ are both isomorphic to $G_{168}$.
\item The group $H_{216}^{SL_3}=\langle S_1,T,V,U\rangle$.  Its projection into $PGL_3(\mathbb{C})$ is isomorphic to the Hessian group $H_{216}$.
\item The group $H_{72}^{SL_3}=\langle S_1,T,V,UVU^{-1}\rangle$. Its projection into $PGL_3(\mathbb{C})$ is the normal subgroup of the Hessian group $H_{72}$.
\item The group $F_{36}^{SL_3}=\langle S_1,T,V\rangle$. Its projection into $PGL_3(\mathbb{C})$ is the normal subgroup of the Hessian group $F_{36}$.
\item The Valentiner group $A_6^{SL_3}=\langle E_1,E_2,E_3,E_4\rangle$. Its projection into $PGL_3(\mathbb{C})$ is isomorphic to $A_6$.
\item The alternating group $A_5=\langle E_1,E_2,E_3\rangle$ and its direct product with the cyclic group of order three $A_5\times C_3=\langle E_1,E_2,E_3,Z\rangle$. Their projection into $PGL_3(\mathbb{C})$ are both isomorphic to $A_5$.
\end{itemize}

In order to simplify the computations involved in obtaining the coefficients $f_i$'s from Theorem \ref{thetheo}, for each maximal finite primitive group we will choose a hypergeometric standard equation that is projectively equivalent to one of the equations from Beuker and Heckman's list \cite{BEUKERS1989, KATO2006}. The equation for $F_{36}^{SL_3}$ is based on Geiselmann and Ulmer's equation \cite{GEISELMAN1997}.

\subsection{The Klein group $G_{168}$ and $G_{168}\times C_3$}

Since the groups $G_{168}$ and $G_{168}\times C_3$ have the same projection into $PGL_3(\mathbb{C})$, projective curves under one group are invariant under the other. Therefore, it suffices to produce standard equations for $G_{168}$.

The invariant subring $\mathbb{C}[X_1,X_2,X_3]^{G_{168}}$ is generated by
\begin{eqnarray*}
F_4 & = & X_1^3X_2+X_2^3X_3+X_3^3X_1,\\
F_6 & = & \frac{1}{54}\det\left[\begin{array}{ccc}
\partial^2 F_4/\partial X_1\partial X_1 & \partial^2 F_4/\partial X_1\partial X_2 & \partial^2 F_4/\partial X_1\partial X_3 \\
\partial^2 F_4/\partial X_2\partial X_1 & \partial^2 F_4/\partial X_2\partial X_2 & \partial^2 F_4/\partial X_2\partial X_3 \\
\partial^2 F_4/\partial X_3\partial X_1 & \partial^2 F_4/\partial X_3\partial X_2 & \partial^2 F_4/\partial X_3\partial X_3 \\
\end{array}\right],\\
F_{14} & =  & \frac{1}{9}\det\left[\begin{array}{cccc}
\partial^2 F_4/\partial X_1\partial X_1 & \partial^2 F_4/\partial X_1\partial X_2 & \partial^2 F_4/\partial X_1\partial X_3 & \partial F_6/\partial X_1\\
\partial^2 F_4/\partial X_2\partial X_1 & \partial^2 F_4/\partial X_2\partial X_2 & \partial^2 F_4/\partial X_2\partial X_3 & \partial F_6/\partial X_2\\
\partial^2 F_4/\partial X_3\partial X_1 & \partial^2 F_4/\partial X_3\partial X_2 & \partial^2 F_4/\partial X_3\partial X_3 & \partial F_6/\partial X_3\\
\partial F_6/\partial X_1 & \partial F_6/\partial X_2 & \partial F_6/\partial X_3 & 0
\end{array}\right],
\end{eqnarray*}
and
\begin{eqnarray*}
F_{21} & = & \frac{1}{14}\det\left[\begin{array}{ccc}
\partial F_4/\partial X_1 & \partial F_4/\partial X_2 & \partial F_4/\partial X_3 \\
\partial F_6/\partial X_1 & \partial F_6/\partial X_2 & \partial F_6/\partial X_3 \\
\partial F_{14}/\partial X_1 & \partial F_{14}/\partial X_2 & \partial F_{14}/\partial X_3
\end{array}\right].
\end{eqnarray*}
As a ring, $\mathbb{C}[X_1,X_2,X_3]^{G_{168}}$ is isomorphic to $\mathbb{C}[Z_4,Z_6,Z_{14},Z_{21}]/(T)$, where
\begin{eqnarray*}
T & = & Z_{21}^2+2048Z_4^9Z_6 -22016Z_4^6Z_6^3 +256Z_{14}Z_4^7 +60032Z_4^3Z_6^5 -1088Z_{14}Z_4^4Z_6^2\\ & & -1728Z_6^7 -1008Z_{14}Z_4Z_6^4 +88Z_{14}^2Z_4^2Z_6 -Z_{14}^3.
\end{eqnarray*}

The equation with solution ${}_3F_2(-1/42,5/42,17/42;1/3,2/3| t)$ is
\begin{eqnarray*}
0 & = & \left(\frac{d}{dt}\right)^3y+\frac{1}{2}\frac{7t-4}{t(t-1)}\left(\frac{d}{dt}\right)^2y+\frac{1}{252}\frac{387t-56}{t^2(t-1)}\left(\frac{d}{dt}\right)y \\
 & & \ -\frac{1}{74088}\frac{85}{t^2(t-1)}y .
\end{eqnarray*}
The image of the Schwarz map describes Klein's quartic, $F_4=0$. If $\mathbf{y}=(y_1,y_2,y_3)$ is a full system of solutions, we have
\begin{eqnarray*}
F_6(\mathbf{y}) & = & 1\\
\dfrac{F_{14}^3}{1728F_6^7}(\mathbf{y}) & = & t
\end{eqnarray*}
and therefore, the equation is standard.

\subsubsection{Example}

Here, we will ilustrate Corollary \ref{coro} using an equation with differential Galois group $G_{168}$. Let $\mathbf{x}=(x_1,x_2,x_3)$ be a full system of solutions to van der Put-Ulmer's equation \cite[Section 8.2.1, \emph{Branch type 2,7,7}]{VANDERPUT2000}
\begin{eqnarray*}
0 & = & \left(\frac{d}{dt}\right)^3x+\frac{7t-2}{t(t-1)}\left(\frac{d}{dt}\right)^2x +\frac{1}{28}\frac{288t^2-161t-7}{t^2(t-1)^2}\left(\frac{d}{dt}\right)x \\
 & & \ +\frac{1}{2744}\frac{6336t^3-5273t^2+343t-686}{t^3(t-1)^3}x.
\end{eqnarray*}
The image of the Schwarz map describes the curve $-\dfrac{7}{53}F_6^3-\dfrac{1}{8}F_4F_{14}+F_{4}^3F_6=0$ \cite{VANDERPUT2020}. For a generic $f_1$, we have
$$ F_4(\mathbf{x}+f_1\mathbf{x}') = \frac {p_1(t)p_2(t)p_3(t)}{38416 \left( t-1
 \right) ^{7}{t}^{6}},$$
where $p_1(t) = 14\,{t}^{2}-14\,t-(19\,t-7)\,f_{{1}}$,
$p_2(t)=14\,{t}^{2}-14\,t-(16\,t-7)\,f_{{1}}$, and
$p_3(t)=1372\,{t}^{4}-2744\,{t}^{3}+1372\,{t}^{2}-(3430\,{t}^{3}f_{{1}}-4802\,{t}^{
2}+1372\,t)f_{{1}}+(2146\,{t}^{2}-1715\,t+343)\,{f_{{1}}}^{2}$.
Let $f_1(t)=(14\,{t}^{2}-14\,t)/(19\,x-7)$. Then, $p_1(t)=0$ and the map $[\mathbf{x}+f_1\mathbf{x}']$ describes Klein's quartic. Let
$$ f(t) = F_6(\mathbf{x}+f_1\mathbf{x}')=\dfrac{3^6t^3}{(t-1)^4(19t-7)^6}$$
and
$$h(t) = \dfrac{F_{14}^3}{1728F_6^7}(\mathbf{x}+f_1\mathbf{x}')=\dfrac{(t+3)^3}{27(t-1)^2}.$$
From Theorem \ref{thetheo}, we obtain $\mathbf{x}(t)+f_1(t)\mathbf{x}'(t)=f(t)^{1/6}\mathbf{y}\left(h(t)\right)$. If we denote
$$F_{G_{168}}(t)={}_3F_2(-1/42,5/42,17/42;1/3,2/3| t),$$
then, by Corollary \ref{coro}, we obtain a solution
\begin{align*}
x(t) & =\dfrac{1932781}{6049137024}\dfrac{1}{t^{1/2}(t-1)^{17/3}}\Bigg\{ t(t-9)^2(t+3)^4F''_{G_{168}}\left(\dfrac{(t+3)^3}{27(t-1)^2}\right)\\
 & \quad +\dfrac{5821200}{113693}(t-1)^2\bigg\{(t+3)(t^3-\dfrac{351}{55}t^2+\dfrac{243}{5}t+\dfrac{567}{55})F'_{G_{168}}\left(\dfrac{(t+3)^3}{27(t-1)^2}\right)\\
 & \qquad +\dfrac{433944}{4675} (t-\dfrac{21}{41})(t-1)^2F_{G_{168}}\left(\dfrac{(t+3)^3}{27(t-1)^2}\right)\bigg\}\Bigg\}.
\end{align*}

\subsection{The Hessian group $H_{216}$ and its subgroups $H_{72}$ and $F_{36}$}

We define $P =X_1X_2X_3$, $S = X_1^3+X_2^3+X_3^3$, $Q=X_1^3X_2^3+X_1^3X_3^3+X_2^3X_3^3$, $R = (X_1^3-X_2^3)(X_1^3-X_3^3)(X_2^3-X_3^3)$,
$F_6 = S^2-12Q$, $\Phi_6 =S^2-18P^2-6PS$, $F_{12} =S^4+216P^3S$, $\Phi_{12} =P(27P^3-S^3)$, and $\Psi_{12}=PS^3+3P^2S^2-18P^3S.$

Then, the ring of invariants for $F_{36}^{SL_3}$, $H_{72}^{SL_3}$,  $H_{216}^{SL_3}$ are
$$\mathbb{C}[X_1,X_2,X_3]^{F_{36}^{SL_3}}  =\mathbb{C}[F_6,\Phi_6,R,F_{12},\Psi_{12}],$$
$$\mathbb{C}[X_1,X_2,X_3]^{H_{72}^{SL_3}}  =\mathbb{C}[F_6,R,F_{12},\Phi_6^2],$$
and
$$\mathbb{C}[X_1,X_2,X_3]^{H_{216}^{SL_3}}  =\mathbb{C}[R,\Phi_{12},F_6F_{12},F_6^3],$$
respectively. As rings, they are isomorphic to
$$\mathbb{C}[Z_6,Y_6,Z_9,Z_{12},Y_{12}]/(T_{18},T_{24}),$$
$$\mathbb{C}[Z_6,Z_9,Z_{12},X_{12}]/(T_{36}),$$
and
$$\mathbb{C}[Z_9,Y_{12},Z_{18},Y_{18}]/(T_{54}),$$
respectively, where
$$ T_{18} =432Z_9^2-Z_6^3+3Z_6Z_{12}-2Y_6^3-36Y_6Y_{12},$$
$$T_{24}  = Z_{12}^2-Y_6(Z_{12}+12Y_{12})+12Y_{12}^2,$$
$$T_{36}  =(432Z_9^2+3Z_6Z_{12}-Z_6^3)^2-4(X_{12}^3-3Z_{12}X_{12}^2+3Z_{12}^2X_{12}),$$
and
$$T_{54}  =Z_{18}^3-\dfrac{1}{4}\big((432Z_9^2-Y_{18}+3Z_{18})^2-4\cdot 1728Y_{12}^3\big)Y_{18}.$$

The equation with solution ${}_3F_2(17/36, 2/9, -1/36;1/3, 2/3| \frac{1}{t})$ is
\begin{eqnarray*}
0 & = & \left(\frac{d}{dt}\right)^3y+\frac{1}{3}\frac{11t-6}{t(t-1)}\left(\frac{d}{dt}\right)^2y+\frac{1}{432}\frac{757t-96}{t^2(t-1)}\left(\frac{d}{dt}\right)y \\
 & & \ -\frac{1}{5832}\frac{17}{t^3(t-1)}y .
\end{eqnarray*}
The image of the Schwarz map describes the curve $F_6=0$. If $\mathbf{y}=(y_1,y_2,y_3)$ is a full system of solutions, then
\begin{eqnarray*}
R(\mathbf{y}) & = & 1\\
\dfrac{6^6R^4}{\Phi_{12}^3}(\mathbf{y}) & = & t
\end{eqnarray*}
and therefore, the equation is standard.

Note that since
$$\Phi_{12}=\dfrac{\Phi_6^2-F_{12}}{12},$$
we can obtain the value of $\Phi_{12}(\mathbf{y})$ from the invariants under the actions of $H_{72}^{SL_3}$ or $F_{36}^{SL_3}$.

\begin{rem} There is no hypergeometric equation with Galois group $H_{72}^{SL_3}$. On the other hand, there exist hypergeometric equations with Galois group $F_{36}^{SL_3}$. For example, the equation with solution ${}_3F_2(-1/12,1/6,5/12;1/4,3/4| \frac{1}{t})$
\begin{eqnarray*}
0 & = & \left(\frac{d}{dt}\right)^3y+\frac{1}{2}\frac{8t-5}{t(t-1)}\left(\frac{d}{dt}\right)^2y+\frac{5}{48}\frac{21t-5}{t^2(t-1)}\left(\frac{d}{dt}\right)y \\
 & & \ -\frac{1}{864}\frac{5}{t^3(t-1)}y .
\end{eqnarray*}
Let $F_3=6\sqrt{3}P+(\sqrt{3}+3)S$ and $\Phi_3=6\sqrt{3}P+(\sqrt{3}-3)S$. The polynomials $F_3$ and $\Phi_3$ are semi-invariants of degree $3$ such that $F_3\Phi_3=\Phi_6$. The image of the Schwarz map describes the elliptic curve $F_3=0$. If $\mathbf{y}=(y_1,y_2,y_3)$ is a full system of solutions, then
\begin{eqnarray*}
F_6(\mathbf{y}) & = & 1\\
\dfrac{F_{6}^3}{F_{6}^3-432R^2}(\mathbf{y}) & = & t
\end{eqnarray*}
and therefore, the equation is standard.
\end{rem}

\subsection{The alternating group $A_6$}

The invariant subring $\mathbb{C}[X_1,X_2,X_3]^{A_{6}^{SL_3}}$ is generated by
\begin{eqnarray*}
F_6 & = & 10X_1^3X_2^3+9X_1^5X_3+9X_2^5X_3-45X_1^2X_2^2X_3^2-135X_1X_2X_3^4+27X_3^6,\\
F_{12} & = & -\frac{1}{20250}\det\left[\begin{array}{ccc}
\partial^2 F_6/\partial X_1\partial X_1 & \partial^2 F_6/\partial X_1\partial X_2 & \partial^2 F_6/\partial X_1\partial X_3 \\
\partial^2 F_6/\partial X_2\partial X_1 & \partial^2 F_6/\partial X_2\partial X_2 & \partial^2 F_6/\partial X_2\partial X_3 \\
\partial^2 F_6/\partial X_3\partial X_1 & \partial^2 F_6/\partial X_3\partial X_2 & \partial^2 F_6/\partial X_3\partial X_3 \\
\end{array}\right],\\
F_{30} & =  & \frac{1}{24300}\det\left[\begin{array}{cccc}
\partial^2 F_6/\partial X_1\partial X_1 & \partial^2 F_6/\partial X_1\partial X_2 & \partial^2 F_6/\partial X_1\partial X_3 & \partial F_{12}/\partial X_1\\
\partial^2 F_6/\partial X_2\partial X_1 & \partial^2 F_6/\partial X_2\partial X_2 & \partial^2 F_6/\partial X_2\partial X_3 & \partial F_{12}/\partial X_2\\
\partial^2 F_6/\partial X_3\partial X_1 & \partial^2 F_6/\partial X_3\partial X_2 & \partial^2 F_6/\partial X_3\partial X_3 & \partial F_{12}/\partial X_3\\
\partial F_{12}/\partial X_1 & \partial F_{12}/\partial X_2 & \partial F_{12}/\partial X_3 & 0
\end{array}\right],
\end{eqnarray*}
and
\begin{eqnarray*}
F_{45} & = & \frac{1}{4860}\det\left[\begin{array}{ccc}
\partial F_6/\partial X_1 & \partial F_6/\partial X_2 & \partial F_6/\partial X_3 \\
\partial F_12/\partial X_1 & \partial F_12/\partial X_2 & \partial F_12/\partial X_3 \\
\partial F_{30}/\partial X_1 & \partial F_{30}/\partial X_2 & \partial F_{30}/\partial X_3
\end{array}\right].
\end{eqnarray*}
As a ring, $\mathbb{C}[X_1,X_2,X_3]^{A_{6}^{SL_3}}$ is isomorphic to $\mathbb{C}[Z_6,Z_{12},Z_{30},Z_{45}]/(T)$, where
\begin{eqnarray*}
T & = & 4Z_{6}^13Z_{12}+80Z_{6}^11Z_{12}^2+816Z_{6}^9Z_{12}^3+18Z_{6}^10Z_{30}+4376Z_{6}^7Z_{12}^4+198Z_{6}^8Z_{12}Z_{30}\\
& & \quad +13084Z_{6}^5Z_{12}^5+954Z_{6}^6Z_{12}^2Z_{30}+12312Z_{6}^3Z_{12}^6-198Z_{6}^4Z_{12}^3Z_{30}+5616Z_{6}Z_{12}^7\\
& & \quad\ -162Z_{6}^5Z_{30}^2-5508Z_{6}^2Z_{12}^4Z_{30}-1944Z_{6}^3Z_{12}Z_{30}^2-1944Z_{12}^5Z_{30}-1458Z_{6}Z_{12}^2Z_{30}^2\\
& & \ \quad\ +729Z_{30}^3-19683Z_{45}^2.\\
\end{eqnarray*}

The equation with solution ${}_3F_2(-1/60,11/60,7/12;1/2,3/4| t)$ is
\begin{eqnarray*}
0 & = & \left(\frac{d}{dt}\right)^3y+\frac{3}{4}\frac{5t-3}{t(t-1)}\left(\frac{d}{dt}\right)^2y+\frac{1}{1200}\frac{2213t-450}{t^2(t-1)}\left(\frac{d}{dt}\right)y \\
 & & \ -\frac{1}{43200}\frac{77}{t^2(t-1)}y .
\end{eqnarray*}
The image of the Schwarz map the curve $F_6=0$. If $\mathbf{y}=(y_1,y_2,y_3)$ is a full system of solutions, we have
\begin{eqnarray*}
F_{12}(\mathbf{y}) & = & 1\\
\dfrac{3F_{30}^2}{8F_{12}^5}(\mathbf{y}) & = & t
\end{eqnarray*}
and therefore, the equation is standard.

\subsection{The alternating group $A_5$ and the group $A_5\times C_3$}

The invariant subring $\mathbb{C}[X_1,X_2,X_3]^{A_{6}^{SL_3}}$ is generated by
\begin{eqnarray*}
F_2 & = & X_1^2+X_2X_3,\\
F_{6} & = & 8X_1^4X_2X_3-2X_1^2X_2^2X_3^2-X_1(X_2^5+X_3^5)+X_2^3X_3^3,\\
F_{10} & =  & 320 X_1^6 X_2^2 X_3^2-160 X_1^4 X_2^3 X_3^3+20 X_1^2 X_2^4 X_3^4+6 X_2^5 X_3^5\\
 & & \quad-4 X_1 (X_2^5+X_3^5) (32 X_1^4-20 X_1^2 X_2 X_3+5 X_2^2 X_3^2)+X_2^{10}+X_3^{10},
\end{eqnarray*}
and
\begin{eqnarray*}
F_{15} & = & (X_2^5-X_3^5) (-1024 X_1^{10}+3840 X_1^8 X_2 X_3-3840 X_1^6 X_2^2 X_3^2+1200 X_1^4 X_2^3 X_3^3\\
& & \quad -100 X_1^2 X_2^4 X_3^4+X_2^{10}+X_3^{10}+2 X_2^5 X_3^5\\
& & \qquad +X_1 (X_2^5+X_3^5) (352 X_1^4-160 X_1^2 X_2 X_3+10 X_2^2 X_3^2)).
\end{eqnarray*}
As a ring, $\mathbb{C}[X_1,X_2,X_3]^{A_{5}}$ is isomorphic to $\mathbb{C}[Z_2,Z_{6},Z_{10},Z_{15}]/(T)$, where
\begin{eqnarray*}
T & = &  Z_{15}^2+1728 Z_{6}^5- Z_{10}^3-720 Z_{2} Z_{6}^3 Z_{10}+80 Z_{2}^2 Z_{6} Z_{10}^2-64 Z_{2}^3(5 Z_{6}^2- Z_{2} Z_{10})^2.\\
\end{eqnarray*}

The equation with solution ${}_3F_2(-1/30,1/6,11/30;1/3,2/3| t)$ is
\begin{eqnarray*}
0 & = & \left(\frac{d}{dt}\right)^3y+\frac{1}{2}\frac{7t-4}{t(t-1)}\left(\frac{d}{dt}\right)^2y+\frac{1}{900}\frac{1389t-200}{t^2(t-1)}\left(\frac{d}{dt}\right)y \\
 & & \ -\frac{1}{5400}\frac{11}{t^2(t-1)}y .
\end{eqnarray*}
It is the symmetric square of the equation with solution ${}_2F_1(-1/60,11/60;2/3| t)$, a standard equation for $A_5^{SL_2}$. In particular the image of the Schwarz map is the rational curve $F_2=0$. If $\mathbf{y}=(y_1,y_2,y_3)$ is a full system of solutions, we have
\begin{eqnarray*}
F_{6}(\mathbf{y}) & = & 1\\
\dfrac{F_{10}^3}{1728F_{6}^5}(\mathbf{y}) & = & t
\end{eqnarray*}
and therefore, the equation is standard.

\bibliographystyle{plain}