Numerical methods for the extended Fisher-Kolmogorov (EFK) equation