Consider a Fuchsian type differential equation written as $$\frac{d^2 y}{dz^2} + \frac{p(z)}{(z - z_1)(z - z_2) \cdots (z - z_m)} \frac{dy}{dz} + \frac{q(z)}{(z - z_1)^2 (z - z_2)^2 \cdots (z - z_m)^2} y = 0, \quad m \geq 2,$$ where $z_1, z_2, \ldots, z_m$ are distinct regular singular points, and $p(z)$ and $q(z)$ are polynomials. How does one relate the solutions to this differential equation to the solutions to the hypergeometric differential equation $$z(1 - z) \frac{d^2 y}{dz^2} + [c - (a + b + 1)z] \frac{dy}{dz} - aby = 0$$ around each regular singular points?
The case $m = 2$ essentially reduces to Riemann's differential equation (if $z = \infty$ is a regular singular point) whose solutions can be written in terms of the hypergeometric functions. Can this be done in general, or for at least four or five regular singular points?
Update: From what I have gathered so far, all homogeneous linear differential equations of the second order having four regular singularities in the extended complex plane, can be transformed into Heun's differential equation whose solutions are related to the hypergeometric functions.