That looks good when I run your code - when I let it rip on a big array of
triangular coordinates, I get a "Domain error in atan2()"

I'll study your code and try to track down the cause.  Probably
degenerate/coincident points or something.

Post a reply to this message