function y=rootfcn(lambda)
global vv s2 c
z = vv./(lambda+s2);
y = z'*z-c;
