Is there an easy way to generate a simple triangular mesh, given an implicit surface defined by

S = {(u,v) : f(u,v) = 0}  where f(u,v) = ± sqrt(1-(v²)-((u-v²)²)) ?

Currently, I have successfully written code to generate the top and bottom half but I am having trouble stitching it into one coherent surface mesh. I feel like there should be a simpler way to directly do it.

I'm working with Julia v1.10.10 and Meshes.jl v0.43.5. I will be working on simulating a pde on this surface.