You have a couple of options:
- Use the PDE toolbox in MATLAB to solve it. Take a look at example 1 here.
- Assuming you know your initial conditions and boundary conditions for your equation, ie what the system looks like at say r=0 and r=R at t=0. Then you could convert all the 1st derivatives (the temporal part) using a forward finite difference and then for the 2nd derivatives (spatial part) use the central difference method. Here's some nice documentation I found online that explains it. Here's also a video that shows this finite difference method but for the classic heat equation instead.
