This project investigated the behavior of a nonlinear harmonic oscillator solver and explained the observed loss of accuracy under certain conditions. It extended a linear harmonic oscillator solver to a nonlinear counterpart using the model 'Method of Manufactured Solutions'.

Testing the implementation of numerical methods for solving the convection diffusion problem with variable coefficients and Neumann boundary conditions