function P = ResolveRede(nv,conec,C)

  [nc dummy]=size(conec);
   A = zeros(nv);
   for k=1:nc 
     p = conec(k,1);
     q = conec(k,2);
     A(p,p) = A(p,p) + C(k);
     A(p,q) = A(p,q) - C(k);
     A(q,p) = A(q,p) - C(k);
     A(q,q) = A(q,q) + C(k);
   end
   b=-0.1*ones(nv,1);
Atilde=A; Atilde(1,:)=0; Atilde(1,1)=1;
b(1,1)=5;
P=Atilde\b;
end