home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / basic / library / trubasic / chaos / xhenon.tru < prev    next >
Encoding:
Text File  |  1980-01-01  |  1.1 KB  |  45 lines
  1. !HENON MAP
  2. LIBRARY "SGLIB.TRC"
  3. CLEAR
  4. PRINT"                               ***HENON MAP***"
  5. PRINT
  6. PRINT"THIS PROGRAM GENERATES THE (X,Y) PHASE DIAGRAM FOR THE HENON MAP."
  7. PRINT "        X(N+1)=1-A*X(N)^2 + Y(N)"
  8. PRINT "        Y(N+1)=B*X(N)"
  9. PRINT" TWO PARAMETERS ARE REQUIRED, A AND B. IF B=1 THE MAP IS "
  10. PRINT"CONSERVATIVE.  IF B < ABS(1) THEN THE MAP IS DISSIPATIVE."
  11. PRINT"TRY A=1.4 AND B=0.3 INITIALLY."
  12. PRINT
  13. !
  14. DIM L(1),M(1)
  15. INPUT PROMPT"INPUT INITIAL X , Y VALUES:":XINT,YINT
  16. INPUT PROMPT"INPUT A,B VALUES:":A,B
  17. INPUT PROMPT"INPUT NUMBER OF THROWAWAY ITERATIONS:":INITNUM
  18. INPUT PROMPT"INPUT NUMBER OF TOTAL ITERATIONS:":FINNUM
  19. !
  20. CALL SETXSCALE(-1.5,1.5)
  21. CALL SETYSCALE(-.5,.5)
  22. CALL SETTEXT("HENON MAP","X","Y")
  23. CALL RESERVELEGEND
  24. !
  25. DATA 0,0
  26. CALL DATAGRAPH(L,M,1,0,"WHITE")
  27. CALL GOTOCANVAS
  28. LET X=XINT
  29. LET Y=YINT
  30. FOR I=INITNUM+1 TO FINNUM+1
  31.     LET XNEW=1-A*(X^2)+Y
  32.     LET YNEW=B*X
  33.     CALL GRAPHPOINT(XNEW,YNEW,1)
  34.     LET X=XNEW
  35.     LET Y=YNEW
  36. NEXT I
  37. CALL ADDLEGEND("A="&STR$(A)&"   B="&STR$(B),0,1,"WHITE")
  38. CALL DRAWLEGEND
  39. get key variable
  40. clear
  41. print"press <esc> key to finish"
  42. END
  43.  
  44.  
  45.