A novel method is described for simulation of an electric stimulation of the nerve system subject to the rate of change of gradient fields. A gradient signal is filtered and a stimulation signal is derived, which is compared with a predetermined stimulation threshold value. An indicator signal is generated if the threshold value is exceeded. Therefore the time dependent and spatially dependent electric fields as defined by the scanning sequence and the gradient coil properties are calculated. A vector combination of said calculated electric field components from each gradient coil axis is performed, which results in a temporal diagram of the total electric field at various spatial locations within the gradient coil. The stimulation probability at each location from said temporal diagram and said stimulation signal is then calculated, and said stimulation probability is compared with the stimulation threshold value at each location within the gradient coil.