This paper presents a new method for designing digital recursive integrators and differentiators by using nonlinear stochastic global optimization based on particle swarm optimization (PSO). A modified PSO is used to optimize the unknown coefficients of second-order and third-order recursive integrators in order to obtain a better magnitude response closer to the ideal integrator. The coefficients of the state-of-theart available filters in the literature are chosen as initial points in the PSO process. By choosing a good starting point, convergence to an optimal solution is greatly facilitated. A dynamic modification of the fitness function used in the PSO process leads to the design of a set of digital integrators each optimized for a specific frequency range. Then, second-order and third-order recursive digital differentiators are designed by inverting and stabilizing the transfer functions of the designed recursive integrators. The obtained stabilized differentiators are further optimized using PSO to further improve their performance. The magnitude responses of the designed filters outperform the existing integrators and differentiators.Additionally Al-Alaoui in [2] proposed an approach to design a digital differentiator by inverting the transfer function of an integrator with similar specifications, stabilizing the resulting transfer function, and compensating the resulting change in magnitude because of the stabilization step that consisted of reflecting the resulting poles, the zeros of the integrator, that lie outside the unit circle at a radius r in the z-plane to 1/r. In [3], he obtained a class of second-order integrators and the corresponding differentiators by interpolating the Simpson and the trapezoidal integration rules and inverting the transfer functions of the obtained integrators.Ngo [4] employed Newton-Cotes integration method to design digital integrators. Tseng-Lee [5] used fractional delay to design digital integrators and Pei and Hsu [6] proposed employing fractional delay filters to obtain new differentiators. Gupta, Jain, and Kumar [7-9] developed digital integrators by interpolating some of the popular digital integration techniques. Varshney, Gupta, and Visweswaran [10] also used the interpolation method to design new digital differentiators. Some recent publications [11][12][13][14] in this field have also been reported in the literature.Recursive digital integrator's design can be formulated as an optimization problem where one searches for the optimal vector of filter coefficients to minimize a certain error function between the designed filter and the ideal integrator given by 1/jw. The previously mentioned methods for designing digital integrators suffer from poor accuracy because only a small fraction of the entire set of candidate vectors of filter coefficients is considered and thus are easily trapped into local optima on the error surface. Papamarkos and Chamzas [15] have used linear programming optimization method to propose the design of digital integrators. Lai, Lin, and ...