Abstract In this manuscript, we implement a spectral collocation method to find the solution of reaction–diffusion equation with some initial and boundary conditions. We approximate by using two-dimensional interpolating polynomial dependent Legendre–Gauss–Lobatto points. fully show that achieved solutions are convergent exact when number points increases. demonstrate capability efficiency prov...