1 files changed, 52 insertions, 0 deletions
diff --git a/buch/chapters/040-rekursion/images/gammaplot.m b/buch/chapters/040-rekursion/images/gammaplot.m
new file mode 100644
index 0000000..3b1d23d
--- /dev/null
+++ b/buch/chapters/040-rekursion/images/gammaplot.m
@@ -0,0 +1,52 @@
+#
+# gammaplot.m
+#
+# (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule
+#
+
+fn = fopen("gammapaths.tex", "w");
+
+function finterval(f, fn, from, to, name, delta)
+	fprintf(fn, "\\def\\gamma%s{", name);
+	x = from + delta;
+	fprintf(fn, "({\\dx*%.4f},{\\dy*%.4f})", x, f(x));
+	x = from + 0.02;
+	for x = (from+0.02:0.02:to-0.02)
+		fprintf(fn, "\n -- ");
+		fprintf(fn, "({\\dx*%.4f},{\\dy*%.4f})", x, f(x));
+	endfor
+	x = to - delta;
+	fprintf(fn, "\n -- ");
+	fprintf(fn, "({\\dx*%.4f},{\\dy*%.4f})", x, f(x));
+	fprintf(fn, "}\n");
+endfunction
+
+function gammainterval(fn, from, to, name, delta)
+	finterval(@gamma, fn, from, to, name, delta)
+endfunction
+
+function retval = gammasin(x) 
+	retval = gamma(x) + sin(x * pi);
+endfunction
+
+function gammasininterval(fn, from, to, name, delta)
+	finterval(@gammasin, fn, from, to, name, delta)
+endfunction
+
+gammainterval(fn,  0,  4.1, "plus",  0.019);
+gammainterval(fn, -1,  0,   "one",   0.019);
+gammainterval(fn, -2, -1,   "two",   0.019);
+gammainterval(fn, -3, -2,   "three", 0.019);
+gammainterval(fn, -4, -3,   "four",  0.005);
+gammainterval(fn, -5, -4,   "five",  0.001);
+gammainterval(fn, -6, -5,   "six",  0.0002);
+
+gammasininterval(fn,  0,  4.1, "sinplus",  0.019);
+gammasininterval(fn, -1,  0,   "sinone",   0.019);
+gammasininterval(fn, -2, -1,   "sintwo",   0.019);
+gammasininterval(fn, -3, -2,   "sinthree", 0.019);
+gammasininterval(fn, -4, -3,   "sinfour",  0.005);
+gammasininterval(fn, -5, -4,   "sinfive",  0.001);
+gammasininterval(fn, -6, -5,   "sinsix",  0.0002);
+
+fclose(fn);