diff --git a/miscellaneous/rankine.py b/miscellaneous/rankine.py
index 4d6706f6..614ee861 100644
--- a/miscellaneous/rankine.py
+++ b/miscellaneous/rankine.py
@@ -118,7 +118,7 @@ def plot_rankine( vlsvReader, point1, point2 ):
    :param point1: The first point of interest along the bow-shock
    :param point2: The second point of interest along the bow-shock
 
-   :returns: pylab figure
+   :returns: matplotlib figure
    '''
    # Get spatial grid sizes:
    xcells = (int)(vlsvReader.read_parameter("xcells_ini"))
diff --git a/pyCalculations/timeevolution.py b/pyCalculations/timeevolution.py
index 26837ac6..1fb9d905 100644
--- a/pyCalculations/timeevolution.py
+++ b/pyCalculations/timeevolution.py
@@ -38,7 +38,7 @@ def cell_time_evolution( vlsvReader_list, variables, cellids, units="" ):
 
        .. code-block:: python
 
-          import pytools as pt; import pylab as pl
+          import pytools as pt; import matplotlib.pyplot as pl
           # Example of usage:
           time_data = pt.calculations.cell_time_evolution( vlsvReader_list=[VlsvReader("bulk.000.vlsv"), VlsvReader("bulk.001.vlsv"), VlsvReader("bulk.002.vlsv")], variables=["rho", "Pressure", "B"], cellids=[2,4], units=["N", "Pascal", "T"] )
 
diff --git a/pyMayaVi/particlepusherinterface.py b/pyMayaVi/particlepusherinterface.py
index 0e779c80..97e4e2ac 100644
--- a/pyMayaVi/particlepusherinterface.py
+++ b/pyMayaVi/particlepusherinterface.py
@@ -32,7 +32,7 @@
 from mayavi.core.ui.mayavi_scene import MayaviScene
 import vlsvfile
 from numpy import mgrid, empty, sin, pi, ravel
-import pylab as pl
+
 from tvtk.api import tvtk
 import traits.api
 import mayavi.api
diff --git a/scripts/magnetopause2d.py b/scripts/magnetopause2d.py
index 62174652..739125e5 100644
--- a/scripts/magnetopause2d.py
+++ b/scripts/magnetopause2d.py
@@ -2,7 +2,6 @@
 Finds the magnetopause position by tracing steamines of the plasma flow for two-dimensional Vlasiator runs. Needs the yt package.
 '''
 
-import matplotlib.pylab as plt
 import numpy as np
 import pytools as pt
 import plot_colormap
diff --git a/scripts/magnetopause3d.py b/scripts/magnetopause3d.py
index a6aa0a9f..a949ae98 100644
--- a/scripts/magnetopause3d.py
+++ b/scripts/magnetopause3d.py
@@ -2,7 +2,7 @@
 Finds the magnetopause position by tracing steamines of the plasma flow for three-dimensional Vlasiator runs. Needs the yt package.
 '''
 from pyCalculations import ids3d
-import matplotlib.pylab as plt
+import matplotlib.pyplot as plt
 import numpy as np
 import pytools as pt
 import plot_colormap3dslice