[V,D] = eig(C, B);
% vv= (eig(C)); 
opt_val = -D(end,end)