Abstract We present the results of an exploratory study of the numerical stochastic perturbation theory (NSPT) applied to the four dimensional twisted Eguchi-Kawai (TEK) model. We employ a Kramers type algorithm based on the Generalized Hybrid Molecular Dynamics (GHMD) algorithm. We have computed the perturbative expansion of square Wilson loops up to O(g 8). https://safeersappliancers.shop/product-category/bosch-dhl785cgb-serie-6-built-in-70cm-4-speeds-c-rated-canopy-cooker-hood/