OPTION ANGLE DEGREES
LET r = 1.5234
PRINT
PRINT "火星の位置： 黄道 ０ 〜 ３６０°で表現"
PRINT "太陽をはさんで反対方向からスタート"
PRINT "地球の公転半径を １ とする距離の単位"
PRINT
FOR i=0 TO 780 STEP 10
   LET f=0
   LET Q1=i*360/365+180
   LET Q2=i*360/687
   IF r*COS(Q2)-COS(Q1)=0 THEN GOTO 10
   LET a=(r*SIN(Q2)-SIN(Q1))/(r*COS(Q2)-COS(Q1))
   LET x=ATN(a)
   IF r*SIN(Q2)<SIN(Q1) THEN
      LET f=360
   ELSEIF r*SIN(Q2)=SIN(Q1) THEN
      IF r*COS(Q2)>COS(Q1) THEN
         LET x=0
      ELSE
         LET x=180
      END IF
   END IF
   IF r*COS(Q2)<COS(Q1) THEN
      LET f=180
   END IF
10    IF r*COS(Q2)-COS(Q1)=0 THEN
         IF r*SIN(Q2)>SIN(Q1) THEN
            LET x=90
         ELSE
            LET x=270
         END IF
      END IF
      LET xx=ROUND(f+x)
      LET yy=ROUND(SQR((r*COS(Q2)-COS(Q1))^2+(r*SIN(Q2)-SIN(Q1))^2)*100)/100
      PRINT USING "####":i;
      PRINT " 日後： ";xx;"度";"    ( 距離:";yy; ")"
   NEXT i
   PRINT
END