DIM r(4)
PRINT
PRINT "人工衛星の高度から周期を導きます。"
PRINT "     （ 静止衛星の高度は 35786 km )"
PRINT "     （ 地球の半径は 6378 km ）"
PRINT
LET r(1)=6378000
LET r(2)=6378000+(42164000-6378000)/2
LET r(3)=42164000
LET r(4)=6378000+(42164000-6378000)*2
FOR i=1 TO 4
   LET w=SQR((6.6743*10^(-11))*(5.9722*10^24)/r(i)^3)    ! 角速度
   LET t=ROUND(2*PI/w,0)                                 ! 周期
   CALL TIME(t,X,Y,Z)
   PRINT "   高度：";
   PRINT USING "#######":(r(i)-6378000)/10^3;
   PRINT " km   → ";
   PRINT USING "####":X;
   PRINT " 時間";
   PRINT USING "####":Y;
   PRINT " 分";
   PRINT USING "####":Z;
   PRINT " 秒"
NEXT i
PRINT
PRINT "  ちなみに、周期の２乗は、地球の中心からの距離の３乗に"
PRINT "比例します。"
PRINT "  また、角速度の２乗は、地球の中心からの距離の３乗に"
PRINT "反比例します。"
PRINT "  また、速度の２乗は、地球の中心からの距離に"
PRINT "反比例します。"
PRINT
END

EXTERNAL SUB TIME(q,A,B,C)
LET A=INT(q/(60*60))
LET B0=MOD(q,60*60)
LET B=INT(B0/60)
LET C=MOD(B0,60)
END SUB