Extensive numerical studies of the dynamo equations due to the global convection are presented to simulate the solar cycle and to open the way to study general stellar magnetic cycles. The dynamo equations which represent the longitudinally-averaged magnetohydrodynamical action (mean magnetohydrodynamics) of the global convection under the influence of the rotation in the solar convection zone are considered here as an initial boundary-value problem. The latitudinal and radial structure of the dynamo action consisting of a generation action due to the differential rotation and a regeneration action due to the global convection is parameterized in accordance with the structure of the rotation and of the global convection. This is done especially in such a way as to represent the presence of the two cells of the regneration action in the radial direction in which the action has opposite signs, which is typical of the regeneration action of the global convection. The effects of the dynamics of the global convection (e.g., the effects of the stratification of the physical conditions in the solar convection zone) are presumed to be all included in those parameters used in the model and they are presumed not to alter the results drastically since these effects are only to change the structure of the regeneration action topologically.