This seems like a homework question. if you can provide the code you've written to solve the problem and ask a specific question about where you're facing difficulty, we may be able to provide some guidance. Please see the guidelines for asking homework questions on MATLAB Answers found here.
For a start, you can look into the pdepe function, to solve 1-D parabolic and elliptic PDEs, PDE toolbox, and this file exchange submission, which might give you some insight.
The submission, CFD-Navier-Stokes, is one of the several submissions in MATLAB File Exchange on MATLAB Central which is a forum for our product users to interact, exchange information and knowledge, without MathWorks' involvement. Feel free to contact the author of this submission directly for specific questions about the implementation.
