The use of theta method in waveform relaxation iterations: convergence analysis