function oe = cutEdge(oe, n1, n2)

oe(n1,find(oe(n1,:)==n2)) = -1;
oe(n2,find(oe(n2,:)==n1)) = -1;