diff --git a/demo/demo_nonlinear-naghdi-clamped-semicylinder.py b/demo/demo_nonlinear-naghdi-clamped-semicylinder.py
index 9ccc1e1..2adda06 100644
--- a/demo/demo_nonlinear-naghdi-clamped-semicylinder.py
+++ b/demo/demo_nonlinear-naghdi-clamped-semicylinder.py
@@ -455,7 +455,7 @@ def gamma(F, d):
 # Calculate the factor alpha as a function of the mesh size h
 h = ufl.CellDiameter(mesh)
 alpha_FS = functionspace(mesh, element("DG", cell, 0))
-alpha_expr = Expression(t**2 / h**2, alpha_FS.element.interpolation_points())
+alpha_expr = Expression(t**2 / h**2, alpha_FS.element.interpolation_points)
 alpha = Function(alpha_FS)
 alpha.interpolate(alpha_expr)
 
@@ -693,7 +693,7 @@ def F(self, x: PETSc.Vec, b: PETSc.Vec) -> None:
 
 # Interpolate phi in the [P2]³ Space
 phi_FS = functionspace(mesh, blocked_element(P2, shape=(3,)))
-phi_expr = Expression(phi0_ufl + u_P2B3, phi_FS.element.interpolation_points())
+phi_expr = Expression(phi0_ufl + u_P2B3, phi_FS.element.interpolation_points)
 phi_func = Function(phi_FS)
 phi_func.interpolate(phi_expr)