The present study is related to design a stochastic framework for the numerical treatment of the Van der Pol heartbeat model (VP-HBM) using the feedforward artificial neural networks (ANNs) under the optimization of particle swarm optimization (PSO) hybridized with the active-set algorithm (ASA), i.e., ANNs-PSO-ASA. The global search PSO scheme and local refinement of ASA are used as an optimization procedure in this study. An error-based merit function is defined using the differential VP-HBM form as well as the initial conditions. The optimization of the merit function is accomplished using the hybrid computing performances of PSO-ASA. The designed performance of ANNs-PSO-ASA is implemented for the numerical treatment of the VP-HBM dynamics by fluctuating the pulse shape adjustment terms, external forcing factor and damping coefficient with fixed ventricular contraction period. To perform the correctness of the present scheme, the obtained numerical results through the designed ANN-PSO-ASA will be compared with the Adams numerical method. The statistical investigations with larger dataset are provided using the "mean absolute deviation", "Theil's inequality coefficient" and "variance account for" operators to perform the applicability, reliability, and effectiveness of the designed ANNs-PSO-ASA scheme for solving the VP-HBM.