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 >
Encoding:
Commodore BASIC  |  1988-01-01  |  547 b   |  24 lines
  1. 10 rem ****************************
  2. 20 rem *   newton-iteration       *
  3. 30 rem * 1987 by stefan vilsmeier *
  4. 40 rem ****************************
  5. 50 :
  6. 60 poke53280,0:poke53281,0
  7. 70 input "[147]startwert x0    ";x0
  8. 80 xm=x0:print
  9. 90 :
  10. 100 rem ****** funktion ************
  11. 110 :
  12. 120 def fnd(x)=x^2-1
  13. 130 def fne(x)=2*x     :rem ableitung
  14. 140 :
  15. 150 rem ***************************
  16. 160 rem *   hauptschleife         *
  17. 170 rem ***************************
  18. 180 :
  19. 190 if fne(xm)=0 then print "geht nicht!":end
  20. 200 x=x-fnd(xm)/fne(xm)
  21. 210 i=i+1:print i,x
  22. 220 if abs(x-xm)>.0005 then xm=x:goto 190
  23. 230 print "fertig!"
  24.