Two types of bifurcation diagrams of cytosolic calcium nonlinear oscillatory systems are presented in rectangular areas determined by two slowly varying parameters. Verification of the periodic dynamics in the two-parameter areas requires solving the underlying model a few hundred thousand or a few million times, depending on the assumed resolution of the desired diagrams (color bifurcation figures). One type of diagram shows period-n oscillations, that is, periodic oscillations having n maximum values in one period. The second type of diagram shows frequency distributions in the rectangular areas. Each of those types of diagrams gives different information regarding the analyzed autonomous systems and they complement each other. In some parts of the considered rectangular areas, the analyzed systems may exhibit non-periodic steady-state solutions, i.e., constant (equilibrium points), oscillatory chaotic or unstable solutions. The identification process distinguishes the later types from the former one (periodic). Our bifurcation diagrams complement other possible two-parameter diagrams one may create for the same autonomous systems, for example, the diagrams of Lyapunov exponents, Ls diagrams for mixed-mode oscillations or the 0–1 test for chaos and sample entropy diagrams. Computing our two-parameter bifurcation diagrams in practice and determining the areas of periodicity is based on using an appropriate numerical solver of the underlying mathematical model (system of differential equations) with an adaptive (or constant) step-size of integration, using parallel computations. The case presented in this paper is illustrated by the diagrams for an autonomous dynamical model for cytosolic calcium oscillations, an interesting nonlinear model with three dynamical variables, sixteen parameters and various nonlinear terms of polynomial and rational types. The identified frequency of oscillations may increase or decrease a few hundred times within the assumed range of parameters, which is a rather unusual property. Such a dynamical model of cytosolic calcium oscillations, with mitochondria included, is an important model in which control of the basic functions of cells is achieved through the Ca2+ signal regulation.