The formula to determine the length of chord intercepted by a parabola $y^2=4ax$ on a line $y=mx+c$, where the variables have their usual meaning is given by the following:
$$L=\left(\frac 4 {m^2}\right)\sqrt{a(1+m^2)(a-mc)}$$
The above formula was given in my book without any proof. I searched the internet to my best, and could not even find this formula, and so its derivation. I tried to derive this formula to find whether it works for all cases. I tried using equation of chord derived from the parametric form of the coordinates, but it gets even more complicated.
So, How to derive this formula?

Substituting $y=mx+c$ in $y^2=4ax$ gives $$m^2x^2+(2mc-4a)x+c^2=0$$ Now, using $$y_1-y_2=m(x_1-x_2),\quad x_1+x_2=\frac{4a-2mc}{m^2},\quad x_1x_2=\frac{c^2}{m^2}$$ where $(x_1,y_1),(x_2,y_2)$ are the intersection points, we get $$\begin{align}L&=\sqrt{(x_1-x_2)^2+(y_1-y_2)^2} \\\\&=\sqrt{(x_1-x_2)^2+m^2(x_1-x_2)^2} \\\\&=\sqrt{(1+m^2)(x_1-x_2)^2} \\\\&=\sqrt{(1+m^2)((x_1+x_2)^2-4x_1x_2)} \\\\&=\sqrt{(1+m^2)\left(\left(\frac{4a-2mc}{m^2}\right)^2-4\cdot \frac{c^2}{m^2}\right)} \\\\&=\sqrt{(1+m^2)\cdot \frac{16a^2-16amc+4m^2c^2-4m^2c^2}{m^4}} \\\\&=\frac{4\sqrt{a(1+m^2)(a-mc)}}{m^2}\end{align}$$