diff --git a/mpython/sparse_array.py b/mpython/sparse_array.py
index bcbc927..a2dae2c 100644
--- a/mpython/sparse_array.py
+++ b/mpython/sparse_array.py
@@ -8,21 +8,21 @@
 
 if sparse:
 
-    class WrappedSparseArray(sparse.sparray, AnyWrappedArray):
+    class WrappedSparseArray(sparse.spmatrix, AnyWrappedArray):
         """Base class for sparse arrays."""
 
         def to_dense(self) -> "Array":
             return Array.from_any(self.todense())
 
-    class SparseArray(sparse.csc_array, _SparseMixin, WrappedSparseArray):
+    class SparseArray(sparse.csc_matrix, _SparseMixin, WrappedSparseArray):
         """
-        Matlab sparse arrays (scipy.sparse backend).
+        Matlab sparse matrices (scipy.sparse backend).
 
         ```python
         # Instantiate from size
-        SparseArray(N, M, ...)
-        SparseArray([N, M, ...])
-        SparseArray.from_shape([N, M, ...])
+        SparseArray(N, M)
+        SparseArray([N, M])
+        SparseArray.from_shape([N, M])
 
         # Instantiate from existing sparse or dense array
         SparseArray(other_array)
@@ -45,7 +45,7 @@ def __init__(self, *args, **kwargs) -> None:
                 ndim = len(arg)
                 return super().__init__(([], [[]] * ndim), shape=arg, **kwargs)
             else:
-                if not isinstance(arg, (np.ndarray, sparse.sparray)):
+                if not isinstance(arg, (np.ndarray, sparse.spmatrix)):
                     arg = np.asanyarray(arg)
                 return super().__init__(arg, **kwargs)
 
@@ -71,7 +71,7 @@ def from_coo(cls, values, indices, shape=None, **kw) -> "SparseArray":
                 New array.
             """
             indices = np.asarray(indices)
-            coo = sparse.coo_array((values, indices), shape=shape, **kw)
+            coo = sparse.coo_matrix((values, indices), shape=shape, **kw)
             return cls.from_any(coo)
 
         @classmethod
@@ -124,7 +124,7 @@ def from_any(cls, other, **kwargs) -> "SparseArray":
             """
             copy = kwargs.pop("copy", None)
             inp = other
-            if not isinstance(other, sparse.sparray):
+            if not isinstance(other, sparse.spmatrix):
                 other = np.asanyarray(other, **kwargs)
             other = cls(other, **kwargs)
             other = _spcopy_if_needed(other, inp, copy)
diff --git a/mpython/utils.py b/mpython/utils.py
index ebd0370..c3d42d2 100644
--- a/mpython/utils.py
+++ b/mpython/utils.py
@@ -108,7 +108,7 @@ def _spcopy_if_needed(out, inp, copy=None):
     """Fallback implementation for asarray(*, copy: bool)"""
     if (
         out is not None
-        and isinstance(inp, sparse.sparray)
+        and isinstance(inp, sparse.spmatrix)
         and out.data.data != inp.data.data
     ):
         if copy:
diff --git a/pyproject.toml b/pyproject.toml
index bf45210..f1c4713 100644
--- a/pyproject.toml
+++ b/pyproject.toml
@@ -32,6 +32,9 @@ dependencies = [
     "numpy"
 ]
 
+[project.optional-dependencies]
+sparse = ["scipy>=0.7"]
+
 [tool.setuptools.dynamic]
 version = {attr = "mpython._version.__version__"}