function fx = felet(x)
global m n nt co T cL ome
ii=x(1:m); v=x(m+1:m+n); iA=x(m+n+1:m+n+nt);
fx=zeros(m+n+nt,1); fx(m+1:m+n)=-cL; ## inic. fx
dv=deltav(ii);             ## funcao G(i)
for k=1:m
 n1=co(k,1); n2=co(k,2);   ## co==conectividade
 fx(k)=dv(k)+v(n1)-v(n2);  ## lei de Ohm
 fx(m+n1)=fx(m+n1)+ii(k);   ## eq. de carga
 fx(m+n2)=fx(m+n2)-ii(k);   ## eq. de carga
endfor
for k=1:nt
 fx(m+T(k))=fx(m+T(k))-iA(k); ## eq. de carga
 fx(m+n+k)=-v(T(k));     ## eq. impos. terra
endfor
endfunction