Create 1D_anharmonic_oscillator.cpp