diff --git a/setup.py b/setup.py
index 894efc8..dce4083 100644
--- a/setup.py
+++ b/setup.py
@@ -28,7 +28,7 @@
           'numpy',
           'h5py',
           'pandas',
-          'matplotlib'
+          'matplotlib>=3.5'
       ],
     include_package_data=True,
     options={"bdist_wheel": {"universal": True}}
diff --git a/src/abg_python/pfh_colormaps.py b/src/abg_python/pfh_colormaps.py
index a3423f1..e04d9d6 100644
--- a/src/abg_python/pfh_colormaps.py
+++ b/src/abg_python/pfh_colormaps.py
@@ -1,4 +1,5 @@
-import matplotlib
+from matplotlib.colors import LinearSegmentedColormap
+import matplotlib as mpl
 
 def load_my_custom_color_tables():
     fna='heat_red'
@@ -7,8 +8,8 @@ def load_my_custom_color_tables():
     'green': ((0., 0., 0.),(0.71, 0.000000, 0.000000),\
             (1., 1.000000, 1.000000),(1.0, 1.0, 1.0)),\
     'blue':  ((0., 0., 0.),(0.71, 0.000000, 0.000000),(1.0, 1.0, 1.0))}
-    my_cmap = matplotlib.colors.LinearSegmentedColormap(fna,cdict_tmp,256);
-    matplotlib.cm.register_cmap(name=fna,cmap=my_cmap);
+    my_cmap = LinearSegmentedColormap(fna,cdict_tmp,256)
+    mpl.colormaps.register(name=fna,cmap=my_cmap)
 
     fna='heat_blue'
     cdict_tmp={\
@@ -16,8 +17,8 @@ def load_my_custom_color_tables():
     'green': ((0., 0., 0.),(0.365079, 0.000000, 0.000000),\
             (0.746032, 1.000000, 1.000000),(1.0, 1.0, 1.0)),\
     'red':  ((0., 0., 0.),(0.746032, 0.000000, 0.000000),(1.0, 1.0, 1.0))}
-    my_cmap = matplotlib.colors.LinearSegmentedColormap(fna,cdict_tmp,256);
-    matplotlib.cm.register_cmap(name=fna,cmap=my_cmap);
+    my_cmap = LinearSegmentedColormap(fna,cdict_tmp,256)
+    mpl.colormaps.register(name=fna,cmap=my_cmap)
 
     fna='heat_green'
     cdict_tmp={\
@@ -25,8 +26,8 @@ def load_my_custom_color_tables():
     'red': ((0., 0., 0.),(0.365079, 0.000000, 0.000000),\
             (1.000000, 1.000000, 1.000000),(1.0, 1.0, 1.0)),\
     'blue':  ((0., 0., 0.),(0.1, 0.000000, 0.000000),(1.0, 1.0, 1.0))}
-    my_cmap = matplotlib.colors.LinearSegmentedColormap(fna,cdict_tmp,256);
-    matplotlib.cm.register_cmap(name=fna,cmap=my_cmap);
+    my_cmap = LinearSegmentedColormap(fna,cdict_tmp,256)
+    mpl.colormaps.register(name=fna,cmap=my_cmap)
 
     fna='heat_redyellow'
     cdict_tmp={\
@@ -34,16 +35,16 @@ def load_my_custom_color_tables():
     'green': ((0., 0., 0.),(0.365079, 0.000000, 0.000000),\
             (0.746032, 1.000000, 1.000000),(1.0, 1.0, 1.0)),\
     'blue':  ((0., 0., 0.),(0.746032, 0.000000, 0.000000),(1.0, 1.0, 1.0))}
-    my_cmap = matplotlib.colors.LinearSegmentedColormap(fna,cdict_tmp,256);
-    matplotlib.cm.register_cmap(name=fna,cmap=my_cmap);
+    my_cmap = LinearSegmentedColormap(fna,cdict_tmp,256)
+    mpl.colormaps.register(name=fna,cmap=my_cmap)
 
     fna='heat_yellow'
     cdict_tmp={\
     'red':   ((0., 0.0416, 0.0416),(0.365079, 1.000000, 1.000000),(1.0, 1.0, 1.0)),\
     'green':   ((0., 0.0416, 0.0416),(0.365079, 1.000000, 1.000000),(1.0, 1.0, 1.0)),\
     'blue':  ((0., 0., 0.),(0.746032, 0.000000, 0.000000),(1.0, 1.0, 1.0))}
-    my_cmap = matplotlib.colors.LinearSegmentedColormap(fna,cdict_tmp,256);
-    matplotlib.cm.register_cmap(name=fna,cmap=my_cmap);
+    my_cmap = LinearSegmentedColormap(fna,cdict_tmp,256)
+    mpl.colormaps.register(name=fna,cmap=my_cmap)
 
     fna='heat_purple'
     cdict_tmp={\
@@ -51,16 +52,16 @@ def load_my_custom_color_tables():
     'green': ((0., 0., 0.),(0.565079, 0.000000, 0.000000),\
             (0.946032, 1.000000, 1.000000),(1.0, 1.0, 1.0)),\
     'blue':   ((0., 0.0416, 0.0416),(0.565079, 1.000000, 1.000000),(1.0, 1.0, 1.0))}
-    my_cmap = matplotlib.colors.LinearSegmentedColormap(fna,cdict_tmp,256);
-    matplotlib.cm.register_cmap(name=fna,cmap=my_cmap);
+    my_cmap = LinearSegmentedColormap(fna,cdict_tmp,256)
+    mpl.colormaps.register(name=fna,cmap=my_cmap)
 
     fna='heat_orange'
     cdict_tmp={\
     'red':   ((0., 0.0416, 0.0416),(0.365079, 1.000000, 1.000000),(1.0, 1.0, 1.0)),\
     'green':   ((0., 0.0416, 0.0416),(0.365079, 0.400000, 0.400000),(1.0, 1.0, 1.0)),\
     'blue':  ((0., 0., 0.),(0.746032, 0.000000, 0.000000),(1.0, 1.0, 1.0))}
-    my_cmap = matplotlib.colors.LinearSegmentedColormap(fna,cdict_tmp,256);
-    matplotlib.cm.register_cmap(name=fna,cmap=my_cmap);
+    my_cmap = LinearSegmentedColormap(fna,cdict_tmp,256)
+    mpl.colormaps.register(name=fna,cmap=my_cmap)
 
     """
     fna='rainbow'
@@ -80,6 +81,6 @@ def load_my_custom_color_tables():
                        (0.4, 0.0, 0.0),
                        (0.6, 20/256., 20/256.),
                        (1.0, 0.0, 0.0))}
-    my_cmap = matplotlib.colors.LinearSegmentedColormap(fna,cdict_tmp,256);
-    matplotlib.cm.register_cmap(name=fna,cmap=my_cmap);
+    my_cmap = LinearSegmentedColormap(fna,cdict_tmp,256)
+    mpl.colormaps.register(name=fna,cmap=my_cmap)
     """