One solution of $$x^2 y''' - 3x^2 y'' + 6xy' -6y = 0$$
for $x>0$ is $\phi_1 = x$. Find a basis for the solutions for $x>0$.
My attempt : I try to find $u(x)$ such that $u\phi_1$ is another solution. $$L(ux)=x^2(ux)'''-3x^2(ux)''+6x(ux)'-6(ux)$$
at this point let $u' = v$. We have the equation $$x^2v''+3x^2v'-3x^3v'=0$$ $$(x^3v')'-3x^3v' = 0$$
This is now a linear homogenous differential equation in $(x^3v')$, and I could go on to solve $v$ from here, which might give me one or two $u$ functions to be used for solutions of the original equations.
Can this method work and how should one approach higher order homogenous equations such as this? Doing all these product rules seems very excessive, and I'm wondering if there is a "smart way" I'm missing.