diff --git a/examples/diffusion/steadyState/mesh20x20/modifiedMeshInput.py b/examples/diffusion/steadyState/mesh20x20/modifiedMeshInput.py
index 5ab0688167..4927ca976c 100755
--- a/examples/diffusion/steadyState/mesh20x20/modifiedMeshInput.py
+++ b/examples/diffusion/steadyState/mesh20x20/modifiedMeshInput.py
@@ -46,7 +46,7 @@
 >>> from fipy import input
 >>> x = mesh.cellCenters[0] # doctest: +GMSH
 >>> analyticalArray = valueLeft + (valueRight - valueLeft) * x / Lx # doctest: +GMSH
->>> print var.allclose(analyticalArray, atol=0.025) # doctest: +GMSH
+>>> print(var.allclose(analyticalArray, atol=0.025)) # doctest: +GMSH
 True
 
 >>> errorVar = abs(var - analyticalArray)
@@ -55,7 +55,7 @@
 >>> NonOrthoVar = CellVariable(name="non-orthogonality",
 ...                            mesh=mesh,
 ...                            value=mesh._nonOrthogonality)
->>> print max(NonOrthoVar) < 0.51 # doctest: +GMSH
+>>> print(max(NonOrthoVar) < 0.51) # doctest: +GMSH
 True
 
 >>> if __name__ == '__main__':