The aim of the paper is to present a procedure for generating service loading for fatigue tests of materials and structures. The generated loading characterizes desired functions of probability distribution and autocorrelation. The proposed numerical procedure uses MATLAB toolboxes and consists of three steps: (a) generation of a sequence of real numbers with the desired autocorrelation function and with any probability distribution function; (b) generation of loading history with the desired probability distribution function; (c) rearrangement of loading history (mentioned in item b) based on a sequence of real numbers with the desired correlation (mentioned in item a).