diff --git a/demo/lbb-stability/Makefile b/demo/lbb-stability/Makefile
index e764fac..022974e 100644
--- a/demo/lbb-stability/Makefile
+++ b/demo/lbb-stability/Makefile
@@ -3,6 +3,23 @@ results-2d.pdf: results-2d.json make_plots.py
 	python make_plots.py --input $< --output $@
 
 results-2d.json: lbb-stability.py
-	python lbb-stability.py --dimension 2 --element mini --output $@
-	python lbb-stability.py --dimension 2 --element taylor-hood --output $@
-	python lbb-stability.py --dimension 2 --element crouzeix-raviart --output $@
+	for element in mini taylor-hood crouzeix-raviart; do \
+		python lbb-stability.py --dimension 2 --element $$element --output $@ ; \
+	done
+
+results-3d.json: lbb-stability.py
+	@for element in mini taylor-hood crouzeix-raviart; do		\
+		for num-layers in 1 2 3; do 				\
+			for vdegree in 1 2 3 4; do 			\
+				python lbb-stability.py 		\
+					--dimension 3 			\
+					--element $$element 		\
+					--num-cells-min 4 		\
+					--num-cells-max 32 		\
+					--num-cells-step 2 		\
+					--num-layers $$num-layers	\
+					--vdegree $$vdegree 		\
+					--output $@ ; 			\
+			done						\
+		done							\
+	done