find border vertex

border_vertex = [];
for i = 1:size(edges,1)
    FI = find(sum(ismember(faces, edges(i,:)),2) == 2);  
    if( dot(normalf(FI(1),:),normalf(FI(2),:)) == -1)
        border_vertex = [border_vertex edges(i,:)];
    end
end

对于每一边, 然后对每个面执行ismember