diff --git a/disba/_dispersion.py b/disba/_dispersion.py
index 8298b59..f572483 100644
--- a/disba/_dispersion.py
+++ b/disba/_dispersion.py
@@ -23,7 +23,7 @@ def __init__(
         velocity_p,
         velocity_s,
         density,
-        algorithm="fast-delta",
+        algorithm="dunkin",
         dc=0.005,
     ):
         """
@@ -39,7 +39,7 @@ def __init__(
             Layer S-wave velocity (in km/s).
         density : array_like
             Layer density (in g/cm3).
-        algorithm : str {'dunkin', 'fast-delta'}, optional, default 'fast-delta'
+        algorithm : str {'dunkin', 'fast-delta'}, optional, default 'dunkin'
             Algorithm to use for computation of Rayleigh-wave dispersion:
              - 'dunkin': Dunkin's matrix (adapted from surf96),
              - 'fast-delta': fast delta matrix (after Buchen and Ben-Hador, 1996).
@@ -98,7 +98,7 @@ def __init__(
         velocity_p,
         velocity_s,
         density,
-        algorithm="fast-delta",
+        algorithm="dunkin",
         dc=0.005,
         dt=0.025,
     ):
@@ -115,7 +115,7 @@ def __init__(
             Layer S-wave velocity (in km/s).
         density : array_like
             Layer density (in g/cm3).
-        algorithm : str {'dunkin', 'fast-delta'}, optional, default 'fast-delta'
+        algorithm : str {'dunkin', 'fast-delta'}, optional, default 'dunkin'
             Algorithm to use for computation of Rayleigh-wave dispersion:
              - 'dunkin': Dunkin's matrix (adapted from surf96),
              - 'fast-delta': fast delta matrix (after Buchen and Ben-Hador, 1996).
diff --git a/disba/_eigen.py b/disba/_eigen.py
index 782a87d..ea0469c 100644
--- a/disba/_eigen.py
+++ b/disba/_eigen.py
@@ -27,7 +27,7 @@ def __init__(
         velocity_p,
         velocity_s,
         density,
-        algorithm="fast-delta",
+        algorithm="dunkin",
         dc=0.005,
     ):
         """
@@ -43,7 +43,7 @@ def __init__(
             Layer S-wave velocity (in km/s).
         density : array_like
             Layer density (in g/cm3).
-        algorithm : str {'dunkin', 'fast-delta'}, optional, default 'fast-delta'
+        algorithm : str {'dunkin', 'fast-delta'}, optional, default 'dunkin'
             Algorithm to use for computation of Rayleigh-wave dispersion:
              - 'dunkin': Dunkin's matrix (adapted from surf96),
              - 'fast-delta': fast delta matrix (after Buchen and Ben-Hador, 1996).
diff --git a/disba/_ellipticity.py b/disba/_ellipticity.py
index 840ab24..95bcfea 100644
--- a/disba/_ellipticity.py
+++ b/disba/_ellipticity.py
@@ -25,7 +25,7 @@ def __init__(
         velocity_p,
         velocity_s,
         density,
-        algorithm="fast-delta",
+        algorithm="dunkin",
         dc=0.005,
     ):
         """
@@ -41,7 +41,7 @@ def __init__(
             Layer S-wave velocity (in km/s).
         density : array_like
             Layer density (in g/cm3).
-        algorithm : str {'dunkin', 'fast-delta'}, optional, default 'fast-delta'
+        algorithm : str {'dunkin', 'fast-delta'}, optional, default 'dunkin'
             Algorithm to use for computation of Rayleigh-wave dispersion:
              - 'dunkin': Dunkin's matrix (adapted from surf96),
              - 'fast-delta': fast delta matrix (after Buchen and Ben-Hador, 1996).
diff --git a/disba/_sensitivity.py b/disba/_sensitivity.py
index 7b074fc..b2a72ee 100644
--- a/disba/_sensitivity.py
+++ b/disba/_sensitivity.py
@@ -27,7 +27,7 @@ def __init__(
         velocity_p,
         velocity_s,
         density,
-        algorithm="fast-delta",
+        algorithm="dunkin",
         dc=0.005,
         dp=0.025,
     ):
@@ -44,7 +44,7 @@ def __init__(
             Layer S-wave velocity (in km/s).
         density : array_like
             Layer density (in g/cm3).
-        algorithm : str {'dunkin', 'fast-delta'}, optional, default 'fast-delta'
+        algorithm : str {'dunkin', 'fast-delta'}, optional, default 'dunkin'
             Algorithm to use for computation of Rayleigh-wave dispersion:
              - 'dunkin': Dunkin's matrix (adapted from surf96),
              - 'fast-delta': fast delta matrix (after Buchen and Ben-Hador, 1996).
@@ -110,7 +110,7 @@ def __init__(
         velocity_p,
         velocity_s,
         density,
-        algorithm="fast-delta",
+        algorithm="dunkin",
         dc=0.005,
         dt=0.025,
         dp=0.025,
@@ -128,7 +128,7 @@ def __init__(
             Layer S-wave velocity (in km/s).
         density : array_like
             Layer density (in g/cm3).
-        algorithm : str {'dunkin', 'fast-delta'}, optional, default 'fast-delta'
+        algorithm : str {'dunkin', 'fast-delta'}, optional, default 'dunkin'
             Algorithm to use for computation of Rayleigh-wave dispersion:
              - 'dunkin': Dunkin's matrix (adapted from surf96),
              - 'fast-delta': fast delta matrix (after Buchen and Ben-Hador, 1996).