After the SBML model is read in via MathSBML, the vector fields are differentiated symbolically to obtain expressions for fx, fp, fxp and fxx. Subsequently, they are numerically evaluated when called by MATCONT. Reparametrization of parameters is done automatically to allow continuation in arbitrary directions. Presently, the toolbox is able to handle bifurcations of saddle-node and Hopf type. Augmentation to include bifurcations of limit cycles is currently underway.