function fx = fhidr2(x)
global m n nr co Ir Hr cL
q=x(1:m); h=x(m+1:m+n); vR=x(m+n+1:m+n+nr);
fx=zeros(m+n+nr,1); fx(m+1:m+n)=-cL; ## inic. fx
dh=deltah2(q);              ## funcao G(q)
for i=1:m
 n1=co(i,1); n2=co(i,2);   ## co==conectividade
 fx(i)=dh(i)+h(n1)-h(n2);  ## eq. de momento
 fx(m+n1)=fx(m+n1)+q(i);   ## eq. de massa
 fx(m+n2)=fx(m+n2)-q(i);   ## eq. de massa
endfor
for i=1:nr
 fx(m+Ir(i))=fx(m+Ir(i))-vR(i); ## eq. de massa
 fx(m+n+i)=-h(Ir(i))+Hr(i);     ## eq. impos. Hr
endfor
endfunction