function  res = reachq(T,HyCase,favor_q)
%this function finds a temperature that achieve the required q
q = calcq(HyCase,T);
res = favor_q - q;