PROGRAM LSYS_KOCH_F_CURVE   !© W.van Duyn 19th May 2004
SET WINDOW 0,1199,0,903
SET background color "white"
SET color 16
PRINT "location: x try 600, y try 450"
INPUT x,y
PRINT"Number of iterations try 4 or 5, size try 100, angle try 4"
INPUT depth,size,ak
CLEAR
SET COLOR "black"
BOX AREA 0,1199,0,903
LET k=360/ak
LET r=2
LET d=size/r^depth
LET a=0
SET COLOR 14
CALL lsys1(depth,d,a,k,x,y)
CALL turnRIGHT(a,k)
CALL turnRIGHT(a,k)
FOR i=1 to ak
SET color i+9
CALL lsys1(depth,d,a,k,x,y)
CALL turnRIGHT(a,k)
NEXT i
END

SUB lsys1(depth,d,a,k,x,y)
IF depth=0 THEN
CALL moveON(d,x,y,a)
EXIT IF
ELSE
CALL turnLEFT(a,k)
CALL lsys2(depth-1,d,a,k,x,y)
CALL push(d,x,y,a)
CALL lsys1(depth-1,d,a,k,x,y)
CALL turnRIGHT(a,k)
CALL lsys1(depth-1,d,a,k,x,y)
CALL pop(d,x,y,a)
CALL turnRIGHT(a,k)
CALL lsys1(depth-1,d,a,k,x,y)
CALL lsys1(depth-1,d,a,k,x,y)
END IF
END SUB

SUB lsys2(depth,d,a,k,x,y)
IF depth=0 THEN
CALL moveON(d,x,y,a)
EXIT IF
ELSE
CALL lsys2(depth-1,d,a,k,x,y)
CALL lsys2(depth-1,d,a,k,x,y)
CALL push(d,x,y,a)
CALL turnRIGHT(a,k)
CALL lsys1(depth-1,d,a,k,x,y)
CALL turnRIGHT(a,k)
CALL lsys1(depth-1,d,a,k,x,y)
CALL pop(d,x,y,a)
CALL turnLEFT(a,k)
CALL lsys1(depth-1,d,a,k,x,y)
CALL turnRIGHT(a,k)
END IF
END SUB

SUB moveON(d,x,y,a)
OPTION ANGLE degrees
LET dx=d*cos(a)
LET dy=d*sin(a)
PLOT x,y;x+dx,y+dy
LET x=x+dx
LET y=y+dy
END SUB

SUB turnRIGHT(a,k)
LET a=a-k
END SUB

SUB turnLEFT(a,k)
LET a=a+k
END SUB

MODULE stack
SHARE state(40,4),stack
LET stack=1

SUB push(d,x,y,a)   !save the graphics state
LET state(stack,1)=d
LET state(stack,2)=x
LET state(stack,3)=y
LET state(stack,4)=a
LET stack=stack+1
END SUB

SUB pop(d,x,y,a)   !restore the graphics state
LET stack=stack-1
LET d=state(stack,1)
LET x=state(stack,2)
LET y=state(stack,3)
LET a=state(stack,4)
END SUB

END MODULE