diff --git a/examples/shfqa_qubit_readout_measurement.md b/examples/shfqa_qubit_readout_measurement.md
index 8b944e54..58c9683e 100644
--- a/examples/shfqa_qubit_readout_measurement.md
+++ b/examples/shfqa_qubit_readout_measurement.md
@@ -51,7 +51,7 @@ device.qachannels[CHANNEL_INDEX].output.on(1)
 In this example the envelope of the readout pulses is a gaussian with a flat top. For each qubit, the envelope is then modulated at the qubit frequency. For illustrative purposes we assume that the frequencies of the qubits are equally spaced in the range [32 MHz, 230 MHz] relative to the center frequency of 5 GHz specified above.
 
 ```python
-from scipy.signal import gaussian
+from scipy.signal.windows import gaussian
 import numpy as np
 
 # Define the parameters for the readout pulses
diff --git a/examples/shfqa_qubit_readout_weights.md b/examples/shfqa_qubit_readout_weights.md
index d54477cb..d0710196 100644
--- a/examples/shfqa_qubit_readout_weights.md
+++ b/examples/shfqa_qubit_readout_weights.md
@@ -73,7 +73,7 @@ Generate the readout pulses with a flat top gaussian envelope. For each qubit th
 
 ```python
 
-from scipy.signal import gaussian
+from scipy.signal.windows import gaussian
 import numpy as np
 
 # Define the parameters for the readout pulses
diff --git a/examples/shfqa_shfqc_power_spectral_density.md b/examples/shfqa_shfqc_power_spectral_density.md
index 8f305145..ac50d015 100644
--- a/examples/shfqa_shfqc_power_spectral_density.md
+++ b/examples/shfqa_shfqc_power_spectral_density.md
@@ -170,7 +170,7 @@ else:
 # Generate an interesting test waveform
 from zhinst.toolkit import Waveforms
 import numpy as np
-from scipy import signal
+from scipy.signal.windows import gaussian
 from zhinst.utils.shfqa import SHFQA_SAMPLING_FREQUENCY
 
 # Modulate a Gaussian with complex sinusoidal
@@ -180,7 +180,7 @@ modulation_freq = 100e6
 waveform_length = np.min(128, int(sweeper.average.integration_time() / SHFQA_SAMPLING_FREQUENCY))
 time_axis = np.arange(waveform_length) / SHFQA_SAMPLING_FREQUENCY
 test_signal = np.exp(1j * 2 * np.pi * modulation_freq * time_axis)
-test_signal *= signal.gaussian(waveform_length, std=waveform_length/8)
+test_signal *= gaussian(waveform_length, std=waveform_length/8)
 
 # Upload waveform to the device
 waveforms = Waveforms()
diff --git a/examples/shfqa_sweeper.md b/examples/shfqa_sweeper.md
index 898a08d6..60ebf86e 100644
--- a/examples/shfqa_sweeper.md
+++ b/examples/shfqa_sweeper.md
@@ -117,7 +117,7 @@ In pulsed spectroscopy the resonator is probed with a signal consisting of an en
 ### Create the envelope
 
 ```python
-from scipy.signal import gaussian
+from scipy.signal.windows import gaussian
 import numpy as np
 
 SHFQA_SAMPLING_FREQUENCY = 2e9
diff --git a/examples/shfqc_helper.py b/examples/shfqc_helper.py
index c25f6063..956296fa 100644
--- a/examples/shfqc_helper.py
+++ b/examples/shfqc_helper.py
@@ -45,7 +45,7 @@ def generate_flat_top_gaussian(
                 unit circle."
         )
 
-    from scipy.signal import gaussian
+    from scipy.signal.windows import gaussian
 
     rise_fall_len = int(rise_fall_time * sampling_rate)
     pulse_len = int(pulse_duration * sampling_rate)