diff --git a/README.md b/README.md
index f0b6f51e..97bd9e82 100644
--- a/README.md
+++ b/README.md
@@ -25,6 +25,7 @@ Here’s a basic example to get you started with QuAM:
 
 ```python
 from quam.components import BasicQuAM, SingleChannel, pulses
+from qm import qua
 
 # Create a root-level QuAM instance
 machine = BasicQuAM()
@@ -34,14 +35,14 @@ channel = SingleChannel(opx_output=("con1", 1))
 machine.channels["output"] = channel
 
 # Add a Gaussian pulse to the channel
-channel.operations["gaussian"] = pulses.Gaussian(
+channel.operations["gaussian"] = pulses.GaussianPulse(
     length=100,  # Pulse length in ns
     amplitude=0.5,  # Peak amplitude of Gaussian pulse
     sigma=20,  # Standard deviation of Guassian pulse
 )
 
 # Play the Gaussian pulse on the channel within a QUA program
-with program() as prog:
+with qua.program() as prog:
     channel.play("gaussian")
 
 # Generate the QUA configuration from QuAM
diff --git a/docs/index.md b/docs/index.md
index 3c7029ef..812146d5 100644
--- a/docs/index.md
+++ b/docs/index.md
@@ -19,6 +19,7 @@ QuAM is not just a tool but a gateway to streamlined and efficient quantum compu
 
 ```python
 from quam.components import *
+from qm import qua
 
 # Create a root-level QuAM instance
 machine = BasicQuAM()
@@ -28,12 +29,12 @@ channel = SingleChannel(opx_output=("con1", 1))
 machine.channels["output"] = channel
 
 # Add a Gaussian pulse to the channel
-channel.operations["gaussian"] = pulses.Gaussian(
+channel.operations["gaussian"] = pulses.GaussianPulse(
     length=100, amplitude=0.5, sigma=20
 )
 
 # Play the Gaussian pulse within a QUA program
-with program() as prog:
+with qua.program() as prog:
     channel.play("gaussian")
 
 # Generate the QUA configuration from QuAM