# How does one find solutions to a polynomial over a finite field?

I'm trying to find the solutions to the polynomial $y^2=x^3+1$ over $\mathbb{F}_5$. I have constructed the correct polynomial ring, but I don't know what the analogous function to .roots() is for the two variable case. Thanks