home *** CD-ROM | disk | FTP | other *** search
/ 64'er 1987 December / 64er_Magazin_87-12_1987_Markt__Technik_de_Side_A.d64 / feigenbaum (.txt) < prev    next >
Encoding:
Commodore BASIC  |  1987-01-01  |  1.1 KB  |  43 lines
  1. 10 rem ****************************
  2. 20 rem *   feigenbaum-diagramm    *
  3. 30 rem * 1987 by stefan vilsmeier *
  4. 40 rem ****************************
  5. 50 :
  6. 60 if a=0 then a=1:load "fractal.obj",8,1
  7. 70 open 1,8,15,"u9":close 1: rem diese zeile kann bei problemen mit dem
  8. 80 rem floppyspeeder weggelassen werden!
  9. 90 poke53280,0:poke53281,0:dimh%(128):print"[147]"
  10. 100 input "anfangswert ";x0:x=.5
  11. 110 rem ************ grafik ein *****
  12. 120 sys 50176,14,0:sys 50179,1:sys 50194
  13. 130 :
  14. 140 :
  15. 150 rem ***************************
  16. 160 rem *     hauptschleife       *
  17. 170 rem ***************************
  18. 180 :
  19. 190 for i=0 to 319
  20. 200 :w=(4-x0)/319:a=x0+w*i
  21. 210 :for t=0 to 100
  22. 220 : x=a*x*(1-x)
  23. 230 :next t
  24. 240 :for t=0 to 100
  25. 250 : x=a*x*(1-x)
  26. 260 : sys50182,i,200-200*x,1
  27. 270 :next t
  28. 280 next i
  29. 300 sys 50185,0,0,319,0,1:sys 50185,319,0,319,199,1
  30. 310 sys 50185,319,199,0,199,1:sys 50185,0,199,0,0,1
  31. 320 get a$:if a$="" goto 320
  32. 330 sys 50179,0:if a$="s" then gosub 400: rem 'grafik speichern'
  33. 340 goto 100: rem neustart
  34. 350 :
  35. 360 :
  36. 370 rem ***************************
  37. 380 rem *   'grafik speichern'    *
  38. 390 rem ***************************
  39. 400 :
  40. 410 input "grafik-name ";n$
  41. 420 open 2,8,2,"pi."+n$+",p,w":sys 50191:close 2
  42. 430 return
  43.