0001 function [res,res_eq] = cba_opt_constraints(u,c_init,network,constraints) 0002 0003 global my_c 0004 0005 network.kinetics.u = u; 0006 [s,v,sdot] = network_steady_state(network,my_c); 0007 0008 res_sign = - v(isfinite(constraints.v_sign)) .* constraints.v_sign(isfinite(constraints.v_sign)); 0009 res_min = constraints.v_min(isfinite(constraints.v_min)) - v(isfinite(constraints.v_min)); 0010 0011 res = [res_sign; res_min]; 0012 0013 res_eq = [];