function y = circle(foo, x)
   y = sqrt(1 - x.*x);
