Predicting primary energy consumption using NDGM(1,1,k,c) model with Simpson formula