home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / basic / library / trubasic / chaos / xreadbas.tru < prev    next >
Encoding:
Text File  |  1980-01-01  |  1.5 KB  |  51 lines
  1.  
  2. !Pendulum Program "XREADBAS"INS OF ATTRACTION
  3. LIBRARY "SGLIB.TRC"
  4. DIM A(1),B(1)
  5. CLEAR
  6. PRINT"                        ***BASINS OF ATTRACTION***"
  7. PRINT"THIS PROGRAM READS TWO DATA FILES CONTAINING 'BASIN OF ATTRACTION' DATA"
  8. PRINT"AND THE CORRESPONDING POINCARE SECTIONS FROM EACH SET OF INITIAL  "
  9. PRINT"CONDITIONS.  THE DATA IS PLOTTED ON THE PHASE PLANE."
  10. PRINT
  11. INPUT PROMPT" INPUT G VALUE:":G
  12. INPUT PROMPT"DAMPING FACTOR:":Q
  13. INPUT PROMPT"WHICH DRIVE FOR DATA DISK:":D$
  14. INPUT PROMPT"NAME OF 'BASINS OF ATTRACTION' FILE:":NAME$
  15. INPUT PROMPT"DO YOU WISH CORRESPONDING POINCARE FILE, YES(1), NO(2) :":P
  16. IF P = 1 THEN
  17.    INPUT PROMPT"NAME OF POINCARE FILE:":POINFILE$
  18.    INPUT PROMPT"INPUT TYPE OF SYMBOL FOR POINCARE PLOT, DOT(1),ASTERISK(3),BOX(10):":STYLE
  19. END IF
  20. !
  21. CALL SETTEXT("PENDULUM - BASINS OF ATTRACTION","INIT. ANGLE","INIT ANG VEL")
  22. CALL SETXSCALE(-3,3)
  23. CALL SETYSCALE(-3,3)
  24. CALL RESERVELEGEND
  25. DATA 0,0
  26. CALL DATAGRAPH(A,B,1,0,"WHITE")
  27. CALL GOTOCANVAS
  28. OPEN #1: NAME D$&":"&NAME$, ORGANIZATION RECORD
  29. SET#1: RECORD 1
  30. WHEN ERROR IN
  31.      DO WHILE MORE #1
  32.         READ #1: XINT,VINT
  33.            CALL GRAPHPOINT(XINT,VINT,4)
  34.      LOOP
  35. USE
  36. END WHEN
  37. IF P = 1 THEN
  38.    OPEN #2:NAME D$&":"&POINFILE$, ORGANIZATION RECORD
  39.    SET#1: RECORD 1
  40.    WHEN ERROR IN
  41.         DO WHILE MORE #1
  42.            READ #2: X,V
  43.            CALL GRAPHPOINT(X,V,STYLE)
  44.         LOOP
  45.    USE
  46.    END WHEN
  47. END IF
  48. CALL ADDLEGEND("G="&STR$(G)&"  Q="&STR$(Q),0,1,"WHITE")
  49. CALL DRAWLEGEND
  50. END
  51.