New derivation of equations, including weak form statement. mshr and spherical mesh only work on python 3.6 for now. Changing initial conditions, time step, grid size can all influence stability of algorithm.