diff --git a/festim/boundary_conditions/dirichlet_bc.py b/festim/boundary_conditions/dirichlet_bc.py
index cd334c483..aea5c2753 100644
--- a/festim/boundary_conditions/dirichlet_bc.py
+++ b/festim/boundary_conditions/dirichlet_bc.py
@@ -119,19 +119,6 @@ def create_value(
                 )
                 self.value_fenics.interpolate(self.bc_expr)
 
-    def create_formulation(self, dofs, function_space):
-        """Applies the boundary condition
-        Args:
-            dofs (numpy.ndarray): the degrees of freedom of surface facets
-            function_space (dolfinx.fem.FunctionSpace): the function space
-        """
-        form = fem.dirichletbc(
-            value=self.value_fenics,
-            dofs=dofs,
-            V=function_space,
-        )
-        return form
-
     def update(self, t):
         """Updates the boundary condition value
 
diff --git a/festim/hydrogen_transport_problem.py b/festim/hydrogen_transport_problem.py
index a2880619d..f7a1db8b9 100644
--- a/festim/hydrogen_transport_problem.py
+++ b/festim/hydrogen_transport_problem.py
@@ -284,7 +284,6 @@ def create_dirichletbc_form(self, bc):
         )
 
         if (len(self.species) == 1) and (isinstance(bc.value_fenics, (fem.Function))):
-            # todo remove DirichletBC.create_form
             return fem.dirichletbc(
                 value=bc.value_fenics,
                 dofs=bc_dofs,