Hi Emeric, This appears to work. Just use RootOf instead of solve: > restart: > f:=(z^2*G^3-z*(2+t*z-2*z)*G^2+(1+t*z^2-2*z^2)*G-1+z); > g:=RootOf(f,G); > series(g,z=0,9); Edwin