home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
64'er 1988 February
/
64er_Magazin_88-02_1988_Markt__Technik_de_Side_A.d64
/
newton-iteration
(
.txt
)
< prev
next >
Wrap
Commodore BASIC
|
2022-10-26
|
547b
|
24 lines
10 rem ****************************
20 rem * newton-iteration *
30 rem * 1987 by stefan vilsmeier *
40 rem ****************************
50 :
60 poke53280,0:poke53281,0
70 input "[147]startwert x0 ";x0
80 xm=x0:print
90 :
100 rem ****** funktion ************
110 :
120 def fnd(x)=x^2-1
130 def fne(x)=2*x :rem ableitung
140 :
150 rem ***************************
160 rem * hauptschleife *
170 rem ***************************
180 :
190 if fne(xm)=0 then print "geht nicht!":end
200 x=x-fnd(xm)/fne(xm)
210 i=i+1:print i,x
220 if abs(x-xm)>.0005 then xm=x:goto 190
230 print "fertig!"