x=[-1:0.01:1];
y_top=sqrt(1 - x.^2);
y_bottom=-sqrt(1 - x.^2);
plot(x,y_top)
hold on
plot (x,y_bottom)