“…A new bi-cubic spline collocation method of fourth order accuracy has been designed by Singh and Singh [57] for a linear EPDE with DBCs. Most recently, new high accuracy approximations in exponential form for solving two-point BVPs on a variable mesh, and 2D nonlinear EPDEs on an unequal uniform mesh have been reported in [51] , [52] , [53] , [54] , [55] . Mohanty and Kumar [56] , and Priyadarshini and Mohanty [ 47 , 49 ] have proposed new high accuracy numerical algorithms for 2D quasilinear EPDEs using unequal mesh.…”