diff options
author | tim30b <tim.toenz@ost.ch> | 2022-05-18 13:53:24 +0200 |
---|---|---|
committer | tim30b <tim.toenz@ost.ch> | 2022-05-18 13:53:24 +0200 |
commit | 8cabb06907dbe5d18df7a83d30edea9477d4e643 (patch) | |
tree | aaf669409ba40b168e76d435300f74a457f8e5ee /buch/papers/fresnel/eulerspirale.m | |
parent | Intro chapters (diff) | |
parent | typos (diff) | |
download | SeminarSpezielleFunktionen-8cabb06907dbe5d18df7a83d30edea9477d4e643.tar.gz SeminarSpezielleFunktionen-8cabb06907dbe5d18df7a83d30edea9477d4e643.zip |
Merge remote-tracking branch 'mueller/master'
Diffstat (limited to '')
-rw-r--r-- | buch/papers/fresnel/eulerspirale.m | 61 |
1 files changed, 61 insertions, 0 deletions
diff --git a/buch/papers/fresnel/eulerspirale.m b/buch/papers/fresnel/eulerspirale.m new file mode 100644 index 0000000..84e3696 --- /dev/null +++ b/buch/papers/fresnel/eulerspirale.m @@ -0,0 +1,61 @@ +# +# eulerspirale.m +# +# (c) 2022 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschlue +# +global n; +n = 1000; +global tmax; +tmax = 10; +global N; +N = round(n*5/tmax); + +function retval = f(x, t) + x = pi * t^2 / 2; + retval = [ cos(x); sin(x) ]; +endfunction + +x0 = [ 0; 0 ]; +t = tmax * (0:n) / n; + +c = lsode(@f, x0, t); + +fn = fopen("eulerpath.tex", "w"); + +fprintf(fn, "\\def\\fresnela{ (0,0)"); +for i = (2:n) + fprintf(fn, "\n\t-- (%.4f,%.4f)", c(i,1), c(i,2)); +end +fprintf(fn, "\n}\n\n"); + +fprintf(fn, "\\def\\fresnelb{ (0,0)"); +for i = (2:n) + fprintf(fn, "\n\t-- (%.4f,%.4f)", -c(i,1), -c(i,2)); +end +fprintf(fn, "\n}\n\n"); + +fprintf(fn, "\\def\\Cplotright{ (0,0)"); +for i = (2:N) + fprintf(fn, "\n\t-- ({%.4f*\\dx},{%.4f*\\dy})", t(i), c(i,1)); +end +fprintf(fn, "\n}\n\n"); + +fprintf(fn, "\\def\\Cplotleft{ (0,0)"); +for i = (2:N) + fprintf(fn, "\n\t-- ({%.4f*\\dx},{%.4f*\\dy})", -t(i), -c(i,1)); +end +fprintf(fn, "\n}\n\n"); + +fprintf(fn, "\\def\\Splotright{ (0,0)"); +for i = (2:N) + fprintf(fn, "\n\t-- ({%.4f*\\dx},{%.4f*\\dy})", t(i), c(i,2)); +end +fprintf(fn, "\n}\n\n"); + +fprintf(fn, "\\def\\Splotleft{ (0,0)"); +for i = (2:N) + fprintf(fn, "\n\t-- ({%.4f*\\dx},{%.4f*\\dy})", -t(i), -c(i,2)); +end +fprintf(fn, "\n}\n\n"); + +fclose(fn); |