function [C Ceq]= col_constr(u1)

global d u NT

u = [u1 ; d(1:3)] ; %L, V, F, zF, qF

x0 = 0.5*ones(1,82)' ;
options = optimset('Display','none','TolFun',1e-10) ;
x = fsolve(@cola_lv,x0,options) ;

xB = x(1) ; 
xD = x(NT) ;

C = [xB-0.01 ;
    0.95-xD] ;

Ceq=[];