This paper describes an implementation of multistep collocation methods, which are applicable to stiff differential problems, singular perturbation problems, and D.A.E.s of index 1 and 2. These methods generalize one-step implicit Runge-Kutta methods as well as multistep one-stage BDF methods. We give numerical comparisons of our code with two representative codes for these methods, RADAU5 and LSODE.