PRINT
PRINT "＜ある直線で反射して点Ｃから点Ｄへ届く光の軌跡＞"
PRINT
DIM S(1,3),Q(1,2),W(1,2)
MAT INPUT PROMPT "A*x+B*y=K    A, B, K  ":S
PRINT
MAT INPUT PROMPT "点Ｃ（ CX, CY )    CX, CY  ":Q
MAT INPUT PROMPT "点Ｄ（ DX, DY )    DX, DY  ":W
LET Aa=S(1,1)
LET Bb=S(1,2)
LET Kk=S(1,3)
LET CX=Q(1,1)
LET CY=Q(1,2)
LET DX=W(1,1)
LET DY=W(1,2)
IF CX<DX THEN    ! もし CX < DX ならば CX と DX を入れ替える
   LET RX=CX
   LET RY=CY
   LET CX=DX
   LET CY=DY
   LET DX=RX
   LET DY=RY
END IF
IF Aa=0 THEN
   IF (CY-Kk/Bb)*(DY-Kk/Bb)<0 THEN
      PRINT  "点Ｃと点Ｄはこの直線に対してお互い反対側に存在するので光は届きません。"
      GOTO 50
   END IF
   LET DYY=2*Kk/Bb-DY
   LET DXX=DX
   ! (x,y)=(DXX,DYY) は 点Ｄのこの直線を対象軸とする鏡面点Ｄ’
   LET V=1
   IF CX=DXX THEN
      LET LA=-1
      LET V=0
      LET LB=CX
      GOTO 10
   END IF   
   LET LA=(CY-DYY)/(CX-DXX)
   LET LB=CY-CX*LA
   ! y = LA * x + LB  は ２点を結ぶ直線    
10    CALL RENRITU(Aa,Bb,Kk,-1*LA,V,LB,XXX,YYY)   ! ２直線の交点を求める
      GOTO 40
   END IF    
   IF Bb=0 THEN
      IF (CX-Kk/Aa)*(DX-Kk/Aa)<0 THEN
         PRINT  "点Ｃと点Ｄはこの直線に対してお互い反対側に存在するので光は届きません。"
         GOTO 50
      END IF
      LET DXX=2*Kk/Aa-DX
      LET DYY=DY  
      ! (x,y)=(DXX,DYY) は 点Ｄのこの直線を対象軸とする鏡面点Ｄ’
      LET HA=(CY-DYY)/(CX-DXX)
      LET HB=CY-CX*HA
      ! y = HA * x + HB  は ２点を結ぶ直線
      CALL RENRITU(Aa,Bb,Kk,-1*HA,1,HB,XXX,YYY)   ! ２直線の交点を求める
      GOTO 40
   END IF
   LET A=-Aa/Bb
   LET B=K/Bb   
   IF CX-DX=0 THEN
      IF (CY-A*CX-B)*(DY-A*DX-B)<O THEN
         PRINT  "点Ｃと点Ｄはこの直線に対してお互い反対側に存在するので光は届きません。"
         GOTO 50
      END IF
      GOTO 20
   END IF
   LET FA=(CY-DY)/(CX-DX)
   LET FB=CY-CX*FA
   ! y = FA * x + FB  は ２点を結ぶ直線
   CALL RENRITU(Aa,Bb,Kk,-1*FA,1,FB,X,Y)  ! ２直線の交点を求める
   IF DX<X AND X<CX THEN
      PRINT
      PRINT "点Ｃと点Ｄはこの直線に対してお互い反対側に存在するので光は届きません。"
      GOTO 50
   END IF
20 PRINT
   LET EA=-1/A
   LET EB=DY-DX*EA
   ! y = EA * x + EB  は この直線に対して垂直で点Ｄを通る直線
   CALL RENRITU(Aa,Bb,Kk,-1*EA,1,EB,XX,YY)   ! ２直線の交点を求める
   LET DXX=2*XX-DX
   LET DYY=2*YY-DY
   LET Z=1
   ! (x,y)=(DXX,DYY) は 点Ｄのこの直線を対象軸とする鏡面点Ｄ’
   IF CX=DXX THEN
      LET GA=-1
      LET Z=0
      LET GB=DXX
      GOTO 30
   END IF
   LET GA=(CY-DYY)/(CX-DXX)
   LET GB=CY-CX*GA
   ! y = GA * x + GB  は 点Ｃと点Ｄ’を結ぶ直線
30 CALL RENRITU(Aa,Bb,Kk,-1*GA,Z,GB,XXX,YYY)   ! ２直線の交点を求める
40 PRINT
   PRINT "この直線 : "; 
   IF Bb=0 THEN
      PRINT "χ ＝ ";Kk/Aa
   ELSEIF Aa=0 THEN
      PRINT "ｙ ＝ ";Kk/Bb
   ELSE
      PRINT "ｙ ＝ ";-1*Aa/Bb;"χ ＋ (";Kk/Bb;")"
   END IF
   PRINT
   PRINT "  点Ｃ（ χ，ｙ ）＝ ( ";Q(1,1);"，";Q(1,2);") から放たれた光は、"
   PRINT "この直線上の点（ χ，ｙ ）＝ ( ";XXX;",";YYY;") で反射し、"
   PRINT "点Ｄ（ χ，ｙ ）＝ ( ";W(1,1);"，";W(1,2);") に至ります。"
   PRINT
   PRINT "光の軌跡の道のりの長さは ";SQR((CX-DXX)^2+(CY-DYY)^2);"です。"
50 PRINT
END

EXTERNAL SUB RENRITU(a,b,c,d,e,f,x,y)
! aχ +  bｙ ＝ c     dχ +  eｙ ＝ f
   IF b*d=a*e THEN
      LET x=10^10
      LET y=10^10
      GOTO 60
   END IF 
   LET x=-(e*c-b*f)/(b*d-a*e)
   LET y=(d*c-a*f)/(b*d-a*e)
60 END SUB