In this paper, an efficient numerical scheme is settled for solving two-dimensional Bratu–Gelfand problem, namely Hybrid Orthonormal Bernstein and Block-Pulse functions wavelet (HOBW) is presented for boundary value problems administered by nonlinear partial differential equations which effectively combines the Orthonormal Bernstein, Block-Pulse functions and the generalized wavelet. Operational Matrix of integration is utilized to provide an approximate result of the BG problems. By using the Operational Matrix, differentiation is changed to the nonlinear system of equations which can be disbanded via the Newton Raphson technique. As per our concentrated inquiry there is no exact solution of the problem and can solve the problem with higher accuracy than the methodologies used to solve this problem. The result is plotted for different values of then compared with the previous numerical results obtained.