function e = genell(a,b,num_pts)
%
a2 = a^2;
b2 = b^2;
x = -a:(1/num_pts):a;
x2 = x.^2;
y_pos = sqrt(b2*(1-(x2/a2)));
y_neg = -y_pos;
e = [x' y_pos'; x' y_neg'];
