diff --git a/demo/lbb-stability/Makefile b/demo/lbb-stability/Makefile
index 022974e..9a78350 100644
--- a/demo/lbb-stability/Makefile
+++ b/demo/lbb-stability/Makefile
@@ -9,7 +9,7 @@ results-2d.json: lbb-stability.py
 
 results-3d.json: lbb-stability.py
 	@for element in mini taylor-hood crouzeix-raviart; do		\
-		for num-layers in 1 2 3; do 				\
+		for num_layers in 1 2 3; do 				\
 			for vdegree in 1 2 3 4; do 			\
 				python lbb-stability.py 		\
 					--dimension 3 			\
@@ -17,7 +17,7 @@ results-3d.json: lbb-stability.py
 					--num-cells-min 4 		\
 					--num-cells-max 32 		\
 					--num-cells-step 2 		\
-					--num-layers $$num-layers	\
+					--num-layers $$num_layers	\
 					--vdegree $$vdegree 		\
 					--output $@ ; 			\
 			done						\
diff --git a/demo/lbb-stability/lbb-stability.py b/demo/lbb-stability/lbb-stability.py
index d5816d2..af316ff 100644
--- a/demo/lbb-stability/lbb-stability.py
+++ b/demo/lbb-stability/lbb-stability.py
@@ -100,6 +100,9 @@ def run(*, dimension, num_cells, element, num_layers=None, vdegree=None):
 
 stability_constants = []
 nmin, nmax, nstep = args.num_cells_min, args.num_cells_max, args.num_cells_step
+print(f"Element: {args.element}")
+print(f"#layers: {args.num_layers}")
+print(f"vdegree: {args.vdegree}")
 for num_cells in range(nmin, nmax + 1, nstep):
     h, λ = run(
         dimension=args.dimension,