The Bloch equations are applied to describe the response of the two-level system, irradiated with bichromatic microwave field under condition of electron paramagnetic resonance in the wide range of saturation condition and relaxation times. The generation of I:he intermodulation effects in such systems has been demonstrated. Numerical procedures of solving a set of stiff, nonlinear differential equations are applied. The line shapes of the generated signals have been simulated for different microwave amplitudes and relaxation rates. The dependence of the intermodulation signals on omega(2) -omega(1), the frequency difference between two incident microwave: fields, is studied. The influence of a significant saturation factor on this dependence is demonstrated and explained. Then an estimation of the relaxation rates from the intermodulation sidebands is examined on the basis of the classical model. In the case when the frequency separation between two alternating fields is less than relaxation rates, analytical formulas for the intermodulation signals have been found.