function [c ceq]=cons_fun(x)
c=[];ceq=norm(x)-1;