In this paper, we present a mathematical model of brain tumor. This model is an extension of a simple twodimensional mathematical model of glioma growth and diffusion which is derived from fractional operator in terms of Caputo which is called the fractional Burgess equations (FBEs). To obtain a solution for this model, a numerical technique is presented which is based on operational matrix. First, we assume the solution of the problem under the study is as an expansion of the Bernoulli polynomials. Then with combination of the operational matrix based on the Bernoulli polynomials and collocation method, the problem under the study is changed to a system of nonlinear algebraic equations. Finally, the proposed technique is simulated and tested on three types of the FBEs to confirm the superiority and accuracy
