Aquí están los primeros términos de la serie de la solución en $0$, según Arce.
$$
f \left( x \right) ={\it a_0}+{\it a_1}\,x+{\it a_2}\,{x}^{2}+{
\se a_3}\,{x}^{3}\\+{\frac {{\it a_3}\,{\it a_2}\,{\it a_1}\,{\it
a_0}}{2}}\,{x}^{4}\\+ \left( {\frac {{{\it a_0}}^{2}{{\it a_1}}^{2}{{
\se a_2}}^{2}{\it a_3}}{5}}+{\frac {3\,{\it a_0}\,{\it a_1}\,{{
\se a_3}}^{2}}{10}}+{\frac {{\it a_0}\,{{\it a_2}}^{2}{\it a_3}}{5
}}+{\frac {{{\it a_1}}^{2}{\it a_2}\,{\it a_3}}{10}} \right) {x}^{5
}\\+ \left( {\frac {{{\it a_0}}^{3}{{\it a_1}}^{3}{{\it a_2}}^{3}{
\se a_3}}{15}}+{\frac {2\,{{\it a_3}}^{2}{\it a_2}\,{{\it a_1}}^{2
}{{\it a_0}}^{2}}{5}}+{\frac {{{\it a_0}}^{2}{\it a_1}\,{{\it a_2}
}^{3}{\it a_3}}{5}}+{\frac {{\it a_0}\,{{\it a_1}}^{3}{{\it a_2}}^
{2}{\it a_3}}{10}}+{\frac {3\,{{\it a_3}}^{2}{\it a_0}\,{\it a_2}
}{10}}+{\frac {{{\it a_1}}^{2}{{\it a_3}}^{2}}{10}}+{\frac {{{\it
a_2}}^{2}{\it a_1}\,{\it a_3}}{10}} \right) {x}^{6}\\+\dots
$$
Con $4$ constantes arbitrarias, como se esperaba. Pero, de nuevo, como era de esperar, no es una combinación lineal de 4 soluciones independientes.
agregó ¿Cómo se consigue esto?
De arce. Las salidas son largas, así que omito.
>DE:=diff(f(x),x,x,x,x)=diff(f(x),x,x,x)*diff(f(x),x,x)*diff(f(x),x)*f(x);
>Order:=7;SS:=dsolve({DE,f(0)=a_0,D(f)(0)=a_1,
(D@@2)(f)(0)=2*a_2,(D@@3)(f)(0)=6*a_3},f(x),series,x=0);
>latex(%);