A numerical solver for a nonlinear Fokker-Planck equation representation of neuronal network dynamics