diff --git a/products/+elltool/+doc/+snip/s_chapter06_section05_snippet01.m b/products/+elltool/+doc/+snip/s_chapter06_section05_snippet01.m
index 2e3f2ef2..a257952e 100644
--- a/products/+elltool/+doc/+snip/s_chapter06_section05_snippet01.m
+++ b/products/+elltool/+doc/+snip/s_chapter06_section05_snippet01.m
@@ -27,7 +27,7 @@
     gMat, vEllObj);
 timeVec = [t1, 0];
 % initial directions:
-nDirs = 5;
+nDirs = 3;
 param1Vec = linspace(0, pi, nDirs);
 param2Vec = linspace(0, 2*pi, nDirs);
 [param1Arr, param2Arr, param3Arr] = meshgrid(param1Vec, param1Vec, param2Vec);
@@ -46,10 +46,3 @@
 
 basisMat = [1 0 0 0; 0 0 1 0].'; % orthogonal basis of (x1, x3) subspace
 psObj = brsObj.projection(basisMat); % reach set projection
-
-% plot projection of reach set internal approximation:
-psObj.plotByIa('b');
-
-%plot backward reach set approximation at time t=0
-psCutObj = psObj.cut(0);
-psCutObj.plotByIa('r');
diff --git a/products/+elltool/+doc/+snip/s_chapter06_section05_snippet02.m b/products/+elltool/+doc/+snip/s_chapter06_section05_snippet02.m
new file mode 100644
index 00000000..4f005ba5
--- /dev/null
+++ b/products/+elltool/+doc/+snip/s_chapter06_section05_snippet02.m
@@ -0,0 +1,3 @@
+% plot projection of reach set internal approximation:
+psObj.plotByIa('b');
+
diff --git a/products/+elltool/+doc/+snip/s_chapter06_section05_snippet03.m b/products/+elltool/+doc/+snip/s_chapter06_section05_snippet03.m
new file mode 100644
index 00000000..3630b5c0
--- /dev/null
+++ b/products/+elltool/+doc/+snip/s_chapter06_section05_snippet03.m
@@ -0,0 +1,3 @@
+%plot backward reach set approximation at time t=0
+psCutObj = psObj.cut(0);
+psCutObj.plotByIa('r');