My early work involved using time-stepping methods to solve differential equations and P.D.E.s.

I began by solving the heat equation (for a scalar) with periodic boundary conditions in one dimension.

I then modified my program to 2d then 3d. The next step was to adapt my program to a vectorial form.

I now have written a program that solves the 3d Heat Equation (aka the Diffusion equation) for a vector

Here's the code for this program.