function [c,n, ceq] = unitdisk(u)
c = lambda-0;
n= c-0;
ceq = integral(u)-c;