N=2000;
[nv nc conec C coord] = RedeHidraBairro(10,10,2,2);

C(1,1) =20; C(2,1) =20; C(3,1) =20;
C(94,1)=20; C(104,1)=20; C(114,1)=20;
C(31,1)=20; C(32,1)=20; C(33,1)=20;

sum = 0;
for i=1:N
	Cn = Cnew(nc,C); P = ResolveRede(nv,conec,Cn); X=min(P);
 if (X<1.15)
  sum = sum + 1;
 end
 ppp(i)=sum/i;
# paux=reshape(P,10,10);
# contourf(paux,15);
#max(P)
#min(P)
# break
end

prob=sum/N