I am working on the following exercise:
Consider the surface given by the equation $x^2+z^2 = y^3(1-y)^3$. Where can a parametrisation $z(x,y)$ be introduced? Calculate $\partial z / \partial_x$ and $\partial z / \partial_y$.
REMARK: I found out that the surface looks like this here:
I reorded the equation to: $$z^2 = y^3(1-y)^3-x^2$$ , which implies
$$z = \pm\sqrt{y^3(1-y)^3-x^2}$$
So I think a parametrisation of the surface by $z(x,y)$ is not possible, every point $(x,y)$ gets mapped to two different $z$ unless $y^3(1-y)^3-x^2 = 0.$ Is there something I am missing in here?

No, you are missing nothing. For each $(x,y,z)$ in your surface, $(x,y,-z)$ belongs to the surface too. And both points have the same first and second coordinates. So, unless $z=0$, there are two points of the surface with the same first and second coordinates and therefore your surface is not of the form$$\left\{\bigl(x,y,z(x,y)\bigr)\,\middle|\,(x,y)\in D\right\}.$$for some subset $D$ of $\mathbb R^2$ and some function $z\colon D\longrightarrow\mathbb R$.