To get to the standard form, you factorize the nominator and denominator polynomials. Then your polynomials will be of the the form \$K_{1}(s - z_1)(s - z_2)\cdots (s - z_n)\$ and \$K_2(s - p_1)(s - p_2)\cdots (s - p_n)\$. Then identify any complex conjugate pairs among the \$z_k\$ and multiply them out. If, for example, \$z_1 = z_2^*\$, then
$$
(s - z_1)(s - z_2) = s^2 - 2 \mathrm{Re}(z_1) s + |z_1|^2 = |z_1|^2\left(1 - \frac{2\mathrm{Re}(z_1)}{|z_1|^2}s + \frac{s^2}{|z_1|^2}\right).
$$
Now identify
$$
\begin{align}
\omega_1^2 &= |z_1|^2\\
1/Q_1 &= -\frac{2 \mathrm{Re}(z_1)}{|z_1|}
\end{align}
$$
and you get the prescribed form of the second order term.
For the remaining roots \$z_k\$, which will be real, extract the factors as
$$
(s - z_k) = -z_k (1 - \frac{s}{z_k}).
$$
and identify \$\omega_k = -z_k\$.
Repeat for the denominator roots \$p_k\$, and gather the constants to the front to get the factor \$K\$.
The roots you got from Wolfram Alpha are, up to the factors of \$i\$ that connect \$s\$ to \$\omega\$, exactly the \$p_k\$. Sometimes they do indeed end up somewhat hairy, but often it's possible to simplify them by identifying common factors (such as paralleled resistors, products RC that always appear together etc).
Finally, if the polynomial has root \$0\$ with multiplicity \$k\$, these will be factors of the form
$$
\left(\frac{s}{\omega_m}\right)^k,
$$
which you can bring to the front. The factors \$\omega_m\$ are now ambigous, as you can in principle include any of them in \$K\$, but often in practice there's some meaningful choice. For example, if you're designing a filter with a certain passband, you take \$K\$ to be the passband gain (and phase), and take the remaining part to be \$\omega_m^k\$.
The roots \$z_k\$ of the nominator are called the zeroes of the transfer function, as those are the complex values of \$s\$ where the transfer function is indeed the value zero. The roots \$p_k\$ of the denominator are the poles, since those are the values of \$s\$ where the transfer function diverges, which indeed looks like pole sticking out of the \$s\$ -plane if you plot it.
Note that factorizing a polynomial (over the complex numbers) requires finding its roots. For a second order polynomial, the quadratic formula gives you the answer immediately. For third and fourth order polynomials there's the cubic and quartic formulas. The cubic formula is already quite long, and the quartic formula is about a full page in small print, so it's often not useful in practice. For orders higher than five, there is no general formula, although special cases can often be solved.
In addition to using the general formulas, the circuit topology often provides considerable simplifications. For example, in the case of two second order sections separated by a buffer, you can analyze the two sections separately using the quadratic, and the standard form of the combined transfer function is directly the product of the standard forms of the individual sections. The same applies for any number of sections separated by buffers, which is one of the main reasons that high order filter are usually designed as series of second order sections.
If, in the end, you cannot find explicitly the roots, or they're too complicated to use, you can still learn about the your circuit by studying the discriminants, which tells you about potential complex conjugate or real roots. In your specific case (assuming you roots are correct, I didn't check), the discriminant is the term inside the square roots,
$$
\Delta = C_f L_2 R_f^2 + C_f L_1 R_f^2 - 4 L_1 L_2.
$$
If this is negative, you have a complex conjugate pair of roots leading to a second order term, and it's positive, you get two real roots. You can divide by \$L_2\$ and \$C_f\$ to get the expression
$$
\tilde{\Delta} \triangleq R_f^2\left(1 + \frac{L_1}{L_2}\right) - 4 \frac{L_1}{C_f},
$$
which has the same sign as the discriminant. From here you see, for example, that if \$C_f\$ is small enough, or \$R_f\$ is small enough, you get a complex conjugate pair.
Best Answer
Once you find \$\dfrac{v_L}{v_1}\$, then just rearranging the terms will give you the standard form. Let me do this example for you.
Step 1: Calculate \$\dfrac{v_L}{v_1}\$
From the circuit, the following three equations can be derived directly. $$v_L = k_3v_3 \frac{R_L}{R_L + R_{s3}} \tag{1}$$ $$v_3 = -g_{m2}v_2(r_{\pi3}||\dfrac{1}{sc_{\pi3}})\tag{2}$$ $$v_2 = -(g_{m2}v_2 - g_{m1}v_1)(r_{\pi2}||\dfrac{1}{sc_{\pi2}})\tag{3.a}$$ $$ \Rightarrow v_2 = \frac{g_{m1}v_1(r_{\pi2}||\dfrac{1}{sc_{\pi2}})}{1+g_{m2}(r_{\pi2}||\dfrac{1}{sc_{\pi2}})}\tag3$$
Step 2: Rearrange the terms
Now from equation (1), (2) and (3), the value of \$v_L\$ can be expressed in terms of \$v_1\$. From this \$\dfrac{v_L}{v_1}\$ can be calculated. The answer will be in the following form:
$$\frac{v_L}{v_1} = K_1 \times \frac{r_{\pi3}}{1+sc_{\pi3}r_{\pi3}}\times \frac{r_{\pi2}}{1+sc_{\pi2}r_{\pi2}+g_{m2}r_{\pi2}}\tag4$$ where $$K_1= -g_{m1}g_{m2}k_3v_3 \frac{R_L}{R_L + R_{s3}}$$
Now its just a matter of rearranging the terms to convert it into standard form. Dividing numerator and denominator of equation (4) by \$(1+g_{m2}r_{\pi2})\$, $$\Rightarrow \frac{v_L}{v_1} = K_1 \times \frac{r_{\pi3}}{1+sc_{\pi3}r_{\pi3}}\times \frac{\frac{r_{\pi2}}{(1+g_{m2}r_{\pi2})}}{1+s\frac{c_{\pi2}r_{\pi2}}{(1+g_{m2}r_{\pi2})}}$$
$$\Rightarrow \frac{v_L}{v_1} = K \times \frac{1}{1+sc_{\pi3}r_{\pi3}}\times \frac{1}{1+s\frac{c_{\pi2}r_{\pi2}}{(1+g_{m2}r_{\pi2})}}\tag5$$ where $$K = K_1 \times r_{\pi3} \times \frac{r_{\pi2}}{(1+g_{m2}r_{\pi2})}$$
Now equation (5) is in standard form and value of poles can be calculated directly.