Add an override to Octave's find() for permutation matrices. Because of find()'s count-limiting and direction arguments, this is slightly more complicated than just copying the permutation vector. I suspect this is a common operation for people who don't know about the 'vector' option to lu().