У Вас резонанс кратности 2 при l=0, поэтому частное решение нужно искать в виде x^2 * ( A*x^3 + B*x^2 + C*x + D )