diff --git a/convert_to_video b/convert_to_video
new file mode 100644
index 0000000..99fb7a8
--- /dev/null
+++ b/convert_to_video
@@ -0,0 +1 @@
+ffmpeg -f image2 -framerate 20 -i test%03d.png -s 1024x1024 -c:v libvpx-vp9 -lossless 1 test.webm
diff --git a/max_slope.plt b/max_slope.plt
new file mode 100644
index 0000000..15482d8
--- /dev/null
+++ b/max_slope.plt
@@ -0,0 +1,18 @@
+set log x
+set y2tics
+set xrange [exp(-1):1]
+set yrange [1.5:2]
+set y2range [1.98:2.1]
+set grid
+
+# set terminal pngcairo enhanced size 1500,1000
+# set output "output/max_slope.png"
+
+plot "output/max_slope_1621.dat" using 1:3 w lp pt 7 ps 0.6 lw 2 t "1621 elements", \
+	 "output/max_slope_24428.dat" using 1:3 w lp pt 7 ps 0.6 lw 2 t "24428 elements", \
+	 "output/max_slope_94252.dat" using 1:3 w lp lw 2 pt 7 ps 0.6 t "94252 elements", \
+	 "output/max_slope_713698.dat" using 1:3 w lp lw 2 pt 7 ps 0.6 t "713698 elements", \
+	 "output/max_slope_1621.dat" using 1:2 w p pt 7 ax x1y2 t "parameter"
+
+pause mouse keypress
+if(MOUSE_KEY != 113) reread
diff --git a/singular_values_movie.plt b/singular_values_movie.plt
new file mode 100644
index 0000000..d585540
--- /dev/null
+++ b/singular_values_movie.plt
@@ -0,0 +1,23 @@
+if(!exists("i")) i = 0
+
+file = sprintf("< ./singular_values %f 1.78", exp((i-50)*0.02))
+
+set samples 1000
+set size square
+set xrange [0:30]
+set yrange [0:30]
+set trange [0:30]
+set grid
+set parametric
+
+set terminal pngcairo enhanced size 1024,1024
+img = sprintf("output/animation/test%03d.png", i);
+print sprintf("write %s", img)
+set output img
+
+plot file using 6:7 w p pt 7 ps 0.5 lc 1 t columnheader, \
+	 t,2*t w l lw 2 t "", \
+	 t,t/2 w l lw 2 t ""
+
+i=i+1
+if(i <= 100) reread