1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
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);
|
