The emergence of multistability in a simple three-dimensional autonomous oscillator is investigated using numerical simulations, calculations of Lyapunov exponents and bifurcation analysis over a broad area of two-dimensional plane of control parameters. Using Neimark–Sacker bifurcation of 1:1 limit cycle as the starting regime, many parameter islands with the coexisting attractors were detected in the phase diagram, including the coexistence of torus, resonant limit cycles and chaos; and transitions between the regimes were considered in detail. The overlapping between resonant limit cycles of different winding numbers, torus and chaos forms the multistability.