diff options
Diffstat (limited to 'vorlesungen/slides/fresnel/eulerspirale.m')
| -rw-r--r-- | vorlesungen/slides/fresnel/eulerspirale.m | 35 |
1 files changed, 31 insertions, 4 deletions
diff --git a/vorlesungen/slides/fresnel/eulerspirale.m b/vorlesungen/slides/fresnel/eulerspirale.m index 312541a..84e3696 100644 --- a/vorlesungen/slides/fresnel/eulerspirale.m +++ b/vorlesungen/slides/fresnel/eulerspirale.m @@ -4,12 +4,15 @@ # (c) 2022 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschlue # global n; -n = 10000; +n = 1000; global tmax; tmax = 10; +global N; +N = round(n*5/tmax); function retval = f(x, t) - retval = [ cos(t*t); sin(t*t) ]; + x = pi * t^2 / 2; + retval = [ cos(x); sin(x) ]; endfunction x0 = [ 0; 0 ]; @@ -23,12 +26,36 @@ 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"); +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"); +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); |