Fast spin echo pulse sequences are adjusted to reduce, or eliminate image artifacts caused by Maxwell terms arising from the linear imaging gradients. The waveforms of the slice selection, phase encoding and readout gradients are adjusted in shape, size or position to eliminate or reduce the phase error caused by the spatially quadratic Maxwell terms.