Point(1)={1/2,1/2,0}; //Centre

c=-7/16;
b=4;
a=-4;

npts=100; //Pair

// Parabola with Spline
For i In {1:npts+1}
x=3/8+1/4*(i-1)/npts;
Point(i+1)={x,a*(x)*(x)+b*(x)+c,0};
EndFor
Spline(1)={2:npts+2};

// Exact ellipse
Ellipse(2)={2,1,2,(npts/2+2)};
Ellipse(3)={(npts/2+2),1,2,npts+2};

// Parabola using BSpline
Point(1000)={3/8,1/2,0};
Point(1001)={1/2,1/2+dy,0};
Point(1002)={5/8,1/2,0};
BSpline(4)={1000:1002};

Mesh.ElementOrder = 2; // P2