問題9 例16
--------------------------------------
LET x=12345678
LET x_max=x
LET x_min=x
PRINT 1,x
FOR k=2 TO 10453
   LET y=x^2
   LET z=INT(y/10000)
   LET x=MOD(z,100000000)
   PRINT k,x
NEXT k
END
---------------------------------------

問題9 例16(擬似一様乱数)
---------------------------------------
OPTION BASE 0
DIM x(1000),y(1000),a(19)
LET x(1)=INT(TIME)
LET x_max=x(1)
LET x_min=x(1)
PRINT 1,x(1)
FOR k=1 TO 999
   LET x(k+1)=MOD(int(x(k)^2/10000),100000000)
   PRINT k+1,x(k+1)
   IF x(k+1)>x_max THEN LET x_max=x(k+1)
   IF x(k+1)<x_min THEN LET x_min=x(k+1)
NEXT k
PRINT "max=";x_max
PRINT "min=";x_min
FOR k=1 TO 1000
   LET y(k)=(x(k)-x_min)/(x_max-x_min+1)
   LET z=INT(y(k)*20)
   LET a(z)=a(z)+1
NEXT k

SET WINDOW -1,21,-1,100
DRAW axes(1,10)
FOR z=0 TO 19
   PRINT USING "%.##〜%.##   ####":0.05*z,0.05*(z+1),a(z)
   PLOT LINES:z,0;z,a(z);z+1,a(z);z+1,0
NEXT z
END
---------------------------------------------

問題9 例17
--------------------------------------
LET a=2341
LET c=7
LET x=1234
LET x_max=x
LET x_min=x
PRINT 1,x
FOR k=1 TO 10000
   LET x=MOD(a*x+c,10000)
   PRINT k+1,x
   IF x>x_max THEN LET x_max=x
   IF x<x_min THEN LET x_min=x
NEXT k
PRINT "max=";x_max
PRINT "min=";x_min
END
---------------------------------------

問題9 例17(擬似一様乱数)
---------------------------------------
OPTION BASE 0
DIM x(1000),y(1000),b(19)
LET a=2341
LET c=7
LET x(1)=MOD(INT(TIME),10000)
LET x_max=x(1)
LET x_min=x(1)
PRINT 1,x(1)
FOR k=1 TO 999
   LET x(k+1)=MOD(a*x(k)+c,10000)
   PRINT k+1,x(k+1)
   IF x(k+1)>x_max THEN LET x_max=x(k+1)
   IF x(k+1)<x_min THEN LET x_min=x(k+1)
NEXT k
PRINT "max=";x_max
PRINT "min=";x_min
FOR k=1 TO 1000
   LET y(k)=(x(k)-x_min)/(x_max-x_min+1)
   LET z=INT(y(k)*20)
   LET b(z)=b(z)+1
NEXT k

SET WINDOW -1,21,-10,100
DRAW axes(1,10)
FOR z=0 TO 19
   PRINT USING "%.##〜%.##   ####":0.05*z,0.05*(z+1),b(z)
   SET AREA COLOR z+1
   PLOT AREA:z,0;z,b(z);z+1,b(z);z+1,0
   PLOT LINES:z,0;z,b(z);z+1,b(z);z+1,0
NEXT z
END
---------------------------------------------


問題10 例16
----------------------------------------------
OPTION BASE 0
DIM x(1000),y(1000),a(19)
LET x(1)=INT(TIME)
LET x_max=x(1)
LET x_min=x(1)
PRINT 1,x(1)
FOR k=1 TO 999
   LET x(k+1)=MOD(int(x(k)^2/10000),100000000)
   PRINT k+1,x(k+1)
   IF x(k+1)>x_max THEN LET x_max=x(k+1)
   IF x(k+1)<x_min THEN LET x_min=x(k+1)
NEXT k
PRINT "max=";x_max
PRINT "min=";x_min
FOR k=1 TO 1000
   LET y(k)=(x(k)-x_min)/(x_max-x_min+1)
   LET z=INT(y(k)*6)+1
   LET a(z)=a(z)+1
NEXT k

SET WINDOW -1,7,-1,200
DRAW axes(1,20)
FOR z=1 TO 6
   PRINT z,a(z)
   PLOT LINES:z-1,0;z-1,a(z);z,a(z);z,0
NEXT z
END
--------------------------------------------------


問題10 例17
------------------------------------
OPTION BASE 0
DIM x(1000),y(1000),b(19)
LET a=2341
LET c=7
LET x(1)=MOD(INT(TIME),10000)
LET x_max=x(1)
LET x_min=x(1)
PRINT 1,x(1)
FOR k=1 TO 999
   LET x(k+1)=MOD(a*x(k)+c,10000)
   PRINT k+1,x(k+1)
   IF x(k+1)>x_max THEN LET x_max=x(k+1)
   IF x(k+1)<x_min THEN LET x_min=x(k+1)
NEXT k
PRINT "max=";x_max
PRINT "min=";x_min
FOR k=1 TO 1000
   LET y(k)=(x(k)-x_min)/(x_max-x_min+1)
   LET z=INT(y(k)*6)+1
   LET b(z)=b(z)+1
NEXT k

SET WINDOW -1,7,-10,200
DRAW axes(1,20)
FOR z=1 TO 6
   PRINT z,b(z)
   SET AREA COLOR z+1
   PLOT AREA:z-1,0;z-1,b(z);z,b(z);z,0
   PLOT LINES:z-1,0;z-1,b(z);z,b(z);z,0
NEXT z
END
--------------------------------------------------------