diff options
author | Andreas Müller <andreas.mueller@ost.ch> | 2021-12-23 22:52:31 +0100 |
---|---|---|
committer | Andreas Müller <andreas.mueller@ost.ch> | 2021-12-23 22:52:31 +0100 |
commit | a97ba371f257ffcb29467142d00c59c4fe7205df (patch) | |
tree | dbcf49f4f56531b33e16c3866ac9bdf9ac916c10 /buch/chapters/030-geometrie/trigo/1.m | |
parent | typo (diff) | |
download | SeminarSpezielleFunktionen-a97ba371f257ffcb29467142d00c59c4fe7205df.tar.gz SeminarSpezielleFunktionen-a97ba371f257ffcb29467142d00c59c4fe7205df.zip |
more trigo stuff
Diffstat (limited to '')
-rw-r--r-- | buch/chapters/030-geometrie/trigo/1.m | 33 |
1 files changed, 33 insertions, 0 deletions
diff --git a/buch/chapters/030-geometrie/trigo/1.m b/buch/chapters/030-geometrie/trigo/1.m new file mode 100644 index 0000000..3eea249 --- /dev/null +++ b/buch/chapters/030-geometrie/trigo/1.m @@ -0,0 +1,33 @@ +# +# 1.m +# +# +n = 10; +format long; + +function retval = snewtonstep(x, s) + corr = (-4*x*x*x+3*x-s) / (-12*x*x + 3); + retval = x - corr; +end + +function retval = cnewtonstep(x, c) + corr = (4*x*x*x-3*x-c) / (12*x*x - 3); + retval = x - corr; +end + +s0 = pi / 180; + +s3 = sind(3) +c3 = cosd(3) + +r = zeros(n+1,2) +r(1,1) = pi / 180; +r(1,2) = sqrt(1-r(1,1)^2); + +for i = (1:n) + r(i+1,1) = snewtonstep(r(i,1), s3); + r(i+1,2) = cnewtonstep(r(i,2), c3); +endfor + +r + |