Keyword for solve_2d for alpha and dt.
Import solve_2d as part of package.
Adjust test for non-unit time step.
Use fixed boundaries for solve_2d.
Removed dt keyword.
Do not reset boundary values.
Rename to reflect we're solving the heat equation.