function [f,dfdx,dfdy] = circle(x,y)

f = x.^2 + y.^2 - 1;
dfdx = 2.*x;
dfdy = 2.*y;