  We design and test a cone finding algorithm to robustly address nonlinear
system analysis through differential positivity. The approach provides a
numerical tool to study multi-stable systems, beyond Lyapunov analysis. The
theory is illustrated on two examples: a consensus problem with some repulsive
interactions and second order agent dynamics, and a controlled duffing
oscillator.
