We present several algorithms for computing minimum energy reaction paths by the method of steepest descents of the potential energy and for interpolating reaction-path potentials. The algorithms are suited for use with global analytic potential surfaces or ab initio calculations of energies and gradients. The efficiencies of the various algorithms are demonstrated and compared by applications to several test cases.