function Joptimum=Jopt(u)


[u_opt,FVAL,EXITFLAG,OUTPUT,LAMBDA,grad,hessian] = fmincon(@col_obj,u0,[],[],[],[],lb,ub,@col_constr,options);

end


