resultant[ resultant[ x^2 - a x + 1, a^2 - a - 1, a ], x^4000 - b x^2000 + 1, x ]

Uh oh! Wolfram|Alpha doesn't run without JavaScript.

Please enable JavaScript. If you don't know how, you can find instructionshere.Once you've done that, refresh this page to start using Wolfram|Alpha.