We study the interacting electromagnetic shock waves with nonconstant profiles in general relativity. It is shown that by modifying the metric functions of the Bell-Szekeres solution, such solutions can be obtained.