Approximation by bivariate Szász–Mirakyan operators preserving e-2(p₁ + p₂)

The present paper is dedicated to the modification of the bivariate generalized Sz?sz?Mirakyan operators while preserving the exponential functions exp(2,2) where exp(τ_1, τ_2) = e-τ_1 p_1 - τ_2 p_2 , τ_1, τ_2 ∈ R_0^+ , and p_1, p_2 ≥ 0 . We thoroughly investigate the weighted approximation properties and also obtain the convergence rate for these operators by utilizing a weighted modulus of continuity. Additionally, we delve into the order of approximation by investigating local approximation results through Peetre's K -functional. Furthermore, we present the GBS (Generalized Boolean Sum) operators of Sz?sz?Mirakyan operators and obtain the order of approximation in terms of the Lipschitz class of B?gel continuous functions and the mixed modulus of smoothness. In order to enhance our theoretical findings and effectively showcase the efficiency of our developed operators, we have included a wide range of numerical examples using various values.