diff --git a/docs/source/graphene_example.md b/docs/source/graphene_example.md
index f53045ebd61496885689b7cc8d9643fb3fcbbabf..61f6eb68de656ca47de6cbd2901b4f9dbaac3668 100644
--- a/docs/source/graphene_example.md
+++ b/docs/source/graphene_example.md
@@ -55,7 +55,7 @@ h_int = kwant_utils.builder_to_tb(int_builder, params)
 model = model.Model(h_0, h_int, filling=2)
 ```
 
-To start the meanfield calculation we also need a starting guess. We will use our random guess generator for this. It creates a random Hermitian hopping dictionary based on the hopping keys provided and the number of degrees of freedom specified. As we don't expect the mean-field solution to contain terms more than the hoppings from the interacting part, we can use the hopping keys from the interacting part. We will use the same numbe of degrees as freedom as both the non-interacting and interacting part, so that they match.
+To start the meanfield calculation we also need a starting guess. We will use our random guess generator for this. It creates a random Hermitian hopping dictionary based on the hopping keys provided and the number of degrees of freedom specified. As we don't expect the mean-field solution to contain terms more than the hoppings from the interacting part, we can use the hopping keys from the interacting part. We will use the same number of degrees as freedom as both the non-interacting and interacting part, so that they match.
 
 ```{code-cell} ipython3
 guess = tb.utils.generate_guess(frozenset(h_int), len(list(h_0.values())[0]))
@@ -67,4 +67,6 @@ After we have defined the guess, we feed it together with the model into the mea
 
 ## Creating a phase diagram
 
-We can now create a phase diagram by varying the Hubbard interaction strength and the nearest neighbor interaction strength. We will use the same hopping dictionary for the non-interacting part as before. We will vary the onsite Hubbard interaction strength from
+We can now create a phase diagram of the gap of the interacting solution. We will use the same hopping dictionary for the non-interacting part as before. We will vary the onsite Hubbard interactio $U$ strength from $0$ to $2$ and the nearest neighbor interaction strength $V$ from $0$ to $1.5$.
+
+```{code-cell} ipython3
diff --git a/docs/source/index.md b/docs/source/index.md
index 3df6663c32a93642efeb4d48c9c4e39058b0dcfc..19853c7ac49277dc216a0401ce3d77cb2bf6febe 100644
--- a/docs/source/index.md
+++ b/docs/source/index.md
@@ -19,6 +19,7 @@ kernelspec:
 :caption: Tutorials
 
 mf_notes.md
+graphene_example.md
 ```
 
 ## What is pymf?
diff --git a/pymf/__init__.py b/pymf/__init__.py
index 680851766e11df163131091ac162b4de364c0d22..a0f2b6e8598cdbc623b27aaeffdd5555fbcb8f96 100644
--- a/pymf/__init__.py
+++ b/pymf/__init__.py
@@ -7,7 +7,6 @@ except ImportError:
     __version_tuple__ = (0, 0, "unknown", "unknown")
 
 from .mf import construct_density_matrix
-from . import mf
 
 __all__ = [
     "construct_density_matrix",