diff --git a/Examples/varia/AccessingSimulationResults.py b/Examples/varia/AccessingSimulationResults.py
index a49b09a389bedfbed6794a4b10203f466e17e29a..930cb06c8e94c52f23a15dbc6caae8d1ddfbb7ea 100755
--- a/Examples/varia/AccessingSimulationResults.py
+++ b/Examples/varia/AccessingSimulationResults.py
@@ -12,37 +12,8 @@ from matplotlib import pyplot as plt
 
 
 def get_sample():
-    """
-    Returns a sample with uncorrelated cylinders on a substrate.
-    """
-
-    # Define materials
-    material_Particle = ba.HomogeneousMaterial("Particle", 0.0006, 2e-08)
-    material_Substrate = ba.HomogeneousMaterial("Substrate", 6e-06, 2e-08)
-    material_Vacuum = ba.HomogeneousMaterial("Vacuum", 0, 0)
-
-    # Define form factors
-    ff = ba.FormFactorCylinder(5*nm, 5*nm)
-
-    # Define particles
-    particle = ba.Particle(material_Particle, ff)
-
-    # Define particle layouts
-    layout = ba.ParticleLayout()
-    layout.addParticle(particle)
-    layout.setTotalParticleSurfaceDensity(0.01)
-
-    # Define layers
-    layer_1 = ba.Layer(material_Vacuum)
-    layer_1.addLayout(layout)
-    layer_2 = ba.Layer(material_Substrate)
-
-    # Define sample
-    sample = ba.MultiLayer()
-    sample.addLayer(layer_1)
-    sample.addLayer(layer_2)
-
-    return sample
+    from bornagain import std_samples
+    return std_samples.cylinders()
 
 
 def get_simulation(sample):
diff --git a/Examples/varia/AxesInDifferentUnits.py b/Examples/varia/AxesInDifferentUnits.py
index 7212c4aeca90562c073dba0a924f43cd2ee54561..7f05efb29adb83665eb4fbfb869b15baba3d71c1 100755
--- a/Examples/varia/AxesInDifferentUnits.py
+++ b/Examples/varia/AxesInDifferentUnits.py
@@ -10,37 +10,8 @@ from matplotlib import rcParams
 
 
 def get_sample():
-    """
-    Returns a sample with uncorrelated cylinders on a substrate.
-    """
-
-    # Define materials
-    material_Particle = ba.HomogeneousMaterial("Particle", 0.0006, 2e-08)
-    material_Substrate = ba.HomogeneousMaterial("Substrate", 6e-06, 2e-08)
-    material_Vacuum = ba.HomogeneousMaterial("Vacuum", 0, 0)
-
-    # Define form factors
-    ff = ba.FormFactorCylinder(5*nm, 5*nm)
-
-    # Define particles
-    particle = ba.Particle(material_Particle, ff)
-
-    # Define particle layouts
-    layout = ba.ParticleLayout()
-    layout.addParticle(particle)
-    layout.setTotalParticleSurfaceDensity(0.01)
-
-    # Define layers
-    layer_1 = ba.Layer(material_Vacuum)
-    layer_1.addLayout(layout)
-    layer_2 = ba.Layer(material_Substrate)
-
-    # Define sample
-    sample = ba.MultiLayer()
-    sample.addLayer(layer_1)
-    sample.addLayer(layer_2)
-
-    return sample
+    from bornagain import std_samples
+    return std_samples.cylinders()
 
 
 def get_simulation(sample):