This thesis considers dynamical systems that have multiple time scales. The focus lies on systems with two fast variables and one slow variable. The twoparameter bifurcation structure of the FitzHugh-Nagumo (FHN) equation is analyzed in detail. A singular bifurcation diagram is constructed and invariant manifolds of the problem are computed. A boundary-value approach to compute slow manifolds of saddle-type is developed. Interactions of classical invariant manifolds and slow manifolds explain the exponentially small turning of a homoclinic bifurcation curve in parameter space. Mixed-mode oscillations and maximal canards are detected in the FHN equation. An asymptotic formula to find maximal canards is proved which is based on the first Lyapunov coefficient at a singular Hopf bifurcation.