diff --git a/README.md b/README.md
index 8ad7ff11..13e42efe 100644
--- a/README.md
+++ b/README.md
@@ -54,13 +54,13 @@ import skgmsh as sg
 We can define the surface using PyVista.
 
 ```python
-edge_source = sg.Polygon(n_sides=4, radius=8, fill=False)
+source = sg.Polygon(n_sides=4, radius=8, fill=False)
 ```
 
 We can then generate a 2D mesh.
 
 ```python
-mesh = sg.frontal_delaunay_2d(edge_source, target_sizes=2.0)
+mesh = source.frontal_delaunay_2d(edge_source=source, target_sizes=2.0)
 ```
 
 To visualize the model, we can use PyVista.
@@ -75,7 +75,7 @@ _ = plotter.add_mesh(
     lighting=False,
     edge_color="gray",
 )
-_ = plotter.add_mesh(edge_source, show_edges=True, line_width=4, color="gray")
+_ = plotter.add_mesh(source, show_edges=True, line_width=4, color="gray")
 plotter.show(cpos="xy")
 ```
 
@@ -86,8 +86,8 @@ plotter.show(cpos="xy")
 We can also generate a 3D mesh.
 
 ```python
-edge_source = sg.Cube()
-mesh = sg.delaunay_3d(edge_source, target_sizes=0.2)
+source = sg.Cube()
+mesh = source.delaunay_3d(edge_source=source, target_sizes=0.2)
 ```
 
 ```python
@@ -100,7 +100,7 @@ _ = plotter.add_mesh(
     lighting=False,
     edge_color="gray",
 )
-_ = plotter.add_mesh(edge_source.extract_all_edges(), line_width=4, color="gray")
+_ = plotter.add_mesh(source.extract_all_edges(), line_width=4, color="gray")
 _ = plotter.add_box_axes()
 plotter.show()
 ```