Rishu Saxena, Anne Gelb and Hans Mittelmann. (2009). A High Order Method for Determining the Edges in the Gradient of a Function.
            <i>Communications in Computational Physics</i>. <i>5</i> (2-4).
            694-711.
         doi: