Use
   ps_cont_gmres           to compute equilibria along a brach
   check_stability_large   to find out the stability of equilibria
   locate_bif_large        to locate a bifurcation point
        Always edit the above codes and change data
        appropriately

To compare the computing times of different techniques
you have the following codes

   some_comparisons
   further_comparisons
   last_comparisons

Remember that oureigs.m is a simplified version of Matlab's
   command eigs programmed by ourselves, where it is
   implemented a simplified version of the
   Implicitly Restarted Arnoldi method by Sorensen & Lehoucq.

ENJOY THE COMPUTATIONS!