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.