In this paper, a bivariate spectral quasilinearization method is used to solve the highly nonlinear two dimensional Bratu problem. The two dimensional Bratu problem is also solved using the Chebyshev spectral collocation ...