diff options
author | Andreas Müller <andreas.mueller@ost.ch> | 2022-05-19 17:28:33 +0200 |
---|---|---|
committer | Andreas Müller <andreas.mueller@ost.ch> | 2022-05-19 17:28:33 +0200 |
commit | 32f1a1d818f0fe28b2ae97071e31a773ee2d028a (patch) | |
tree | 4615672911c7c79223bd8c4819c7ed9c31ef6d67 /buch/papers/fresnel/eulerspirale.m | |
parent | typos (diff) | |
download | SeminarSpezielleFunktionen-32f1a1d818f0fe28b2ae97071e31a773ee2d028a.tar.gz SeminarSpezielleFunktionen-32f1a1d818f0fe28b2ae97071e31a773ee2d028a.zip |
some local changes
Diffstat (limited to 'buch/papers/fresnel/eulerspirale.m')
-rw-r--r-- | buch/papers/fresnel/eulerspirale.m | 61 |
1 files changed, 0 insertions, 61 deletions
diff --git a/buch/papers/fresnel/eulerspirale.m b/buch/papers/fresnel/eulerspirale.m deleted file mode 100644 index 84e3696..0000000 --- a/buch/papers/fresnel/eulerspirale.m +++ /dev/null @@ -1,61 +0,0 @@ -# -# 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); |