In this paper we provide a simple MATLAB routine which computes the moving median with trend and seasonality. This approach is linear and for this reason has its disadvantages. So this routine can be improved by combining Monte-Carlo simulations, genetic algorithms simulations and wavelets decomposition.