Create solution_diffusion.reaction.equation.py

Requires side_conditions.py, physical_constants.py, and make_plot.py