Euler type \(x^{3}y^{\prime \prime \prime }+x^{2}y^{\prime \prime }+xy^{\prime }+y=f\left ( x\right ) \)

ode internal name "higher_order_ODE_non_constant_coefficients_of_type_Euler"

This uses same algorithm as for second order Euler type ode but for higher order.