[eng] This paper implements the approach introduced by MacKinnon (J Bus Econ Stat 12:167-176, 1994, J Appl Econom 11:601-618, 1996) to estimate the response surface of the test statistics of seasonal unit root tests with OLS and GLS detrending for quarterly and monthly time series. The Gauss code that is available in the supplementary material of the paper produces p values for five test statistics depending on the sample size, deterministic terms and frequency of the data. A comparison with previous studies is undertaken, and an empirical example using airport passenger arrivals to a tourist destination is carried out. Quantile function coefficients are reported for simple computation of critical values for tests at 1, 5 and 10 % significance levels.