Summary
We discuss recent developments of the toolbox MatCont for numerical bifurcation analysis. We have added two ways to define models: importing from SBML and with a pseudospectral approximation. The latter allows for investigating systems with delays and renewal equations. Next, we highlight a new interface within MatCont to compute Lyapunov exponents, making it easier for users to check for multi-dimensional tori and chaos, which we demonstrate for a nonlinear forced oscillator. Finally, we discuss homoclinic orbits for a laser model. We highlight the workflow of getting bifurcation curves and how to obtain reproducible data for publications.