% Returns y>0 for any x in the unit circle
function y = circle(x)
y = sqrt(1 - x.*x); % Note that this works when x is an array.