Re: マンデルブロート - しばっち

2026/07/19 (Sun) 07:51:54

軌道トラップによる描画


FOR I=0 TO 7
SET COLOR MIX(I) BITAND(I,2)/2,BITAND(I,4)/4,BITAND(I,1)
NEXT I
CLEAR
LET XS=-2
LET YS=-2
LET XE=2
LET YE=2
SET WINDOW XS,XE,YE,YS
SET POINT STYLE 1
ASK BITMAP SIZE XSIZE,YSIZE
LET KMAX=50
LET N=2
FOR Y=0 TO YSIZE
FOR X=0 TO XSIZE
LET A=COMPLEX(WORLDX(X),WORLDY(Y))
LET Z=0
LET MINDIST=1E+8
FOR K=1 TO KMAX
LET Z=Z^N+A
LET MINDIST=MIN(MINDIST,ABS(Z))
IF ABS(Z)>2 THEN
LET COL=INT(255*(1-MINDIST))
SET POINT COLOR MAX(0,MIN(COL,255))
PLOT POINTS :WORLDX(X),WORLDY(Y)
EXIT FOR
END IF
NEXT K
NEXT X
NEXT Y
END

Re: マンデルブロート - しばっち

2026/07/19 (Sun) 07:53:09

Zの距離を利用したもの


FOR I=0 TO 7
SET COLOR MIX(I) BITAND(I,2)/2,BITAND(I,4)/4,BITAND(I,1)
NEXT I
CLEAR
LET XS=-2
LET YS=-2
LET XE=2
LET YE=2
SET WINDOW XS,XE,YE,YS
SET POINT STYLE 1
ASK BITMAP SIZE XSIZE,YSIZE
LET KMAX=50
LET N=2
FOR Y=0 TO YSIZE
FOR X=0 TO XSIZE
LET A=COMPLEX(WORLDX(X),WORLDY(Y))
LET Z=0
LET DZ=1
FOR K=1 TO KMAX
LET Z=Z^N+A
LET DZ=N*Z^(N-1)*DZ
IF ABS(Z)>2 THEN
LET R=ABS(Z)
LET DR=ABS(DZ)
IF DR=0 THEN
LET COL=0
ELSE
LET D=R*LOG(R)/DR
LET COL=INT(D*200)
END IF
SET POINT COLOR MAX(0,MIN(COL,255))
PLOT POINTS :WORLDX(X),WORLDY(Y)
EXIT FOR
END IF
NEXT K
NEXT X
NEXT Y
END

Re: マンデルブロート - しばっち

2026/07/19 (Sun) 07:54:37

マンデルブロートの着色法には繰り返し数以外にもある。
下記プログラムはZの偏角を利用している。


SET COLOR MIX(0) 0,0,0
SET COLOR MODE "NATIVE"
CLEAR
LET XS=-2
LET YS=-2
LET XE=2
LET YE=2
SET WINDOW XS,XE,YE,YS
SET POINT STYLE 1
ASK BITMAP SIZE XSIZE,YSIZE
LET KMAX=50
LET N=2
FOR Y=0 TO YSIZE
FOR X=0 TO XSIZE
LET A=COMPLEX(WORLDX(X),WORLDY(Y))
LET Z=0
FOR K=1 TO KMAX
LET Z=Z^N+A
IF ABS(Z)>2 THEN
LET H=MOD(DEG(ANGLE(RE(Z),IM(Z)))+360,360)
CALL HSV2RGB(RR,GG,BB,H,255,255)
SET COLOR COLORINDEX(RR/255,GG/255,BB/255)
PLOT POINTS :WORLDX(X),WORLDY(Y)
EXIT FOR
END IF
NEXT K
NEXT X
NEXT Y
END

EXTERNAL SUB HSV2RGB(R,G,B,H,S,V)
IF S=0 THEN
LET R=V
LET G=V
LET B=V
EXIT SUB
END IF
LET T=V-S*V/255
LET HH=H
IF H>=300 OR H<60 THEN
IF H>=300 THEN LET HH=360-HH
IF H<60 THEN LET HH=-HH
LET HH=HH/60
LET RR=0
IF HH<0 THEN
LET BB=1
LET GG=HH+BB
ELSE
LET GG=1
LET BB=GG-HH
END IF
ELSEIF H>=60 AND H<180 THEN
LET HH=HH-120
LET HH=HH/60
LET GG=0
IF HH<0 THEN
LET BB=1
LET RR=HH+BB
ELSE
LET RR=1
LET BB=RR-HH
END IF
ELSEIF H>=180 AND H<300 THEN
LET HH=HH-240
LET HH=HH/60
LET BB=0
IF HH<0 THEN
LET RR=1
LET GG=HH+RR
ELSE
LET GG=1
LET RR=GG-HH
END IF
END IF
LET R=-RR*(V-T)+V
LET R=INT(R)
LET G=-GG*(V-T)+V
LET G=INT(G)
LET B=-BB*(V-T)+V
LET B=INT(B)
END SUB

Re: マンデルブロート - knoike

2026/07/20 (Mon) 06:21:54

本題からは離れた余談なのですが、「Mandelbrot」の t は発音せず、「マンデルブロ（マンデルブロー）」と読みます。
「ballet（バレエ）」とか、
「depot（デポ）」とか、
「buffet（ビュッフェ）」などと同じです。

Re: 4次元グラフィックス - しばっち

2026/07/05 (Sun) 07:47:08

八胞体を表示します。


https://ja.wikipedia.org/wiki/正八胞体


OPTION ANGLE DEGREES
LET LMIN=1E+10
LET LMAX=-1E+10
DIM M(4,4),POINT(4),ROTXY(4,4),ROTXZ(4,4),ROTYW(4,4),ROTYZ(4,4),ROTXW(4,4),ROTZW(4,4)
DIM X(50),Y(50),Z(50),W(50),MESH(50,10)
FOR I=1 TO 16
READ X(I),Y(I),Z(I),W(I)
LET XMAX=MAX(X(I),XMAX)
LET XMIN=MIN(X(I),XMIN)
LET YMAX=MAX(Y(I),YMAX)
LET YMIN=MIN(Y(I),YMIN)
LET ZMAX=MAX(Z(I),ZMAX)
LET ZMIN=MIN(Z(I),ZMIN)
LET WMAX=MAX(W(I),WMAX)
LET WMIN=MIN(W(I),WMIN)
NEXT I
DATA -0.5,-0.5,-0.5,-0.5
DATA -0.5,-0.5,-0.5,0.5
DATA -0.5,-0.5,0.5,-0.5
DATA -0.5,-0.5,0.5,0.5
DATA -0.5,0.5,-0.5,-0.5
DATA -0.5,0.5,-0.5,0.5
DATA -0.5,0.5,0.5,-0.5
DATA -0.5,0.5,0.5,0.5
DATA 0.5,-0.5,-0.5,-0.5
DATA 0.5,-0.5,-0.5,0.5
DATA 0.5,-0.5,0.5,-0.5
DATA 0.5,-0.5,0.5,0.5
DATA 0.5,0.5,-0.5,-0.5
DATA 0.5,0.5,-0.5,0.5
DATA 0.5,0.5,0.5,-0.5
DATA 0.5,0.5,0.5,0.5

DO
READ IF MISSING THEN EXIT DO:N
LET NUM=NUM+1
FOR I=1 TO N
READ MESH(NUM,I)
NEXT I
LOOP
DATA 4,0,2,3,1
DATA 4,0,4,5,1
DATA 4,0,4,6,2
DATA 4,0,8,9,1
DATA 4,0,8,10,2
DATA 4,0,8,12,4
DATA 4,4,6,7,5
DATA 4,2,6,7,3
DATA 4,1,5,7,3
DATA 4,2,10,11,3
DATA 4,1,9,11,3
DATA 4,1,9,13,5
DATA 4,8,10,11,9
DATA 4,8,12,13,9
DATA 4,8,12,14,10
DATA 4,4,12,13,5
DATA 4,4,12,14,6
DATA 4,2,10,14,6
DATA 4,12,14,15,13
DATA 4,10,14,15,11
DATA 4,9,13,15,11
DATA 4,6,14,15,7
DATA 4,5,13,15,7
DATA 4,3,11,15,7
LET RANGE=MAX(XMAX-XMIN,MAX(YMAX-YMIN,MAX(ZMAX-ZMIN,WMAX-WMIN)))
LOCATE VALUE NOWAIT(1),RANGE 0 TO 3,AT 1:SCALE
LOCATE VALUE NOWAIT(2),RANGE -RANGE TO RANGE,AT 0:XMOVE
LOCATE VALUE NOWAIT(3),RANGE -RANGE TO RANGE,AT 0:YMOVE
LOCATE VALUE NOWAIT(4),RANGE -RANGE TO RANGE,AT 0:ZMOVE
LOCATE VALUE NOWAIT(5),RANGE -RANGE TO RANGE,AT 0:WMOVE
LOCATE VALUE NOWAIT(6),RANGE -360 TO 360,AT 0:XY
LOCATE VALUE NOWAIT(7),RANGE -360 TO 360,AT 0:XZ
LOCATE VALUE NOWAIT(8),RANGE -360 TO 360,AT 0:YW
LOCATE VALUE NOWAIT(9),RANGE -360 TO 360,AT 0:YZ
LOCATE VALUE NOWAIT(10),RANGE -360 TO 360,AT 0:XW
LOCATE VALUE NOWAIT(11),RANGE -360 TO 360,AT 0:ZW
DO
LOCATE VALUE NOWAIT(1):SCALE
LOCATE VALUE NOWAIT(2):XMOVE
LOCATE VALUE NOWAIT(3):YMOVE
LOCATE VALUE NOWAIT(4):ZMOVE
LOCATE VALUE NOWAIT(5):WMOVE
LOCATE VALUE NOWAIT(6):XY
LOCATE VALUE NOWAIT(7):XZ
LOCATE VALUE NOWAIT(8):YW
LOCATE VALUE NOWAIT(9):YZ
LOCATE VALUE NOWAIT(10):XW
LOCATE VALUE NOWAIT(11):ZW
MAT ROTXY=IDN !'XY平面上の回転(ZW平面を軸とする回転)
LET ROTXY(1,1)=COS(XY)
LET ROTXY(1,2)=SIN(XY)
LET ROTXY(2,1)=-SIN(XY)
LET ROTXY(2,2)=COS(XY)
MAT ROTXZ=IDN !'XZ平面上の回転(YW平面を軸とする回転)
LET ROTXZ(1,1)=COS(XZ)
LET ROTXZ(1,3)=-SIN(XZ)
LET ROTXZ(3,1)=SIN(XZ)
LET ROTXZ(3,3)=COS(XZ)
MAT ROTYW=IDN !'YW平面上の回転(XZ平面を軸とする回転)
LET ROTYW(2,2)=COS(YW)
LET ROTYW(2,4)=-SIN(YW)
LET ROTYW(4,2)=SIN(YW)
LET ROTYW(4,4)=COS(YW)
MAT ROTYZ=IDN !'YZ平面上の回転(XW平面を軸とする回転)
LET ROTYZ(2,2)=COS(YZ)
LET ROTYZ(2,3)=SIN(YZ)
LET ROTYZ(3,2)=-SIN(YZ)
LET ROTYZ(3,3)=COS(YZ)
MAT ROTXW=IDN !'XW平面上の回転(YZ平面を軸とする回転)
LET ROTXW(1,1)=COS(XW)
LET ROTXW(1,4)=SIN(XW)
LET ROTXW(4,1)=-SIN(XW)
LET ROTXW(4,4)=COS(XW)
MAT ROTZW=IDN !'ZW平面上の回転(XY平面を軸とする回転)
LET ROTZW(3,3)=COS(ZW)
LET ROTZW(3,4)=-SIN(ZW)
LET ROTZW(4,3)=SIN(ZW)
LET ROTZW(4,4)=COS(ZW)
MAT M=ROTXY*ROTXZ*ROTYW*ROTYZ*ROTXW*ROTZW ! 4次元回転
SET DRAW MODE HIDDEN
CLEAR
FOR I=1 TO NUM
CALL PLOT(X(MESH(I,1)+1),Y(MESH(I,1)+1),Z(MESH(I,1)+1),W(MESH(I,1)+1))
CALL PLOT(X(MESH(I,2)+1),Y(MESH(I,2)+1),Z(MESH(I,2)+1),W(MESH(I,2)+1))
CALL PLOT(X(MESH(I,3)+1),Y(MESH(I,3)+1),Z(MESH(I,3)+1),W(MESH(I,3)+1))
CALL PLOT(X(MESH(I,4)+1),Y(MESH(I,4)+1),Z(MESH(I,4)+1),W(MESH(I,4)+1))
CALL PLOT(X(MESH(I,1)+1),Y(MESH(I,1)+1),Z(MESH(I,1)+1),W(MESH(I,1)+1))
PLOT LINES
NEXT I
IF FL=0 THEN
SET WINDOW -LMAX*1.2,LMAX*1.2,-LMAX*1.2,LMAX*1.2
LET FL=1
END IF
SET DRAW MODE EXPLICIT
LOOP

SUB PLOT(X,Y,Z,W)
LET LW=2
LET POINT(1)=X-XMOVE
LET POINT(2)=Y-YMOVE
LET POINT(3)=Z-ZMOVE
LET POINT(4)=W-WMOVE
MAT POINT=POINT*M
! FOR K=1 TO 3
! LET POINT(K)=POINT(K)/(LW-POINT(4))
! NEXT
IF FL=0 THEN
LET LMIN=MIN(LMIN,POINT(1))
LET LMAX=MAX(LMAX,POINT(1))
LET LMIN=MIN(LMIN,POINT(2))
LET LMAX=MAX(LMAX,POINT(2))
ELSE
PLOT LINES:POINT(1)*SCALE,POINT(2)*SCALE;
END IF
END SUB
END

Re: 4次元グラフィックス - しばっち

2026/07/05 (Sun) 07:48:22

十六胞体

https://ja.wikipedia.org/wiki/正十六胞体

以下データのみ



DATA -1,0,0,0
DATA 0,-1,0,0
DATA 0,0,-1,0
DATA 0,0,0,-1
DATA 1,0,0,0
DATA 0,1,0,0
DATA 0,0,1,0
DATA 0,0,0,1



DATA 3,5,6,7
DATA 3,4,6,7
DATA 3,4,5,7
DATA 3,4,5,6
DATA 3,0,6,7
DATA 3,0,5,7
DATA 3,0,5,6
DATA 3,1,6,7
DATA 3,1,4,7
DATA 3,1,4,6
DATA 3,0,1,7
DATA 3,0,1,6
DATA 3,2,5,7
DATA 3,2,4,7
DATA 3,2,4,5
DATA 3,0,2,7
DATA 3,0,2,5
DATA 3,1,2,7
DATA 3,1,2,4
DATA 3,0,1,2
DATA 3,3,5,6
DATA 3,3,4,6
DATA 3,3,4,5
DATA 3,0,3,6
DATA 3,0,3,5
DATA 3,1,3,6
DATA 3,1,3,4
DATA 3,0,1,3
DATA 3,2,3,5
DATA 3,2,3,4
DATA 3,0,2,3
DATA 3,1,2,3

Re: 4次元グラフィックス - しばっち

2026/07/05 (Sun) 07:49:44

二十四胞体


https://ja.wikipedia.org/wiki/正二十四胞体

以下データのみ


DATA -1,0,0,0
DATA 0,-1,0,0
DATA 0,0,-1,0
DATA 0,0,0,-1
DATA 1,0,0,0
DATA 0,1,0,0
DATA 0,0,1,0
DATA 0,0,0,1
DATA -0.5,-0.5,-0.5,-0.5
DATA -0.5,-0.5,-0.5,0.5
DATA -0.5,-0.5,0.5,-0.5
DATA -0.5,-0.5,0.5,0.5
DATA -0.5,0.5,-0.5,-0.5
DATA -0.5,0.5,-0.5,0.5
DATA -0.5,0.5,0.5,-0.5
DATA -0.5,0.5,0.5,0.5
DATA 0.5,-0.5,-0.5,-0.5
DATA 0.5,-0.5,-0.5,0.5
DATA 0.5,-0.5,0.5,-0.5
DATA 0.5,-0.5,0.5,0.5
DATA 0.5,0.5,-0.5,-0.5
DATA 0.5,0.5,-0.5,0.5
DATA 0.5,0.5,0.5,-0.5
DATA 0.5,0.5,0.5,0.5


DATA 3,0,8,9
DATA 3,0,8,10
DATA 3,1,8,9
DATA 3,1,8,10
DATA 3,0,10,11
DATA 3,0,9,11
DATA 3,1,10,11
DATA 3,1,9,11
DATA 3,0,8,12
DATA 3,2,8,9
DATA 3,2,8,12
DATA 3,0,12,13
DATA 3,0,9,13
DATA 3,2,9,13
DATA 3,2,12,13
DATA 3,3,8,10
DATA 3,3,8,12
DATA 3,0,12,14
DATA 3,0,10,14
DATA 3,3,12,14
DATA 3,3,10,14
DATA 3,1,8,16
DATA 3,2,8,16
DATA 3,1,9,17
DATA 3,2,9,17
DATA 3,1,16,17
DATA 3,2,16,17
DATA 3,3,8,16
DATA 3,1,10,18
DATA 3,1,16,18
DATA 3,3,16,18
DATA 3,3,10,18
DATA 3,2,16,20
DATA 3,3,16,20
DATA 3,2,12,20
DATA 3,3,12,20
DATA 3,0,14,15
DATA 3,0,13,15
DATA 3,5,12,13
DATA 3,5,12,14
DATA 3,5,14,15
DATA 3,5,13,15
DATA 3,0,11,15
DATA 3,6,10,11
DATA 3,6,10,14
DATA 3,6,14,15
DATA 3,6,11,15
DATA 3,7,9,11
DATA 3,7,9,13
DATA 3,7,13,15
DATA 3,7,11,15
DATA 3,1,11,19
DATA 3,1,18,19
DATA 3,6,10,18
DATA 3,6,18,19
DATA 3,6,11,19
DATA 3,7,9,17
DATA 3,1,17,19
DATA 3,7,17,19
DATA 3,7,11,19
DATA 3,2,17,21
DATA 3,2,13,21
DATA 3,7,17,21
DATA 3,7,13,21
DATA 3,4,16,17
DATA 3,4,16,18
DATA 3,4,18,19
DATA 3,4,17,19
DATA 3,4,16,20
DATA 3,2,20,21
DATA 3,4,20,21
DATA 3,4,17,21
DATA 3,3,20,22
DATA 3,3,18,22
DATA 3,4,20,22
DATA 3,4,18,22
DATA 3,5,12,20
DATA 3,5,20,21
DATA 3,5,13,21
DATA 3,3,14,22
DATA 3,5,20,22
DATA 3,5,14,22
DATA 3,6,18,22
DATA 3,6,14,22
DATA 3,4,22,23
DATA 3,4,21,23
DATA 3,5,22,23
DATA 3,5,21,23
DATA 3,4,19,23
DATA 3,6,22,23
DATA 3,6,19,23
DATA 3,7,21,23
DATA 3,7,19,23
DATA 3,5,15,23
DATA 3,6,15,23
DATA 3,7,15,23

Re: 4次元グラフィックス - しばっち

2026/07/05 (Sun) 07:51:04

4次元トーラス



OPTION ANGLE DEGREES
LET LMIN=1E+10
LET LMAX=-1E+10
LET NN=40 !'分割数
LET MM=40
DIM M(4,4),POINT(4),ROTXY(4,4),ROTXZ(4,4),ROTYU(4,4),ROTYZ(4,4),ROTXU(4,4),ROTZU(4,4)
DIM XX(0 TO NN,0 TO MM),YY(0 TO NN,0 TO MM),ZZ(0 TO NN,0 TO MM),WW(0 TO NN,0 TO MM)
SET POINT STYLE 1
FOR J=0 TO MM
FOR I=0 TO NN
LET ALPHA=I*360/NN
LET BETA=J*360/MM
LET XX(I,J)=COS(ALPHA)
LET YY(I,J)=SIN(ALPHA)
LET ZZ(I,J)=COS(BETA)
LET WW(I,J)=SIN(BETA)
LET XMIN=MIN(XMIN,XX(I,J))
LET XMAX=MAX(XMAX,XX(I,J))
LET YMIN=MIN(YMIN,YY(I,J))
LET YMAX=MAX(YMAX,YY(I,J))
LET ZMIN=MIN(ZMIN,ZZ(I,J))
LET ZMAX=MAX(ZMAX,ZZ(I,J))
LET WMIN=MIN(WMIN,WW(I,J))
LET WMAX=MAX(WMAX,WW(I,J))
NEXT I
NEXT J
LET RANGE=MAX(XMAX-XMIN,MAX(YMAX-YMIN,MAX(ZMAX-ZMIN,WMAX-WMIN)))
LOCATE VALUE NOWAIT(1),RANGE 0 TO 3,AT 1:SCALE
LOCATE VALUE NOWAIT(2),RANGE -RANGE TO RANGE,AT 0:XMOVE
LOCATE VALUE NOWAIT(3),RANGE -RANGE TO RANGE,AT 0:YMOVE
LOCATE VALUE NOWAIT(4),RANGE -RANGE TO RANGE,AT 0:ZMOVE
LOCATE VALUE NOWAIT(5),RANGE -RANGE TO RANGE,AT 0:WMOVE
LOCATE VALUE NOWAIT(6),RANGE -360 TO 360,AT 0:XY
LOCATE VALUE NOWAIT(7),RANGE -360 TO 360,AT 0:XZ
LOCATE VALUE NOWAIT(8),RANGE -360 TO 360,AT 0:YU
LOCATE VALUE NOWAIT(9),RANGE -360 TO 360,AT 0:YZ
LOCATE VALUE NOWAIT(10),RANGE -360 TO 360,AT 0:XU
LOCATE VALUE NOWAIT(11),RANGE -360 TO 360,AT 0:ZU
DO
LOCATE VALUE NOWAIT(1):SCALE
LOCATE VALUE NOWAIT(2):XMOVE
LOCATE VALUE NOWAIT(3):YMOVE
LOCATE VALUE NOWAIT(4):ZMOVE
LOCATE VALUE NOWAIT(5):WMOVE
LOCATE VALUE NOWAIT(6):XY
LOCATE VALUE NOWAIT(7):XZ
LOCATE VALUE NOWAIT(8):YU
LOCATE VALUE NOWAIT(9):YZ
LOCATE VALUE NOWAIT(10):XU
LOCATE VALUE NOWAIT(11):ZU
MAT ROTXY=IDN !'XY平面上の回転(ZU平面を軸とする回転)
LET ROTXY(1,1)=COS(XY)
LET ROTXY(1,2)=SIN(XY)
LET ROTXY(2,1)=-SIN(XY)
LET ROTXY(2,2)=COS(XY)
MAT ROTXZ=IDN !'XZ平面上の回転(YU平面を軸とする回転)
LET ROTXZ(1,1)=COS(XZ)
LET ROTXZ(1,3)=-SIN(XZ)
LET ROTXZ(3,1)=SIN(XZ)
LET ROTXZ(3,3)=COS(XZ)
MAT ROTYU=IDN !'YU平面上の回転(XZ平面を軸とする回転)
LET ROTYU(2,2)=COS(YU)
LET ROTYU(2,4)=-SIN(YU)
LET ROTYU(4,2)=SIN(YU)
LET ROTYU(4,4)=COS(YU)
MAT ROTYZ=IDN !'YZ平面上の回転(XU平面を軸とする回転)
LET ROTYZ(2,2)=COS(YZ)
LET ROTYZ(2,3)=SIN(YZ)
LET ROTYZ(3,2)=-SIN(YZ)
LET ROTYZ(3,3)=COS(YZ)
MAT ROTXU=IDN !'XU平面上の回転(YZ平面を軸とする回転)
LET ROTXU(1,1)=COS(XU)
LET ROTXU(1,4)=SIN(XU)
LET ROTXU(4,1)=-SIN(XU)
LET ROTXU(4,4)=COS(XU)
MAT ROTZU=IDN !'ZU平面上の回転(XY平面を軸とする回転)
LET ROTZU(3,3)=COS(ZU)
LET ROTZU(3,4)=-SIN(ZU)
LET ROTZU(4,3)=SIN(ZU)
LET ROTZU(4,4)=COS(ZU)
MAT M=ROTXY*ROTXZ*ROTYU*ROTYZ*ROTXU*ROTZU
SET DRAW MODE HIDDEN
CLEAR
FOR J=0 TO MM-1
FOR I=0 TO NN-1
CALL PLOT(XX(I,J),YY(I,J),ZZ(I,J),WW(I,J))
CALL PLOT(XX(I,J+1),YY(I,J+1),ZZ(I,J+1),WW(I,J+1))
CALL PLOT(XX(I+1,J+1),YY(I+1,J+1),ZZ(I+1,J+1),WW(I+1,J+1))
CALL PLOT(XX(I+1,J),YY(I+1,J),ZZ(I+1,J),WW(I+1,J))
CALL PLOT(XX(I,J),YY(I,J),ZZ(I,J),WW(I,J))
PLOT LINES
NEXT I
NEXT J
IF FL=0 THEN
SET WINDOW -LMAX*1.2,LMAX*1.2,-LMAX*1.2,LMAX*1.2
LET FL=1
END IF
SET DRAW MODE EXPLICIT
LOOP

SUB PLOT(X,Y,Z,W)
LET LW=2
LET POINT(1)=X-XMOVE
LET POINT(2)=Y-YMOVE
LET POINT(3)=Z-ZMOVE
LET POINT(4)=W-WMOVE
MAT POINT=POINT*M
! FOR K=1 TO 3
! LET POINT(K)=POINT(K)/(LW-POINT(4))
! NEXT
IF FL=0 THEN
LET LMIN=MIN(LMIN,POINT(1))
LET LMAX=MAX(LMAX,POINT(1))
LET LMIN=MIN(LMIN,POINT(2))
LET LMAX=MAX(LMAX,POINT(2))
ELSE
PLOT LINES:POINT(1)*SCALE,POINT(2)*SCALE;
END IF
END SUB
END

Re: PLOT LINESの使い方 - SHIRAISHI Kazuo

2026/06/02 (Tue) 08:06:12

LET nMax=100
LET dwt=2*PI/nMax
FOR k=0 TO nMax
LET wt=k*dwt
PRINT k, wt-PI
NEXT K
END
を実行してみるとわかりますが，wt=PIになることがないのが原因のような気がします。
wtの値で分類するのでなく，kの値で分類すればうまくいくと思います。
(補足)
組込関数PIは16桁を超える精度を持っていますが，数値変数wtは15桁の精度で丸められます。
具体的には，PI=3.141592653589793238ですが，
k=50のとき，wt=3.1415926535898です。
だから，k=50のとき，wt>PIです。

Re: PLOT LINESの使い方 - 島村1243

2026/06/02 (Tue) 10:38:53

SHRAISHI先生、早速のご教示有難う御座いました。
ご教示に従って、For Next を下記の様に修正したら、目的を達成できました!!!

FOR k=0 TO nMax
LET wt=k*dwt
IF k<nMax/2 THEN PLOT LINES:wt,wt/PI;
IF k=nMax/2 THEN !---縦垂直線の処理
PLOT LINES:wt,1;
PLOT LINES:wt,0;
END IF
IF k>nMax/2 THEN PLOT LINES:wt,(wt-PI)/PI;
NEXT K

Re: 文字探しゲーム - しばっち

2026/05/31 (Sun) 07:42:41

指定された文字を探し出しクリックします。
全部見つけたと思ったら「完了」をクリックすると判定されます。

難易度を上げる為に回転文字や反転文字等も含まれています。


OPTION ANGLE DEGREES
RANDOMIZE
LET COL=15
LET ROW=10
DIM A$(200),M$(COL,ROW),R(COL,ROW),MOJI$(5),SX(COL,ROW),SY(COL,ROW)
DO
READ IF MISSING THEN EXIT DO:D$
FOR J=1 TO LEN(D$)
LET K=K+1
LET A$(K)=MID$(D$,J,1)
NEXT J
LOOP
CALL GINIT(600,600)
LET N=INT(RND*4+2)
FOR I=1 TO N
LET MOJI$(I)=A$(RND*K+1)
NEXT I
DATA "あいうえおかがきぎくぐけげこごさざしじすずせぜそぞただちぢつづてでとどなにぬねのはばぱひびぴふぶぷへべぺほぼぽまみむめもやゆよらりるれろわゐゑをん"
DATA "アイウヴエオカガキギクグケゲコゴサザシジスズセゼソゾタダチヂツヅテデトドナニヌネノハバパヒビピフブプヘベペホボポマミムメモヤユヨラリルレロワヰヱヲン"
LET H=80
SET TEXT JUSTIFY "LEFT","TOP"
SET TEXT COLOR 7
SET TEXT HEIGHT H*.8
FOR I=0 TO N-1
PLOT TEXT ,AT 50+100*I,500:MOJI$(I+1)
NEXT I
LET H=MIN(600/COL,600/ROW)
SET TEXT HEIGHT H*.8
PLOT TEXT ,AT 50,200:"を探してください"
DO
MOUSE POLL X,Y,LEFT,RIGHT
LOOP UNTIL RIGHT=1 OR LEFT=1
DO
MOUSE POLL X,Y,LEFT,RIGHT
LOOP WHILE RIGHT=1 OR LEFT=1
CLEAR
SET TEXT JUSTIFY "CENTER","HALF"
FOR J=0 TO ROW-1
FOR I=0 TO COL-1
IF RND<.4 THEN LET S$=MOJI$(INT(RND*N)+1) ELSE LET S$=A$(INT(RND*K)+1)
LET M$(I+1,J+1)=S$
LET R(I+1,J+1)=90*INT(RND*4)
LET SX(I+1,J+1)=(-1)^INT(2*RND)
LET SY(I+1,J+1)=(-1)^INT(2*RND)
DRAW WORD(S$) WITH ROTATE(R(I+1,J+1))*SCALE(SX(I+1,J+1),SY(I+1,J+1))*SHIFT(H/2+I*H,H/2+J*H+200)
NEXT I
NEXT J
SET TEXT JUSTIFY "LEFT","TOP"
CALL BOX(450,100,450+H*2,50,7)
PLOT TEXT,AT 450,100:"完了"
SET TEXT BACKGROUND "OPAQUE"
LET T=TIME
DO
SET TEXT JUSTIFY "LEFT","TOP"
SET TEXT COLOR 7
PLOT TEXT,AT 0,100:"MISS "&USING$("%%%",MISS)
PLOT TEXT,AT 0,50:"SCORE "&USING$("%%%",HIT)
SET TEXT JUSTIFY "CENTER","HALF"
DO
MOUSE POLL X,Y,LEFT,RIGHT
LOOP UNTIL RIGHT=1 OR LEFT=1
DO
MOUSE POLL X,Y,LEFT,RIGHT
LOOP WHILE RIGHT=1 OR LEFT=1
IF X>450 AND X<530 AND Y>50 AND Y<100 THEN ! 完了をクリックしたら
FOR J=0 TO ROW-1
FOR I=0 TO COL-1
FOR L=1 TO N
IF M$(I+1,J+1)=MOJI$(L) THEN
SET TEXT COLOR 4
DRAW WORD(M$(I+1,J+1)) WITH ROTATE(R(I+1,J+1))*SCALE(SX(I+1,J+1),SY(I+1,J+1))*SHIFT(H/2+I*H,H/2+200+J*H)
LET MISS=MISS+1
END IF
NEXT L
NEXT I
NEXT J
SET TEXT JUSTIFY "LEFT","TOP"
SET TEXT COLOR 7
PLOT TEXT,AT 0,150:"TIME "&USING$("%%",INT((TIME-T)/60)) &":"&USING$("%%",MOD(INT(TIME-T),60))
PLOT TEXT,AT 0,100:"MISS "&USING$("%%%",MISS)
PLOT TEXT,AT 0,50:"SCORE "&USING$("%%%",HIT)
STOP
END IF
LET XX=INT(X/H)
LET YY=INT((Y-200)/H)
LET FL=0
FOR L=1 TO N
IF M$(XX+1,YY+1)=MOJI$(L) THEN LET FL=1
NEXT L
IF FL=1 THEN
SET TEXT COLOR 0
DRAW WORD(M$(XX+1,YY+1)) WITH ROTATE(R(XX+1,YY+1))*SCALE(SX(XX+1,YY+1),SY(XX+1,YY+1))*SHIFT(H/2+XX*H,200+H/2+YY*H)
LET M$(XX+1,YY+1)=""
LET HIT=HIT+1
ELSE
SET TEXT COLOR 2
DRAW WORD(M$(XX+1,YY+1)) WITH ROTATE(R(XX+1,YY+1))*SCALE(SX(XX+1,YY+1),SY(XX+1,YY+1))*SHIFT(H/2+XX*H,200+H/2+YY*H)
LET M$(XX+1,YY+1)=""
LET MISS=MISS+1
END IF
LOOP
END

EXTERNAL SUB GINIT(XSIZE,YSIZE)
SET BITMAP SIZE XSIZE,YSIZE
SET WINDOW 0,XSIZE-1,0,YSIZE-1
SET POINT STYLE 1
SET COLOR MODE "REGULAR"
SET COLOR MIX(0) 0,0,0
SET COLOR MIX(1) 0,0,1
SET COLOR MIX(2) 1,0,0
SET COLOR MIX(3) 1,0,1
SET COLOR MIX(4) 0,1,0
SET COLOR MIX(5) 0,1,1
SET COLOR MIX(6) 1,1,0
SET COLOR MIX(7) 1,1,1
CLEAR
END SUB

EXTERNAL SUB BOXFULL(X1,Y1,X2,Y2,C)
SET COLOR C
PLOT AREA:X1,Y1;X2,Y1;X2,Y2;X1,Y2;X1,Y1
END SUB

EXTERNAL SUB LINE(XS,YS,XE,YE,C)
SET COLOR C
PLOT LINES
PLOT LINES:XS,YS;XE,YE
END SUB

EXTERNAL SUB BOX(XS,YS,XE,YE,C)
CALL LINE(XS,YS,XE,YS,C)
CALL LINE(XE,YS,XE,YE,C)
CALL LINE(XE,YE,XS,YE,C)
CALL LINE(XS,YE,XS,YS,C)
END SUB

EXTERNAL PICTURE WORD(A$)
PLOT LETTERS ,AT 0,0:A$
END PICTURE

Re: 三角関数値が異常です - SHIRAISHI Kazuo

2026/04/12 (Sun) 08:53:30

ご報告ありがとうございした。
原因は，Intel FPUのFSIN，FCOS命令をそのまま使っていたことです。
マニュアルをよく読むと，引数は -2^63 ～ 2^63 の範囲でなければならないようです。
倍精度数の精度は53ビットなので，その範囲外の数値に対して三角関数の値を求めることは無意味ですが，どう扱うのが適切か，検討します。

Re: 三角関数値が異常です - しばっち

2026/04/19 (Sun) 07:55:10

周期性があるのでMOD関数を利用すれば計算自体は可能ですが
誤差が累積していくようです

https://ja.wikipedia.org/wiki/ワイエルシュトラス関数

LET B=23
LET X=.4
FOR I=1 TO 15
LET C=B^I*PI*X
LET X=COS(C)
LET Y=COS(MOD(C,2*PI))
PRINT I;X;Y;X-Y
NEXT I
END


実行結果


1 -.809016994374948 -.809016994374948 -5.55111512312578E-16
2 .99555901481267 .995559014812675 -4.88498130835069E-15
3 -.994477987872487 -.994477987872643 1.55653268052447E-13
4 .624271813435692 .624271813462338 -2.66459077025161E-11
5 5.18260078207372E-2 .051826008312573 -4.91835808191787E-10
6 -.903399056291523 -.903399055889035 -4.02487820849728E-10
7 .47391871819713 .473919049744119 -3.31546989085219E-7
8 -.922909753831513 -.922908005104853 -1.74872665914538E-6
9 .791479622428268 .791603979074663 -1.24356646394985E-4
10 -.556013013484207 -.552670804109104 -3.34220937510299E-3
11 .999592310690031 .999339079455446 2.53231234584206E-4
12 .774903379838499 -.974661220596496 1.74956460043499
13 -.526222536776585 .960170286650366 -1.48639282342695
14 -1.9165008703538E+19 -.653643620863612 -1.9165008703538E+19
15 -1.6053747716131E+40 1 -1.6053747716131E+40

Re: 三角関数値が異常です - SHIRAISHI Kazuo

2026/04/19 (Sun) 09:21:09

Delphi7のsin, cosは引数の絶対値が2^63を超えるとどちらも値が0になります。（十進BASIC Ver.7 はアセンブラでFSIN, FCOSを使用していました）
lazarusで用いられるFPC (Free Pascal) のsin, cosは，Intel FPUを使う環境ではintel FPUの仕様そのままでした。
十進BASIC Ver.7，Ver.8の修正版では，2πで割った余りを求めてから三角関数の値を求めることにしました。
ただし，ここでのπは，64ビット精度の拡張精度数です。倍精度より少し高精度ですが，真の値ではありません。
引数の絶対値が大きいときでも計算結果は-1～1の範囲に納まりますが，数値自体はほぼ無意味です。
BASICAccとParact BASICは，当面，手を加えずFPCの計算に委ねることにしました。
なお，2πでの剰余を精密に計算したいときは，十進1000桁モードを使用してください。十進1000桁モードだとPIは1000桁超の精度を持ちます。

Re: 三角関数値が異常です - SHIRAISHI Kazuo

2026/04/19 (Sun) 09:30:19

十進1000桁モードでの実行結果です。
OPTION ARITHMETIC DECIMAL_HIGH
LET B=15
LET X=1
FOR N=1 TO 30
PRINT N;COS(B^N*X*PI);SIN(B^N*X*PI)
NEXT N
END

1 -1 -.00000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000092857116
2 -1 -.0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000139285674
3 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000208928511
4 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000003133927665
5 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000047008914975
6 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000705133724625
7 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000010577005869375
8 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000158655088040625
9 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000002379826320609375
10 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000035697394809140625
11 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000535460922137109375
12 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000008031913832056640625
13 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000120478707480849609375
14 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001807180612212744140625
15 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000027107708454675537109375
16 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000406615637747867431640625
17 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000006099234566218011474609375
18 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000091488518493270172119140625
19 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001372327777399052581787109375
20 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000020584916660985788726806640625
21 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000308773749914786830902099609375
22 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000004631606248721802463531494140625
23 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000069474093594296325474456787109375
24 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001042111405862226673116851806640625
25 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000015631671089231921290752777099609375
26 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000234475066341508702147291656494140625
27 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000003517125995123063372607374847412109375
28 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000052756889926846599849707622711181640625
29 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000791353348902701594788002340667724609375
30 -1 -.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000011870300233540525653181627110015869140625

この結果が物語ることは，B^N*X*PIの計算結果と，sinを計算するときに十進BASIC内部で用いるπの精度の違いです。
十進BASICサンプルプログラムのQ&A\SIN_PI.BASも参照してください。
要するに，PI自体がπの近似値であること，そして，SINの計算精度が保証される範囲自体も狭いことです。
十進1000桁モードの機械最小値を1E-1017に設定しているのも，精度保証との関係です。

Re: 三角関数値が異常です - SHIRAISHI Kazuo

2026/04/19 (Sun) 10:47:08

2進モードでのPIは，およそ3.141592653589793115999782457です。
この値を元にπの整数倍を計算していくと，10^16倍したとき，およそ-1.2ほどの誤差がでてきます。三角関数の値を求めるのには不適切な誤差です。
10進1000桁モードで計算してみてください。

OPTION ARITHMETIC Decimal_High
LET p=3.141592653589793115999782457
FOR i=0 TO 16
PRINT 10^i*p - 10^i*PI
NEXT i
END

Re: 三角関数値が異常です - SHIRAISHI Kazuo

2026/04/19 (Sun) 14:55:24

三角関数の正確な値を求める基本は，xを2πで割った余りを正確に求めることです。
xを2*PIで割った商をq，余りをrとして，PI=π+eだとすると，
x=q*2*PI+r
PI=π+e
なので，
x=q*2*π＋q*2*e+r
だから，xを2πで割った余りはq*2*e+rです。

Re: 三角関数値が異常です - SHIRAISHI Kazuo

2026/04/22 (Wed) 07:57:17

PI≒πの近似誤差を気にせず計算したいときは，度を単位に計算することが推奨できます。
OPTION ANGLE DEGREES
FOR k=50 TO 60
LET x=180*2^k
PRINT k,x,COS(x),SIN(x)
NEXT k
END
2進モードでの実行結果は
50 2.02661983231672E+17 1 0
51 4.05323966463345E+17 1 0
52 8.10647932926689E+17 1 0
53 1.62129586585338E+18 1 0
54 3.24259173170676E+18 1 0
55 6.48518346341351E+18 1 0
56 1.2970366926827E+19 1 0
57 2.59407338536541E+19 1 0
58 5.18814677073081E+19 .999390827019096 0
59 1.03762935414616E+20 .999390827019096 0
60 2.07525870829232E+20 .994521895368273 0
30度，45度など，度を単位に切りのいい数に対する三角関数の値を求めたいときはこちらが推奨できます。
2進モードだと，2^53までの整数は正確に表現できます。

Re: 三角関数値が異常です - SHIRAISHI Kazuo

2026/04/22 (Wed) 08:06:24

上のプログラムで，K=58以降でCOS(x)が1でない値になり，さらにその先まで計算していくとk=64以降で0になります。これは，十進BASICのバグ(?)です。
COS(x)=SIN(x+90)として計算しているのが原因で，x+90の計算が誤差を持つような範囲で正しい値が求まらなくなっています。

Re: 三角関数値が異常です - SHIRAISHI Kazuo

2026/04/24 (Fri) 17:24:18

https://ja.wikipedia.org/wiki/ワイエルシュトラス関数
で，
COS(b^n*x*π)を正確に計算したいのであれば，b^n*xを2で割った余りをrとおくと，COS(b^n*x*π)=COS(r*π)
0≦r＜2なので，COS(r*π)をCOS(r*PI)で近似してもさほど問題ない。


OPTION ARITHMETIC NATIVE
LET a=1/2
LET b=11
FUNCTION w(x)
LET t=0
LET n=0
FOR n=200 TO 0 STEP -1
LET r=MOD(b^n*x,2)
LET t=t+a^n*COS(r*PI)
NEXT n
LET w=t
END FUNCTION
SET POINT STYLE 1
SET WINDOW -2,2,-2,2
DRAW grid
FOR x=-2 TO 2 STEP 0.00001
PLOT LINES:x,w(x);
NEXT x
END

Re: 三角関数値が異常です - SHIRAISHI Kazuo

2026/04/25 (Sat) 07:25:56

上掲のプログラムで，w(x)≒0付近の相対誤差ではなく，絶対誤差のみを要求するとき，FOR文の始値n=200は過大です。
COSの値は―1～1の範囲になるのでw(x)の和は公比a=1/2の等比級数の和で抑えられます。
たとえば，n=54に選べば，計算しない残余項の和は2^(-54)以下になります。
なお，n≦54であっても，nが少し大きな数だと，b^n*xはかなり大きな数になって誤差をはらみ，r=MOD(b^n*x,2)は多くの場合，0になります。
そのあたりまで考えて正確に計算したいのであれば，MOD(b^n*x,2)を正確に導き出す手法を編み出す必要があります。

Re: 三角関数値が異常です - SHIRAISHI Kazuo

2026/04/25 (Sat) 09:45:10

MOD(b^n*x,2) = MOD((b^(n-1)*x)*b,2) = MOD(MOD(b^(n-1)*x,2)*b,2) を使えばよさそうです。

LET b=11
FOR x=-2 TO 2 STEP 2^(-10)
LET r=MOD(x,2)
FOR n=1 TO 60
LET r=MOD(r*b,2)
PRINT x,n,r
NEXT n
NEXT x
END

xに2進小数として切りのいい数を指定しないと，カオスになります。

Re: 三角関数値が異常です - SHIRAISHI Kazuo

2026/05/03 (Sun) 08:36:57

上掲プログラムで，xに2進小数として切りのいい数を指定しないとカオスになるのは，2進モードのMOD関数が十分な精度を持たないのが影響していました。修正します。

Re: 三角関数値が異常です - SHIRAISHI Kazuo

2026/05/16 (Sat) 09:10:26

ワイエルシュトラス関数のグラフを描くプログラムです。
グラフを描く目的ならさほど多くの項の計算は不要です。
2進倍精度で計算するのであれば，140行で求めているNmaxまで計算することが目安になります。

100 REM ワイエルシュトラス関数
110 OPTION ARITHMETIC NATIVE
120 LET a=0.5
130 LET b=13
140 LET Nmax=-54*LOG2(a)
150 FUNCTION w(x)
160 LET t=0
170 LET r=MOD(x,2)
180 FOR n=0 TO NMax
190 LET t=t+a^n*COS(r*PI)
200 LET r=MOD(r*b,2)
210 NEXT n
220 LET w=t
230 END FUNCTION
240 SET POINT STYLE 1
250 SET WINDOW -2,2,-2,2
260 DRAW GRID
270 FOR x=-2 TO 2 STEP 2^(-18)
280 PLOT POINTS:x,w(x)
290 NEXT x
300 END


https://ja.wikipedia.org/wiki/ワイエルシュトラス関数
にある画像は，a=0.5，b=3 程度で作成されたもののように思えます。

Re: ガチャ - しばっち

2026/05/03 (Sun) 07:43:01

https://dskjal.com/statistics/chance-calculator-jp
https://mikami3345.cloudfree.jp/Probability/Probability.html

! スマホガチャ　当選率P%をN回ガチャ引いて1個以上当たる確率
!LET P=1/100 ! 当選率
!LET N=100 ! ガチャを引く回数
!PRINT (1-(1-P)^N)*100;"%"

SET WINDOW -10,150,-10,100
DRAW GRID(10,10)
LET P=1/100 ! 当選率
FOR N=0 TO 150 ! ガチャを引く回数
PLOT LINES:N,(1-(1-P)^N)*100;
NEXT N
END
-----------------------------------------------------------------------
上記のシュミレーションプログラム


RANDOMIZE
DIM GACHA(100)
INPUT PROMPT "当選確率 (%)(1-99) ":L
INPUT PROMPT "ガチャを引く回数 ":N
LET M=500 ! 試行回数
FOR I=1 TO L ! 100個のうちL個が当たり
LET GACHA(I)=1
NEXT I
FOR I=1 TO M
FOR J=1 TO 100 ! シャッフル
SWAP GACHA(J),GACHA(INT(RND*100+1))
NEXT J
LET FLG=0
FOR J=1 TO N
LET P=INT(RND*100+1) ! ガチャを引く
IF GACHA(P)=1 THEN ! 当たりを引いた
LET FLG=1
EXIT FOR
END IF
NEXT J
IF FLG=0 THEN LET C=C+1 ! 全部外れたらカウント
NEXT I
PRINT "当選確率";L;"%の時、";N;"回ガチャを引いて当たる確率";(M-C)/M*100;"%"
END

　　　　　　　　実行結果


当選確率 (%)(1-99) 1
ガチャを引く回数 100
当選確率 1 %の時、 100 回ガチャを引いて当たる確率 62.2 %
-----------------------------------------------------------------------
当選確率L%の時、N回ガチャを引いて当たる回数の確率


OPTION BASE 0
RANDOMIZE
DIM GACHA(100),C(100)
INPUT PROMPT "当選確率 (%)(1-99) ":L ! L=1
INPUT PROMPT "ガチャを引く回数 ":N ! N=100
LET M=500 ! 試行回数
FOR I=1 TO L ! 100個のうちL個が当たり
LET GACHA(I)=1
NEXT I
FOR I=1 TO M
FOR J=1 TO 100 ! シャッフル
SWAP GACHA(J),GACHA(INT(RND*100+1))
NEXT J
LET HIT=0
FOR J=1 TO N
LET P=INT(RND*100+1) ! ガチャを引く
IF GACHA(P)=1 THEN ! 当たりを引いた
LET HIT=HIT+1
END IF
NEXT J
LET C(HIT)=C(HIT)+1
NEXT I
FOR I=1 TO 100
IF C(I)>0 THEN PRINT "当選確率";L;"%の時、";N;"回ガチャを引いて";I;"回当たる確率";C(I)/M*100;"%"
NEXT I
END

　　　　　　　　　　実行結果


当選確率 (%)(1-99) 1
ガチャを引く回数 100
当選確率 1 %の時、 100 回ガチャを引いて 1 回当たる確率 34.2 %
当選確率 1 %の時、 100 回ガチャを引いて 2 回当たる確率 21.4 %
当選確率 1 %の時、 100 回ガチャを引いて 3 回当たる確率 5.8 %
当選確率 1 %の時、 100 回ガチャを引いて 4 回当たる確率 1.4 %
当選確率 1 %の時、 100 回ガチャを引いて 5 回当たる確率 .2 %
当選確率 1 %の時、 100 回ガチャを引いて 9 回当たる確率 .2 %
-----------------------------------------------------------------------
入手確率 X%を達成するために必要な試行回数


LET P=1/100 ! 当選率
LET N=99/100 ! 入手確率
FOR X=1 TO 10000
IF 1-(1-P)^X>N THEN EXIT FOR
NEXT X
PRINT "当選確率";P*100;"% 入手確率";N*100;"% 達するために必要な試行回数";X;"回"
END
-----------------------------------------------------------------------
N回ガチャ引いても当たらない確率


LET N=100 ! 試行回数
LET P=.5/100 ! 当選確率
PRINT "当選確率";P*100;"%の時、";N;"回ガチャを引いても当たらない確率";(1-P)^N*100;"%"
END
-----------------------------------------------------------------------
当たりが2個ある場合の確率

但し、当選確率1%という時、合わせての合計確率1/200+1/200で1/100つまり計1%と
それぞれが1%、つまり1/100+1/100で計2%となる場合の2通りがある。


RANDOMIZE
DIM GACHA(200)
LET M=500 ! 試行回数
LET MODE=1
SELECT CASE MODE
CASE 1
INPUT PROMPT "合わせての当選確率 (%)(1-99) ":L
FOR I=1 TO L ! 200個のうちL個が当たり L/200+L/200=L/100(L%)
LET C=C+1
LET GACHA(C)=1
LET C=C+1
LET GACHA(C)=2
NEXT I
CASE 2
INPUT PROMPT "個別の当選確率 (%)(1-99) ":L
FOR I=1 TO L*2 ! 200個のうち2*L個が当たり L/100(L%)+L/100(L%)=2*L/100 (2L%)
LET C=C+1
LET GACHA(C)=1
LET C=C+1
LET GACHA(C)=2
NEXT I
!CASE 3
! INPUT PROMPT "1,2合わせての当選確率 (%)(1-99) ":L1,L2
! FOR I=1 TO L1 ! 200個のうちL1個とL2個が当たり L1/200+L2/200
! LET C=C+1
! LET GACHA(C)=1
! NEXT I
! FOR I=1 TO L2
! LET C=C+1
! LET GACHA(C)=2
! NEXT I
!CASE 4
! INPUT PROMPT "1,2個別の当選確率 (%)(1-99) ":L1,L2
! FOR I=1 TO L1*2 ! 200個のうちL1個とL2個が当たり L1/100+L2/100
! LET C=C+1
! LET GACHA(C)=1
! NEXT I
! FOR I=1 TO L2*2
! LET C=C+1
! LET GACHA(C)=2
! NEXT I
END SELECT
INPUT PROMPT "ガチャを引く回数 ":N
FOR I=1 TO M
FOR J=1 TO 200 ! シャッフル
SWAP GACHA(J),GACHA(INT(RND*200+1))
NEXT J
LET HIT1=0
LET HIT2=0
FOR J=1 TO N
LET P=INT(RND*200+1) ! ガチャを引く
IF GACHA(P)=1 THEN LET HIT1=HIT1+1
IF GACHA(P)=2 THEN LET HIT2=HIT2+1
NEXT J
IF (HIT1>0 AND HIT2=0) OR (HIT1=0 AND HIT2>0) THEN LET EITHER=EITHER+1
IF HIT1>0 AND HIT2>0 THEN LET BOTH=BOTH+1
IF HIT1=0 AND HIT2=0 THEN LET MISS=MISS+1
NEXT I
IF L>0 THEN
PRINT "当選確率";L;"%の時、";N;"回ガチャを引いてどちらかが当たる確率";EITHER/M*100;"%"
PRINT "当選確率";L;"%の時、";N;"回ガチャを引いてどちらも当たる確率";BOTH/M*100;"%"
PRINT "当選確率";L;"%の時、";N;"回ガチャを引いてどちらも外れる確率";MISS/M*100;"%"
ELSE
PRINT "当選確率";L1;"%、";L2;"%の時、";N;"回ガチャを引いてどちらかが当たる確率";EITHER/M*100;"%"
PRINT "当選確率";L1;"%、";L2;"%の時、";N;"回ガチャを引いてどちらも当たる確率";BOTH/M*100;"%"
PRINT "当選確率";L1;"%、";L2;"%の時、";N;"回ガチャを引いてどちらも外れる確率";MISS/M*100;"%"
END IF
END

　　　　　　　　　実行結果


合わせての当選確率 (%)(1-99) 1
ガチャを引く回数 100
当選確率 1 %の時、 100 回ガチャを引いてどちらかが当たる確率 46.2 %
当選確率 1 %の時、 100 回ガチャを引いてどちらも当たる確率 15 %
当選確率 1 %の時、 100 回ガチャを引いてどちらも外れる確率 38.8 %

-----------------------------------------------------------------------
当たりが3個ある場合の確率

RANDOMIZE
DIM GACHA(300)
LET M=500 ! 試行回数
LET MODE=1
SELECT CASE MODE
CASE 1
INPUT PROMPT "合わせての当選確率 (%)(1-99) ":L
FOR I=1 TO L ! 300個のうちL個が当たり L/300+L/300+L/300=L/100(L%)
LET C=C+1
LET GACHA(C)=1
LET C=C+1
LET GACHA(C)=2
LET C=C+1
LET GACHA(C)=3
NEXT I
CASE 2
INPUT PROMPT "個別の当選確率 (%)(1-99) ":L
FOR I=1 TO L*3 ! 300個のうち3*L個が当たり L/100(L%)+L/100(L%)+L/100(L%)=L/50 (3L%)
LET C=C+1
LET GACHA(C)=1
LET C=C+1
LET GACHA(C)=2
LET C=C+1
LET GACHA(C)=3
NEXT I
!CASE 3
! INPUT PROMPT "1,2,3合わせての当選確率 (%)(1-99) ":L1,L2,L3
! FOR I=1 TO L1 ! 300個のうちL個が当たり L1/300+L2/300+L3/300
! LET C=C+1
! LET GACHA(C)=1
! NEXT I
! FOR I=1 TO L2
! LET C=C+1
! LET GACHA(C)=2
! NEXT I
! FOR I=1 TO L3
! LET C=C+1
! LET GACHA(C)=3
! NEXT I
!CASE 4
! INPUT PROMPT "1,2,3個別の当選確率 (%)(1-99) ":L1,L2,L3
! FOR I=1 TO L1*3 ! 300個のうちL個が当たり L1/100+L2/100+L3/100
! LET C=C+1
! LET GACHA(C)=1
! NEXT I
! FOR I=1 TO L2*3
! LET C=C+1
! LET GACHA(C)=2
! NEXT I
! FOR I=1 TO L3*3
! LET C=C+1
! LET GACHA(C)=3
! NEXT I
END SELECT
INPUT PROMPT "ガチャを引く回数 ":N
FOR I=1 TO M
FOR J=1 TO 300 ! シャッフル
SWAP GACHA(J),GACHA(INT(RND*300+1))
NEXT J
LET HIT1=0
LET HIT2=0
LET HIT3=0
FOR J=1 TO N
LET P=INT(RND*300+1) ! ガチャを引く
IF GACHA(P)=1 THEN LET HIT1=HIT1+1
IF GACHA(P)=2 THEN LET HIT2=HIT2+1
IF GACHA(P)=3 THEN LET HIT3=HIT3+1
NEXT J
IF HIT1>0 AND HIT2=0 AND HIT3=0 THEN LET EITHER1=EITHER1+1
IF HIT1=0 AND HIT2>0 AND HIT3=0 THEN LET EITHER1=EITHER1+1
IF HIT1=0 AND HIT2=0 AND HIT3>0 THEN LET EITHER1=EITHER1+1

IF HIT1>0 AND HIT2>0 AND HIT3=0 THEN LET EITHER2=EITHER2+1 ! 1,2
IF HIT1>0 AND HIT3>0 AND HIT2=0 THEN LET EITHER2=EITHER2+1 ! 1,3
IF HIT2>0 AND HIT3>0 AND HIT1=0 THEN LET EITHER2=EITHER2+1 ! 2,3

IF HIT1>0 AND HIT2>0 AND HIT3>0 THEN LET ALL=ALL+1
IF HIT1=0 AND HIT2=0 AND HIT3=0 THEN LET MISS=MISS+1
NEXT I
IF L>0 THEN
PRINT "当選確率";L;"%の時、";N;"回ガチャを引いてどれか1種類が当たる確率";EITHER1/M*100;"%"
PRINT "当選確率";L;"%の時、";N;"回ガチャを引いてどれか2種類が当たる確率";EITHER2/M*100;"%"
PRINT "当選確率";L;"%の時、";N;"回ガチャを引いて全種類が当たる確率";ALL/M*100;"%"
PRINT "当選確率";L;"%の時、";N;"回ガチャを引いてどれも外れる確率";MISS/M*100;"%"
ELSE
PRINT "当選確率";L1;"%、";L2;"%、";L3;"%の時、";N;"回ガチャを引いてどれか1種類が当たる確率";EITHER1/M*100;"%"
PRINT "当選確率";L1;"%、";L2;"%、";L3;"%の時、";N;"回ガチャを引いてどれか2種類が当たる確率";EITHER2/M*100;"%"
PRINT "当選確率";L1;"%、";L2;"%、";L3;"%の時、";N;"回ガチャを引いて全種類が当たる確率";ALL/M*100;"%"
PRINT "当選確率";L1;"%、";L2;"%、";L3;"%の時、";N;"回ガチャを引いてどれも外れる確率";MISS/M*100;"%"
END IF
END

　　　　　　　　実行結果


合わせての当選確率 (%)(1-99) 1
ガチャを引く回数 100
当選確率 1 %の時、 100 回ガチャを引いてどれか1種類が当たる確率 44.4 %
当選確率 1 %の時、 100 回ガチャを引いてどれか2種類が当たる確率 17.4 %
当選確率 1 %の時、 100 回ガチャを引いて全種類が当たる確率 2.2 %
当選確率 1 %の時、 100 回ガチャを引いてどれも外れる確率 36 %
-----------------------------------------------------------------------
当たりが複数の場合


OPTION BASE 0
RANDOMIZE
INPUT PROMPT "当たりの種類 ":K
INPUT PROMPT "ガチャを引く回数 ":N
DIM GACHA(100*K),C(K,N),HIT(K),LL(K)
LET M=500 ! 試行回数
LET MODE=1
SELECT CASE MODE
CASE 1
FOR J=1 TO K
INPUT PROMPT STR$(J)&"の当選確率 (%)(1-99) ":LL(J)
NEXT J
FOR I=1 TO K
FOR J=1 TO LL(I)
LET CC=CC+1
LET GACHA(CC)=I
NEXT J
NEXT I
CASE 2
FOR J=1 TO K
INPUT PROMPT STR$(J)&"の当選確率 (%)(1-99) ":LL(J)
NEXT J
FOR I=1 TO K
FOR J=1 TO LL(I)*K
LET CC=CC+1
LET GACHA(CC)=I
NEXT J
NEXT I
END SELECT
FOR I=1 TO M
FOR J=1 TO 100*K ! シャッフル
SWAP GACHA(J),GACHA(INT(RND*100*K+1))
NEXT J
MAT HIT=ZER
FOR J=1 TO N
LET P=INT(RND*100*K+1) ! ガチャを引く
FOR F=1 TO K
IF GACHA(P)=F THEN
LET HIT(F)=HIT(F)+1 ! F種目の当たり引いた
EXIT FOR
END IF
NEXT F
NEXT J
FOR F=1 TO K
LET C(F,HIT(F))=C(F,HIT(F))+1
NEXT F
NEXT I
FOR I=1 TO N
FOR F=1 TO K
IF C(F,I)>0 THEN
PRINT "当選確率";LL(F);"%の時、";N;"回ガチャを引いて";F;"種目が";I;"回当たる確率";C(F,I)/M*100;"%"
END IF
NEXT F
IF C(K,I)>0 THEN PRINT
NEXT I
END

　　　　　　　　　　　　実行結果

当たりの種類 3
ガチャを引く回数 100
1の当選確率 (%)(1-99) 1
2の当選確率 (%)(1-99) 2
3の当選確率 (%)(1-99) 3
当選確率 1 %の時、 100 回ガチャを引いて 1 種目が 1 回当たる確率 23.8 %
当選確率 2 %の時、 100 回ガチャを引いて 2 種目が 1 回当たる確率 32.4 %
当選確率 3 %の時、 100 回ガチャを引いて 3 種目が 1 回当たる確率 38 %

当選確率 1 %の時、 100 回ガチャを引いて 1 種目が 2 回当たる確率 4 %
当選確率 2 %の時、 100 回ガチャを引いて 2 種目が 2 回当たる確率 11.8 %
当選確率 3 %の時、 100 回ガチャを引いて 3 種目が 2 回当たる確率 17.2 %

当選確率 2 %の時、 100 回ガチャを引いて 2 種目が 3 回当たる確率 2.6 %
当選確率 3 %の時、 100 回ガチャを引いて 3 種目が 3 回当たる確率 7 %

当選確率 3 %の時、 100 回ガチャを引いて 3 種目が 4 回当たる確率 1.4 %

当選確率 3 %の時、 100 回ガチャを引いて 3 種目が 5 回当たる確率 .6 %
-----------------------------------------------------------------------
これ以下より当たりがK種で内2つ、3つ、4つが当たりの場合のプログラムを挙げているが
当たり判定を微妙に変えているので注意すること

当たりがK種ある場合の内、少なくとも2つが当たる確率を個別に求める



A,B,Cの3種類当たりがあるとして(A>0 を当たり、A=0をハズレとする)

少なくともAかBが当たり A>0 OR B>0
少なくともAとBが当たり A>0 AND B>0
AかBのみが当たり (A>0 OR B>0) AND C=0
AとBのみが当たり A>0 AND B>0 AND C=0



また、IF GACHA(P)=JJ THEN LET HIT(JJ)=HIT(JJ)+1 として
当たる確率を求めているが

IF GACHA(P)<>JJ THEN LET MISS(JJ)=MISS(JJ)+1 と変更すれば
外れの確率を求めることもできる


RANDOMIZE
INPUT PROMPT "当たりの種類":K
DIM GACHA(100*K),HIT(K),EITHER(DEFARRAY2(K,K)),LL(K)
LET M=500 ! 試行回数
LET MODE=3
SELECT CASE MODE
!CASE 1
! INPUT PROMPT "合わせての当選確率 (%)(1-99) ":L
! FOR I=1 TO L
! FOR J=1 TO K
! LET C=C+1
! LET GACHA(C)=J
! NEXT J
! NEXT I
!CASE 2
! INPUT PROMPT "個別の当選確率 (%)(1-99) ":L
! FOR I=1 TO L*K
! FOR J=1 TO K
! LET C=C+1
! LET GACHA(C)=J
! NEXT J
! NEXT I
CASE 3
FOR J=1 TO K
INPUT PROMPT STR$(J)&"の当選確率 (%)(1-99) ":LL(J)
NEXT J
FOR I=1 TO K
FOR J=1 TO LL(I)
LET C=C+1
LET GACHA(C)=I
NEXT J
NEXT I
CASE 4
FOR J=1 TO K
INPUT PROMPT STR$(J)&"の当選確率 (%)(1-99) ":LL(J)
NEXT J
FOR I=1 TO K
FOR J=1 TO LL(I)*K
LET C=C+1
LET GACHA(C)=I
NEXT J
NEXT I
END SELECT
INPUT PROMPT "ガチャを引く回数 ":N
FOR I=1 TO M
FOR J=1 TO 100*K ! シャッフル
SWAP GACHA(J),GACHA(INT(RND*100*K+1))
NEXT J
MAT HIT=ZER
FOR J=1 TO N
LET P=INT(RND*100*K+1) ! ガチャを引く
FOR JJ=1 TO K
IF GACHA(P)=JJ THEN LET HIT(JJ)=HIT(JJ)+1 ! 当たりなら
NEXT JJ
NEXT J
FOR J1=1 TO K-1
FOR J2=J1+1 TO K
IF HIT(J1)>0 OR HIT(J2)>0 THEN
LET EITHER(ARRAY2(J1,J2))=EITHER(ARRAY2(J1,J2))+1
END IF
NEXT J2
NEXT J1
NEXT I
FOR J1=1 TO K-1
FOR J2=J1+1 TO K
PRINT N;"回ガチャを引いて少なくとも";J1;"か";J2;"が当たる確率";EITHER(ARRAY2(J1,J2))/M*100;"%"
NEXT J2
NEXT J1

FUNCTION DEFARRAY2(AMAX,BMAX)
LET MAX1=AMAX
LET MAX2=BMAX
LET DEFARRAY2=MAX2*MAX1
END FUNCTION

FUNCTION ARRAY2(X,Y)
LET X=INT(X)
LET Y=INT(Y)
IF X>MAX1 THEN LET X=MAX1
IF Y>MAX2 THEN LET Y=MAX2
IF X<1 THEN LET X=1
IF Y<1 THEN LET Y=1
LET SUM= MAX1*(Y-1)
LET SUM=SUM+X
LET ARRAY2=SUM
END FUNCTION
END

　　　　　　　　　　実行結果


当たりの種類3
1の当選確率 (%)(1-99) 1
2の当選確率 (%)(1-99) 2
3の当選確率 (%)(1-99) 3
ガチャを引く回数 100
100 回ガチャを引いて少なくとも 1 か 2 が当たる確率 63.2 %
100 回ガチャを引いて少なくとも 1 か 3 が当たる確率 76 %
100 回ガチャを引いて少なくとも 2 か 3 が当たる確率 81.2 %
-----------------------------------------------------------------------
当たりがK種ある場合の内3つのみが当たる確率を個別に求める

RANDOMIZE
DO
INPUT PROMPT "当たりの種類":K ! K>=4
LOOP WHILE K<4
DIM GACHA(100*K),HIT(K),EITHER(DEFARRAY3(K,K,K)),LL(K)
LET M=500 ! 試行回数
LET MODE=4
SELECT CASE MODE
!CASE 1
! INPUT PROMPT "合わせての当選確率 (%)(1-99) ":L
! FOR I=1 TO L
! FOR J=1 TO K
! LET C=C+1
! LET GACHA(C)=J
! NEXT J
! NEXT I
!CASE 2
! INPUT PROMPT "個別の当選確率 (%)(1-99) ":L
! FOR I=1 TO L*K
! FOR J=1 TO K
! LET C=C+1
! LET GACHA(C)=J
! NEXT J
! NEXT I
CASE 3
FOR J=1 TO K
INPUT PROMPT STR$(J)&"の当選確率 (%)(1-99) ":LL(J)
NEXT J
FOR I=1 TO K
FOR J=1 TO LL(I)
LET C=C+1
LET GACHA(C)=I
NEXT J
NEXT I
CASE 4
FOR J=1 TO K
INPUT PROMPT STR$(J)&"の当選確率 (%)(1-99) ":LL(J)
NEXT J
FOR I=1 TO K
FOR J=1 TO LL(I)*K
LET C=C+1
LET GACHA(C)=I
NEXT J
NEXT I
END SELECT
INPUT PROMPT "ガチャを引く回数 ":N
FOR I=1 TO M
FOR J=1 TO 100*K ! シャッフル
SWAP GACHA(J),GACHA(INT(RND*100*K+1))
NEXT J
MAT HIT=ZER
FOR J=1 TO N
LET P=INT(RND*100*K+1) ! ガチャを引く
FOR JJ=1 TO K
IF GACHA(P)=JJ THEN LET HIT(JJ)=HIT(JJ)+1 ! 当たりなら
NEXT JJ
NEXT J
FOR J1=1 TO K-2
FOR J2=J1+1 TO K-1
FOR J3=J2+1 TO K
IF HIT(J1)>0 AND HIT(J2)>0 AND HIT(J3)>0 THEN
LET FLG=0
FOR JJ=1 TO K
IF J1<>JJ AND J2<>JJ AND J3<>JJ AND HIT(JJ)<>0 THEN LET FLG=1
NEXT JJ
IF FLG=0 THEN LET EITHER(ARRAY3(J1,J2,J3))=EITHER(ARRAY3(J1,J2,J3))+1
END IF
NEXT J3
NEXT J2
NEXT J1
NEXT I
FOR J1=1 TO K-2
FOR J2=J1+1 TO K-1
FOR J3=J2+1 TO K
PRINT N;"回ガチャを引いて";J1;"と";J2;"と";J3;"のみが当たる確率";EITHER(ARRAY3(J1,J2,J3))/M*100;"%"
NEXT J3
NEXT J2
NEXT J1

FUNCTION DEFARRAY3(AMAX,BMAX,CMAX)
LET MAX1=AMAX
LET MAX2=BMAX
LET MAX3=CMAX
LET DEFARRAY3=MAX3*MAX2*MAX1
END FUNCTION

FUNCTION ARRAY3(X,Y,Z)
LET X=INT(X)
LET Y=INT(Y)
LET Z=INT(Z)
IF X>MAX1 THEN LET X=MAX1
IF Y>MAX2 THEN LET Y=MAX2
IF Z>MAX3 THEN LET Z=MAX3
IF X<1 THEN LET X=1
IF Y<1 THEN LET Y=1
IF Z<1 THEN LET Z=1
LET SUM= MAX2*MAX1*(Z-1)
LET SUM=SUM+MAX1*(Y-1)
LET SUM=SUM+X
LET ARRAY3=SUM
END FUNCTION
END

　　　　　　　　　実行結果


当たりの種類6
1の当選確率 (%)(1-99) 1
2の当選確率 (%)(1-99) 2
3の当選確率 (%)(1-99) 3
4の当選確率 (%)(1-99) 4
5の当選確率 (%)(1-99) 5
6の当選確率 (%)(1-99) 6
ガチャを引く回数 100
100 回ガチャを引いて 1 と 2 と 3 のみが当たる確率 0 %
100 回ガチャを引いて 1 と 2 と 4 のみが当たる確率 0 %
100 回ガチャを引いて 1 と 2 と 5 のみが当たる確率 0 %
100 回ガチャを引いて 1 と 2 と 6 のみが当たる確率 0 %
100 回ガチャを引いて 1 と 3 と 4 のみが当たる確率 0 %
100 回ガチャを引いて 1 と 3 と 5 のみが当たる確率 0 %
100 回ガチャを引いて 1 と 3 と 6 のみが当たる確率 0 %
100 回ガチャを引いて 1 と 4 と 5 のみが当たる確率 0 %
100 回ガチャを引いて 1 と 4 と 6 のみが当たる確率 0 %
100 回ガチャを引いて 1 と 5 と 6 のみが当たる確率 0 %
100 回ガチャを引いて 2 と 3 と 4 のみが当たる確率 0 %
100 回ガチャを引いて 2 と 3 と 5 のみが当たる確率 0 %
100 回ガチャを引いて 2 と 3 と 6 のみが当たる確率 0 %
100 回ガチャを引いて 2 と 4 と 5 のみが当たる確率 0 %
100 回ガチャを引いて 2 と 4 と 6 のみが当たる確率 0 %
100 回ガチャを引いて 2 と 5 と 6 のみが当たる確率 0 %
100 回ガチャを引いて 3 と 4 と 5 のみが当たる確率 .2 %
100 回ガチャを引いて 3 と 4 と 6 のみが当たる確率 0 %
100 回ガチャを引いて 3 と 5 と 6 のみが当たる確率 0 %
100 回ガチャを引いて 4 と 5 と 6 のみが当たる確率 .2 %
-----------------------------------------------------------------------
当たりがK種ある場合の内4つが当たる確率を個別に求める


RANDOMIZE
DO
INPUT PROMPT "当たりの種類":K ! K>=5
LOOP WHILE K<5
DIM GACHA(100*K),HIT(K),EITHER(DEFARRAY4(K,K,K,K)),LL(K)
LET M=500 ! 試行回数
LET MODE=4
SELECT CASE MODE
!CASE 1
! INPUT PROMPT "合わせての当選確率 (%)(1-99) ":L
! FOR I=1 TO L
! FOR J=1 TO K
! LET C=C+1
! LET GACHA(C)=J
! NEXT J
! NEXT I
!CASE 2
! INPUT PROMPT "個別の当選確率 (%)(1-99) ":L
! FOR I=1 TO L*K
! FOR J=1 TO K
! LET C=C+1
! LET GACHA(C)=J
! NEXT J
! NEXT I
CASE 3
FOR J=1 TO K
INPUT PROMPT STR$(J)&"の当選確率 (%)(1-99) ":LL(J)
NEXT J
FOR I=1 TO K
FOR J=1 TO LL(I)
LET C=C+1
LET GACHA(C)=I
NEXT J
NEXT I
CASE 4
FOR J=1 TO K
INPUT PROMPT STR$(J)&"の当選確率 (%)(1-99) ":LL(J)
NEXT J
FOR I=1 TO K
FOR J=1 TO LL(I)*K
LET C=C+1
LET GACHA(C)=I
NEXT J
NEXT I
END SELECT
INPUT PROMPT "ガチャを引く回数 ":N
FOR I=1 TO M
FOR J=1 TO 100*K ! シャッフル
SWAP GACHA(J),GACHA(INT(RND*100*K+1))
NEXT J
MAT HIT=ZER
FOR J=1 TO N
LET P=INT(RND*100*K+1) ! ガチャを引く
FOR JJ=1 TO K
IF GACHA(P)=JJ THEN LET HIT(JJ)=HIT(JJ)+1 ! 当たりなら
NEXT JJ
NEXT J
FOR J1=1 TO K-3
FOR J2=J1+1 TO K-2
FOR J3=J2+1 TO K-1
FOR J4=J3+1 TO K
IF HIT(J1)>0 OR HIT(J2)>0 OR HIT(J3)>0 OR HIT(J4)>0 THEN
LET FLG=0
FOR JJ=1 TO K
IF J1<>JJ AND J2<>JJ AND J3<>JJ AND J4<>JJ AND HIT(JJ)<>0 THEN LET FLG=1
NEXT JJ
IF FLG=0 THEN LET EITHER(ARRAY4(J1,J2,J3,J4))=EITHER(ARRAY4(J1,J2,J3,J4))+1
END IF
NEXT J4
NEXT J3
NEXT J2
NEXT J1
NEXT I
FOR J1=1 TO K-3
FOR J2=J1+1 TO K-2
FOR J3=J2+1 TO K-1
FOR J4=J3+1 TO K
PRINT N;"回ガチャを引いて";J1;"か";J2;"か";J3;"か";J4;"のみが当たる確率";EITHER(ARRAY4(J1,J2,J3,J4))/M*100;"%"
NEXT J4
NEXT J3
NEXT J2
NEXT J1

FUNCTION DEFARRAY4(AMAX,BMAX,CMAX,DMAX)
LET MAX1=AMAX
LET MAX2=BMAX
LET MAX3=CMAX
LET MAX4=DMAX
LET DEFARRAY4=MAX4*MAX3*MAX2*MAX1
END FUNCTION

FUNCTION ARRAY4(X,Y,Z,W)
LET X=INT(X)
LET Y=INT(Y)
LET Z=INT(Z)
LET W=INT(W)
IF X>MAX1 THEN LET X=MAX1
IF Y>MAX2 THEN LET Y=MAX2
IF Z>MAX3 THEN LET Z=MAX3
IF W>MAX4 THEN LET W=MAX4
IF X<1 THEN LET X=1
IF Y<1 THEN LET Y=1
IF Z<1 THEN LET Z=1
IF W<1 THEN LET W=1
LET SUM= MAX3*MAX2*MAX1*(W-1)
LET SUM=SUM+MAX2*MAX1*(Z-1)
LET SUM=SUM+MAX1*(Y-1)
LET SUM=SUM+X
LET ARRAY4=SUM
END FUNCTION
END

　　　　　　　　　実行結果

当たりの種類5
1の当選確率 (%)(1-99) 1
2の当選確率 (%)(1-99) 2
3の当選確率 (%)(1-99) 3
4の当選確率 (%)(1-99) 4
5の当選確率 (%)(1-99) 5
ガチャを引く回数 100
100 回ガチャを引いて 1 か 2 か 3 か 4 のみが当たる確率 .6 %
100 回ガチャを引いて 1 か 2 か 3 か 5 のみが当たる確率 1.4 %
100 回ガチャを引いて 1 か 2 か 4 か 5 のみが当たる確率 5.2 %
100 回ガチャを引いて 1 か 3 か 4 か 5 のみが当たる確率 12.6 %
100 回ガチャを引いて 2 か 3 か 4 か 5 のみが当たる確率 33 %

Re: 実行ファイルの作り方 - SHIRAISHI Kazuo

2026/05/02 (Sat) 07:46:01

2進モード，複素数モードで動くプログラムは，たいていの場合，BASIC Accelerator で実行ファイル化できます。
https://decimalbasic.web.fc2.com/BASICAccJa.htm
実行ファイルは，outoutフォルダにNoName.exeとして生成されるので，適宜，名称を変更してください。
特別なDLLに依存せずに動作します。

Re: 実行ファイルの作り方 - 永野護

2026/05/02 (Sat) 08:20:14

ありがとうございました。

Re: ver7881がダウンロードできません - SHIRAISHI Kazuo

2026/04/22 (Wed) 11:42:08

ご報告ありがとうございました。
BASIC7881.zipはBASICW95のアーカイブになっていました。お手数をお掛けしますが，再度，ダウンロードをお願いします。

Re: じゃんけん - しばっち

2026/04/12 (Sun) 07:46:03

じゃんけんをして一人だけが勝ち抜けするまでの回数
但し、途中で勝っても抜けず一人だけが勝つ迄の回数


RANDOMIZE
INPUT PROMPT "人数=":N
DIM A(N),J$(0 TO 2),P(0 TO 2)
LET J$(0)="グー"
LET J$(1)="チョキ"
LET J$(2)="パー"
DO
LET COUNT=COUNT+1
MAT A=ZER
MAT P=ZER
! PRINT COUNT;"回"
FOR I=1 TO N
LET A(I)=INT(RND*3)
! PRINT J$(A(I)),
NEXT I
! PRINT
FOR I=1 TO N
LET P(A(I))=P(A(I))+1
NEXT I
IF P(0)=1 AND P(1)>0 AND P(2)=0 THEN EXIT DO ! グーの一人勝ち
IF P(0)>0 AND P(1)=0 AND P(2)=1 THEN EXIT DO ! パーの一人勝ち
IF P(0)=0 AND P(1)=1 AND P(2)>0 THEN EXIT DO ! チョキの一人勝ち
LOOP
PRINT COUNT;"回のジャンケンで一人勝ちになりました"
END
-------------------------------------------------------------------
20回試行して一人勝ちする迄の平均回数を求めてみた。


RANDOMIZE
DIM A(15),P(0 TO 2)
LET L=20
FOR NN=2 TO 15 ! 人数
LET SUM=0
FOR K=1 TO L ! 試行回数
LET N=NN
LET COUNT=0
DO
LET COUNT=COUNT+1
MAT P=ZER
FOR I=1 TO N
LET A(I)=INT(RND*3)
NEXT I
FOR I=1 TO N
LET P(A(I))=P(A(I))+1
NEXT I
IF P(0)=1 AND P(1)>0 AND P(2)=0 THEN EXIT DO ! グーの一人勝ち
IF P(0)>0 AND P(1)=0 AND P(2)=1 THEN EXIT DO ! パーの一人勝ち
IF P(0)=0 AND P(1)=1 AND P(2)>0 THEN EXIT DO ! チョキの一人勝ち
LOOP
LET SUM=SUM+COUNT
NEXT K
PRINT NN;"人の時";SUM/L;"回のジャンケンで一人勝ちになりました"
NEXT NN
END

　　　　　　　　　実行結果


2 人の時 1.7 回のジャンケンで一人勝ちになりました
3 人の時 2.4 回のジャンケンで一人勝ちになりました
4 人の時 3.95 回のジャンケンで一人勝ちになりました
5 人の時 14.65 回のジャンケンで一人勝ちになりました
6 人の時 43.25 回のジャンケンで一人勝ちになりました
7 人の時 107.45 回のジャンケンで一人勝ちになりました
8 人の時 217.75 回のジャンケンで一人勝ちになりました
9 人の時 740.9 回のジャンケンで一人勝ちになりました
10 人の時 2258.8 回のジャンケンで一人勝ちになりました
11 人の時 3811.85 回のジャンケンで一人勝ちになりました
12 人の時 22032.85 回のジャンケンで一人勝ちになりました
13 人の時 40635.15 回のジャンケンで一人勝ちになりました
14 人の時 85129.2 回のジャンケンで一人勝ちになりました
15 人の時 242720.1 回のジャンケンで一人勝ちになりました

Re: 同じ誕生日の人 - しばっち

2026/04/12 (Sun) 07:43:12

1クラス内に自分と同じ誕生日の人がいる確率


RANDOMIZE
LET M=300 ! 試行回数
INPUT PROMPT "１クラスの人数=":N
DIM A(N)
LET P=INT(RND*365+1) ! 自分の誕生日
FOR K=1 TO M
FOR I=1 TO N
LET A(I)=INT(RND*365+1) ! クラスの誕生日
NEXT I
FOR I=1 TO N
IF A(I)=P THEN
LET C=C+1 ! 自分と同じならカウント
EXIT FOR
END IF
NEXT I
NEXT K
PRINT "自分と同じ誕生日の人がいる確率";C/M*100;"%"
END
--------------------------------------------------------------
300回試行して平均求めてみた


RANDOMIZE
LET M=300 ! 試行回数
DIM A(50)
FOR N=2 TO 50 ! 人数
LET C=0
LET P=INT(RND*365+1) ! 自分の誕生日
FOR K=1 TO M
FOR I=1 TO N
LET A(I)=INT(RND*365+1) ! クラスの誕生日
NEXT I
FOR I=1 TO N
IF A(I)=P THEN
LET C=C+1 ! 自分と同じならカウント
EXIT FOR
END IF
NEXT I
NEXT K
PRINT N;"人のクラスの中に自分と同じ誕生の日の人がいる確率";C/M*100;"%"
NEXT N
END


　　　　　　　　　　実行結果


2 人のクラスの中に自分と同じ誕生の日の人がいる確率 .333333333333333 %
3 人のクラスの中に自分と同じ誕生の日の人がいる確率 .666666666666667 %
4 人のクラスの中に自分と同じ誕生の日の人がいる確率 .333333333333333 %
5 人のクラスの中に自分と同じ誕生の日の人がいる確率 1.33333333333333 %
6 人のクラスの中に自分と同じ誕生の日の人がいる確率 2.33333333333333 %
7 人のクラスの中に自分と同じ誕生の日の人がいる確率 2.66666666666667 %
8 人のクラスの中に自分と同じ誕生の日の人がいる確率 2.66666666666667 %
9 人のクラスの中に自分と同じ誕生の日の人がいる確率 2.33333333333333 %
10 人のクラスの中に自分と同じ誕生の日の人がいる確率 6 %
11 人のクラスの中に自分と同じ誕生の日の人がいる確率 2.66666666666667 %
12 人のクラスの中に自分と同じ誕生の日の人がいる確率 3 %
13 人のクラスの中に自分と同じ誕生の日の人がいる確率 3.66666666666667 %
14 人のクラスの中に自分と同じ誕生の日の人がいる確率 5 %
15 人のクラスの中に自分と同じ誕生の日の人がいる確率 3.66666666666667 %
16 人のクラスの中に自分と同じ誕生の日の人がいる確率 3.33333333333333 %
17 人のクラスの中に自分と同じ誕生の日の人がいる確率 6 %
18 人のクラスの中に自分と同じ誕生の日の人がいる確率 3.66666666666667 %
19 人のクラスの中に自分と同じ誕生の日の人がいる確率 5.33333333333333 %
20 人のクラスの中に自分と同じ誕生の日の人がいる確率 6.33333333333333 %
21 人のクラスの中に自分と同じ誕生の日の人がいる確率 6.33333333333333 %
22 人のクラスの中に自分と同じ誕生の日の人がいる確率 6.66666666666667 %
23 人のクラスの中に自分と同じ誕生の日の人がいる確率 5.33333333333333 %
24 人のクラスの中に自分と同じ誕生の日の人がいる確率 6.33333333333333 %
25 人のクラスの中に自分と同じ誕生の日の人がいる確率 6.66666666666667 %
26 人のクラスの中に自分と同じ誕生の日の人がいる確率 6.66666666666667 %
27 人のクラスの中に自分と同じ誕生の日の人がいる確率 3.66666666666667 %
28 人のクラスの中に自分と同じ誕生の日の人がいる確率 7.33333333333333 %
29 人のクラスの中に自分と同じ誕生の日の人がいる確率 9.33333333333333 %
30 人のクラスの中に自分と同じ誕生の日の人がいる確率 5.33333333333333 %
31 人のクラスの中に自分と同じ誕生の日の人がいる確率 9 %
32 人のクラスの中に自分と同じ誕生の日の人がいる確率 7.66666666666667 %
33 人のクラスの中に自分と同じ誕生の日の人がいる確率 8.66666666666667 %
34 人のクラスの中に自分と同じ誕生の日の人がいる確率 8.66666666666667 %
35 人のクラスの中に自分と同じ誕生の日の人がいる確率 6.33333333333333 %
36 人のクラスの中に自分と同じ誕生の日の人がいる確率 8.66666666666667 %
37 人のクラスの中に自分と同じ誕生の日の人がいる確率 10.6666666666667 %
38 人のクラスの中に自分と同じ誕生の日の人がいる確率 7.66666666666667 %
39 人のクラスの中に自分と同じ誕生の日の人がいる確率 11.6666666666667 %
40 人のクラスの中に自分と同じ誕生の日の人がいる確率 7.66666666666667 %
41 人のクラスの中に自分と同じ誕生の日の人がいる確率 10.6666666666667 %
42 人のクラスの中に自分と同じ誕生の日の人がいる確率 10 %
43 人のクラスの中に自分と同じ誕生の日の人がいる確率 13.3333333333333 %
44 人のクラスの中に自分と同じ誕生の日の人がいる確率 12.3333333333333 %
45 人のクラスの中に自分と同じ誕生の日の人がいる確率 7.66666666666667 %
46 人のクラスの中に自分と同じ誕生の日の人がいる確率 11.6666666666667 %
47 人のクラスの中に自分と同じ誕生の日の人がいる確率 12.3333333333333 %
48 人のクラスの中に自分と同じ誕生の日の人がいる確率 9.66666666666667 %
49 人のクラスの中に自分と同じ誕生の日の人がいる確率 10.6666666666667 %
50 人のクラスの中に自分と同じ誕生の日の人がいる確率 13 %
--------------------------------------------------------------
確率から人数求めてみた


RANDOMIZE
LET F=300 ! 試行回数
DIM A(10000)
DO
INPUT PROMPT "自分と同じ誕生日の人がいる確率 (%)":Z
LOOP UNTIL Z<100 AND Z>0
LET Z=Z/100
LET L=1
LET R=10000 ! 人数
LET P=INT(RND*365+1)
DO
LET N=INT((L+R)/2)
LET C=0
FOR K=1 TO F
FOR I=1 TO N
LET A(I)=INT(RND*365+1) ! クラスの誕生日
NEXT I
FOR I=1 TO N
IF A(I)=P THEN
LET C=C+1
EXIT FOR
END IF
NEXT I
NEXT K
IF C/F<Z THEN LET L=N ELSE LET R=N
LOOP UNTIL L=R
PRINT L;"人の時、自分と同じ誕生日の人がいる確率";Z*100;"%になります"
END

Re: 十進BASIC版：素数計数関数 pi(x)正確値計算プログラム - しばっち

2026/04/05 (Sun) 07:54:18

私も素数個数関数をルジャンドル法で作ってみた。
とりあえずやっと動くようになりました。

まだ改良の余地はあるかとは思いますが
参考程度にはなるかと思います。


https://mathworld.wolfram.com/LegendresFormula.html
https://www.hs.chuo-u.ac.jp/contents/wp-content/themes/chu-fu/pdf/bulletin/issue30/issue30_pdf10.pdf
https://rosettacode.org/wiki/Legendre_prime_counting_function


OPTION BASE 0
LET N=1000000
PUBLIC NUMERIC PICACHE(10000),PRIME(10000),CACHE(10000,100),COUNT
MAT CACHE=(-1)*CON
CALL PRIMESIEVE(INT(SQR(N)))
PRINT π(N)
END

EXTERNAL SUB PRIMESIEVE(N)
MAT PICACHE=CON
LET PICACHE(1)=0
LET COUNT=0
FOR I=2 TO N
IF PICACHE(I)=1 THEN
LET COUNT=COUNT+1
LET PRIME(COUNT)=I
FOR J=I*I TO N STEP I
LET PICACHE(J)=0
NEXT J
END IF
NEXT I
FOR I=1 TO N-1
LET PICACHE(I+1)=PICACHE(I+1)+PICACHE(I)
NEXT I
END SUB

EXTERNAL FUNCTION φ(X,A)
IF X=0 THEN
LET φ=0
EXIT FUNCTION
END IF
IF X=1 THEN
LET φ=1
EXIT FUNCTION
END IF
IF A=1 THEN
LET φ=IP(X/2)
LET CACHE(X,1)=IP(X/2)
EXIT FUNCTION
END IF
IF A=2 THEN
LET φ=IP(X/3)
LET CACHE(X,2)=IP(X/3)
EXIT FUNCTION
END IF
IF A=3 THEN ! COMB(3,1)+COMB(3,2)+COMB(3,3) 2,3,5
LET S=X-IP(X/2)-IP(X/3)-IP(X/5)+IP(X/(2*3))+IP(X/(2*5))+IP(X/(3*5))-IP(X/(2*3*5))
LET φ=S
IF X<10000 THEN LET CACHE(X,3)=S
EXIT FUNCTION
END IF
IF A=4 THEN ! COMB(4,1)+COMB(4,2)+COMB(4,3)+COMB(4,4) 2,3,5,7
LET S=X-IP(X/2)-IP(X/3)-IP(X/5)-IP(X/7)+IP(X/(2*3))+IP(X/(2*5))+IP(X/(2*7))+IP(X/(3*5))+IP(X/(3*7))+IP(X/(5*7))-IP(X/(2*3*5))-IP(X/(2*3*7))-IP(X/(2*5*7))-IP(X/(3*5*7))+IP(X/(2*3*5*7))
LET φ=S
IF X<10000 THEN LET CACHE(X,4)=S
EXIT FUNCTION
END IF
IF A=5 THEN ! COMB(5,1)+COMB(5,2)+COMB(5,3)+COMB(5,4)+COMB(5,5) 2,3,5,7,11
LET S=X-IP(X/2)-IP(X/3)-IP(X/5)-IP(X/7)-IP(X/11)+IP(X/(2*3))+IP(X/(2*5))+IP(X/(2*7))+IP(X/(2*11))+IP(X/(3*5))+IP(X/(3*7))+IP(X/(3*11))+IP(X/(5*7))+IP(X/(5*11))+IP(X/(7*11))
LET S=S-IP(X/(2*3*5))-IP(X/(2*3*7))-IP(X/(2*3*11))-IP(X/(2*5*7))-IP(X/(2*5*11))-IP(X/(2*7*11))-IP(X/(3*5*7))-IP(X/(3*5*11))-IP(X/(3*7*11))-IP(X/(5*7*11))
LET S=S+IP(X/(2*3*5*7))+IP(X/(2*3*5*11))+IP(X/(2*3*7*11))+IP(X/(2*5*7*11))+IP(X/(3*5*7*11))-IP(X/(2*3*5*7*11))
LET φ=S
IF X<10000 THEN LET CACHE(X,5)=S
EXIT FUNCTION
END IF
IF X<10000 AND A<100 AND CACHE(X,A)<>-1 THEN
LET S=CACHE(X,A)
ELSE
LET S=φ(X,A-1)
IF X<10000 AND A<100 THEN LET CACHE(X,A)=S
END IF
LET XX=IP(X/PRIME(A))
IF XX<10000 AND A<100 AND CACHE(XX,A)<>-1 THEN
LET T=CACHE(XX,A)
ELSE
LET T=φ(XX,A-1)
IF XX<10000 AND A<100 THEN LET CACHE(XX,A)=T
END IF
LET φ=S-T
END FUNCTION

EXTERNAL FUNCTION π(X)
IF X<=PRIME(COUNT) THEN
LET π=PICACHE(X)
EXIT FUNCTION
END IF
!LET T=π(INT(SQR(X)))
LET T=COUNT
LET π=T+φ(X,T)-1
END FUNCTION

Re: 要望 - SHIRAISHI Kazuo

2026/03/15 (Sun) 13:16:47

X^N = EXP(N*LOG(X)) と定義してしまえば可能ですが，一般には負数の非整数乗は未定義とします。
たとえば，
https://detail.chiebukuro.yahoo.co.jp/qa/question_detail/q14250831054
要するに，値が１つに定まらないからで，0÷0を未定義にするのと同様の理屈です。

Re: 十進BASIC7.8.8の文字化け - SHIRAISHI Kazuo

2026/03/07 (Sat) 10:17:04

7.8.8以前のバージョンでは問題ないのでしょうか。
ver. 7.8.7.7と7.8.7.6が　https://decimalbasic.web.fc2.com/basicw32.htm　から入手可能です。

なお，Windows 11 Home 25H2 26200.7840　でテストする範囲では正常です。

Re: 十進BASIC7.8.8の文字化け - knoike

2026/03/07 (Sat) 12:39:21

knoike です。ここでは初めて書き込みます。

設定 → 時刻と言語 → 言語と地域
で、
[Windows の表示言語] の右端の v をクリックして展開される項目の
[ベータ版: 世界中の言語に対応するために Unicode UTF-8 を使用します]
が
[オン] になっていると化けます。
デフォルトは [オフ] で、[オフ] なら化けません。正常に表示されます。

「OSのローケルや表示を日本語優先にすることである程度は直りましたが、」という書き込みから、
言語パックをいくつかインストールしてあるのだと想像していますが、
たいていのアプリはここの設定が [オフ] のままで正常に表示できるような多言語対応がなされていると思いますので、
特別な事情がない限りは [オフ] にしておくとよいのではないでしょうか。

Re: 十進BASIC7.8.8の文字化け - た

2026/03/07 (Sat) 20:44:35

ありがとうございます。

>SHIRAISHI Kazuo さん
旧バージョンを試してみましたが、結果は同じでした。

>knoikeさん
ご指摘の部分がONでしたので、OFFにしてみましたが、改善は見られませんでした。

みなさん、ご指摘の通り、私のwindows 11の環境が怪しそうなので、まず、OS周りで不具合がないか、一度確認してみます。

ありがとうございました。

Re: 十進BASIC7.8.8の文字化け - た

2026/03/07 (Sat) 20:53:31

今、試しで、「十進BASIC ver. 8」を使ってみたら、全く文字化けがありません。


文字コードがutf-8に対応した、ということですが、Version 7の中の別の文字コードが化けている可能性があるのでは、と推定しています。

現時点で、Version 8であれば文字化けは発生しませんので、当面はこちらを使いたいと思います。

ありがとうございました。

Re: 十進BASIC7.8.8の文字化け - SHIRAISHI Kazuo

2026/03/08 (Sun) 07:55:50

設定 → 時刻と言語 → 言語と地域 [ベータ版: 世界中の言語に対応するために Unicode UTF-8 を使用します] を試してみました。
投稿された通りの結果になります。
Ver.7は内部でshift-JISを採用しています。また，「機能語挿入」の内容部分は　BASIC.kwFとBASIC.kwSから読み込んでいます。
BASIC.kwFとBASIC.kwSはshift-JISで書かれているので文字化けして当然ですが，それ以外の部分はβ版で未対応なのだろうと思います。
なお，Ver.7 日本語版の翻訳系は，文字コードがShift-JISであることを前提にしています。UTF-8の文字を含むプログラムは，おそらく正常に実行できません。

Re: 十進BASIC7.8.8の文字化け - SHIRAISHI Kazuo

2026/03/08 (Sun) 11:59:40

十進BASIC Ver.8は，内部文字コードがUTF-8ですが，Windows APIは16ビットwide文字対応版を使うので，Windowsのlocaleの影響を受けません。
BASICAccなども同様です。

Re: 十進BASIC7.8.8の文字化け - SHIRAISHI Kazuo

2026/03/09 (Mon) 09:21:33

十進BASIC ver.8（日本語版）はWindowsのコードページ設定に関係なく日本語で動作します。
Ver.8でもWindows APIの呼び出しが使えます。
https://decimalbasic.web.fc2.com/ExtDLL.htm
WindowsAPIのMessageBoxAは，ANSI文字でメッセージボックスを表示します。
WindowsのANSI文字というのは，Windowsのコードページごとに定まる文字のことで，日本語版WindowsではShift-JISのことです。

DECLARE EXTERNAL FUNCTION MesBox
LET n=MesBox(0,"日本語","BASIC",3)
PRINT n
END
EXTERNAL FUNCTION MesBox(owner,text$,caption$,flag)
ASSIGN "user32.dll","MessageBoxA"
END FUNCTION

をVer.8で実行すると文字化けします(Ver.8で"日本語"がUTF-8文字なので)。
"日本語" のところを ANSI$("日本語") に変えると正常に機能します。
ANSI$関数は，コードページの指定にしたがってUTF-8文字列をANSI文字列に変える関数です。
MessageBoxAの説明は
https://learn.microsoft.com/ja-jp/windows/win32/api/winuser/nf-winuser-messageboxa
にあります。

「世界中の言語に対応するために Unicode UTF-8 を使用します」の意味は，ANSI文字としてUTF-8を採用することのように思われます。
Ver.8はWindowsの言語設定を変えても動作するので，テストしてみていただけないでしょうか。

Re: 十進BASIC7.8.8の文字化け - た

2026/03/14 (Sat) 19:05:46

返信が遅くなり、申し訳ございません。

承知しました。
WindowsAPIの利用含め、Version.8を一通り試して、何か気が付いた点がございましたら、また報告させていただきます。

お時間いただきまして、ありがとうございました。

よろしくお願いいたします。

Re: arm版Windows11 - SHIRAISHI Kazuo

2026/02/16 (Mon) 08:45:54

テストプログラムです。
10 OPTION ARITHMETIC NATIVE
20 LET x=0
30 WHEN EXCEPTION IN
40 PRINT 1/x
50 USE
60 PRINT "正常"
70 END WHEN
80 END
Ver.7 をarm版Winows11で実行すると，このテストプログラムでフリーズします。
Ver.8 だとこのテストは通過しますが，複素数モードのサンプルプログラムMU_COSH.BASやMANDELBL.BASを実行するとすべて黒で塗りつぶされてしまいます。
テスト環境は，MAC (Apple Silicon) 上のVrtualBox上のWindows11 Canary です。
実機だと異なる結果になるかもしれません。実機 Windows on arm での動作報告を期待（歓迎）します。

Re: arm版Windows11 - SHIRAISHI Kazuo

2026/02/16 (Mon) 20:24:22

MAC (Apple Silicon) 上のVrtualBox上のWindows11 Canaryにおいて，十進BASIC Ver.8（Windows版)で
10 OPTION ARITHMETIC NATIVE
20 LET x=0
30 WHEN EXCEPTION IN
40 PRINT 1/x
50 USE
60 PRINT EXTYPE
70 END WHEN
80 END
を実行すると，結果が1002になります。3001が得られるのが正しいので，Ver.8に対しても，FPU例外のエミュレーションに不具合があるように思います。
ただし，Apple SiliconとWindows機のarmとでは浮動小数点例外の扱いが異なることが原因の可能性もあります。実機版 Windows on arm での動作報告を期待します。

Re: 動作報告です。 - しばっち

2026/02/15 (Sun) 08:03:13

私も作ってみました。
9桁区切りで計算させ無駄なループを省き高速化してみました。

2進モードで実行してください。

OPTION BASE 0
DIM A(4000)
LET A(0)=1
LET T=TIME
FOR I=2 TO 10000
MAT A=(I)*A
FOR J=0 TO N
IF A(J)>=1000000000 THEN
LET R=INT(A(J)/1000000000)
LET A(J)=A(J)-R*1000000000
LET A(J+1)=A(J+1)+R
IF J+1>N THEN LET N=N+1
END IF
NEXT J
NEXT I
PRINT STR$(A(N));
FOR J=N-1 TO 0 STEP-1
PRINT USING "%%%%%%%%%":A(J);
NEXT J
PRINT
PRINT TIME-T;"秒"
END


階乗の計算には1,2,3,4と順に掛けていけばいいが
Nの階乗は1～N/2までを掛けたものとN/2+1～Nまでを掛けたものとを掛けたものと定義することもできます。(2分割)
更に、1～N/3までを掛けたものとN/3+1～2*N/3までを掛けたものと2*N/3+1～Nまでを掛けたものとを掛けた合わせたものと定義することもできます。(3分割)


また100!を下記のように書いてみた。

100!=(1*3*5*7*9*11*13*15*17*19*21*23*25*27*29*31*33*35*37*39*41*43*45*47*49*51*53*55*57*59*61*63*65*67*69*71*73*75*77*79*81*83*85*87*89*91*93*95*97*99)*
2^25*(1*3*5*7*9*11*13*15*17*19*21*23*25*27*29*31*33*35*37*39*41*43*45*47*49)*
4^13*(1*3*5*7*9*11*13*15*17*19*21*23*25)*
8^6 *(1*3*5*7*9*11)*
16^3 *(1*3*5)*
32^2 *(1*3)*
64^1 *(1)

ペキ乗でまとめていけばもう少し最適化できそうですが、すると素因数分解して求めるのと変わらなくなるのでこのままにしておく。



もし、プログラミングを続けていくなら、ぜひ自作(自前)のライブラリーの
作成を勧めます。ライブラリーがあれば毎回1から作る必要がなくなり
2や3から始められます。
ライブラリーが充実していけば7や8からだって始められます。

他言語ならば他の誰か(又は団体や企業等)が作ったライブラリーをダウンロードしてそれを利用するだけですが
十進BASICにもこの掲示板に豊富なプログラム(ライブラリー)がありますので
過去ログを読みなおすこともお薦めします。

まずプログラムの仕様設計を決めることが必要です。
それにはルーチンの呼び出し方を決めルーチン名や使用する変数や個数、配列変数のサイズを決めます。

そして実際にサブルーチン、関数を作成していきます。

多倍長計算でいえばまず足し算、引き算ルーチンを作ります。
表示ルーチン等も併せて作ります。
次に掛け算、割り算ルーチンを作ります。

ルーチンのバージョンアップ(機能拡張)も行っていきます。
負数に対応させるとか、多倍長整数と多倍長小数の両方に対応させる等です。
バージョンアップに伴い仕様設計も変更していきます。

但し、バージョンアップ(更新)していく場合は、バージョンの新しいものと古いものを使用したルーチンが混在することになりますので
更新作業は慎重に行わないとプログラムが動かなくなったりします。

必要なら検索したり管理する作業用のツールプログラムも作ったりします。


四則計算ができたら次は関数ルーチンです。SIN,COS,SQR,EXP,LOG等を作っていきます。
関数ルーチンを基に円周率πなどを求めたり1/πやSQR(π),π^πなどのその派生する定数値を求めたり
するのもいいかと思います。

また、ルーチンを一般化していくのも大事かと思います。

2乗、3乗ルーチンができたらN乗ルーチンを
平方根、立法根ルーチンができたらN乗根を求めるルーチンを作っていきます。

2次、3次とできればN次に、2階、3階とくればN階に対応させていく。
共通項を見つけ出し変数化していくことは大変な作業かもしれませんが
一般化するのは利用価値がより高いからです。

N乗の次はNを整数から実数へと拡張していきます。
実数の次は複素数に対応させたりします。(複素数の次は4元数か!?)

その為にはアルゴリズムの策定や研究が必要になります。
ネットだけでなく書籍(専門書)も参考にしたりします。

テストプログラムも作り動作検証や研究したりも必要になります。

うまく動かなかったり未完成未完動のプログラムも増えていきますが
消さずに残しておきます。後で完成できるかもしれません。

その次のステージは高速化や高機能多機能化(2変数や多変数対応等)などより高度なプログラムへと移行していきます。

高速化は高度な知識、経験等を必要とする高等テクニックです。
単なる思い付きでできるほど簡単な作業ではありません。

その為には日頃からの情報収集や研究などが大切です。



以上は私事の経験則ですがお役に立てれば幸いに思います。

Re: 動作報告です。 - Tarosa

2026/02/15 (Sun) 23:10:14

しばっちさま毎度お世話になります。

!Google　AI例
!10000階乗/123456（十進BASIC）
OPTION BASE 0
DIM A(4000), Q(4000)
LET A(0)=1
LET N=0
LET RADIX = 1000000000

! --- 階乗計算 (A = 10000!) ---
FOR I=2 TO 10000
MAT A=(I)*A
FOR J=0 TO N
IF A(J)>=RADIX THEN
LET R=INT(A(J)/RADIX)
LET A(J)=A(J)-R*RADIX
LET A(J+1)=A(J+1)+R
IF J+1>N THEN LET N=N+1
END IF
NEXT J
NEXT I

! --- 多倍長除算 (Q = A / B) ---
LET B = 123456 ! 除数（割る数）
LET rrem = 0 ! 余り
LET QN = 0 ! 商の最高次
FOR J = N TO 0 STEP -1
! 筆算の考え方：前の桁の余りを10億倍して現在の桁に足す
LET current = rrem * RADIX + A(J)
LET Q(J) = INT(current / B)
LET rrem = MOD(current, B)
IF Q(J) > 0 AND J > QN THEN LET QN = J
NEXT J

! --- 結果表示 ---
PRINT "10000! を "; B; " で割った商の先頭付近:"
PRINT STR$(Q(QN));
FOR J=QN-1 TO MAX(0, QN-5) STEP-1
PRINT RIGHT$("000000000" & STR$(Q(J)), 9);
NEXT J
PRINT "..."
PRINT "余り:"; rrem
END


Google　AIの方がチャットさんよりもいい感じです。
しかし、今回も有理数モードの方が高速

!' ===========================================
!' 有理数モード　10000階乗/123456（十進BASIC）
!' ===========================================
OPTION ARITHMETIC RATIONAL
LET t0=TIME
PRINT FACT(10000)/123456
PRINT TIME-t0;"seconds elapsed."
END

私の場合は多倍長はAIが頼りです。しかし、本当は

!チュドノフスキー(Chudnovsky)級数
OPTION ARITHMETIC DECIMAL_HIGH
! （Google　AI例: 1000桁モード）

LET t0 = TIME

! 定数の設定
LET C = 640320
LET C3_over_24 = (C^3) / 24 ! 漸化式用の定数
LET S = SQR(C) ! √640320

! 初期値 (k=0 の項)
LET k = 0
LET M = 1 ! 階乗部分の積
LET L = 13591409 ! 分子の 545140134*k + 13591409
LET X = C * S ! 分母の C^(3k+1.5)
LET Z = L / X ! 初項を合計に代入

! ループ計算 (1000桁精度)
FOR k = 1 TO 70 !1回のループで何桁? 精度向上
! 階乗部分の漸化式更新 (factor(6k) / (factor(3k)*factor(k)^3) を更新)
! 前のMに対し、(6k-5)*(6k-4)*(6k-3)*(6k-2)*(6k-1)*(6k) / (k^3 * (3k-2)*(3k-1)*(3k)) を掛けるのに相当
! これを整理すると以下の係数になる
LET M = M * (-(6*k-5)*(2*k-1)*(6*k-1)) / (k^3 * C3_over_24)

LET L = 545140134 * k + 13591409

! 和の更新
LET Z = Z + (M * L) / (C * S)
NEXT k

! 最終計算: 1/pi = 12 * Z
LET my_pi = 1 / (12 * Z)

PRINT my_pi
PRINT PI-my_pi
PRINT TIME - t0; "seconds elapsed."
END

チュドノフスキー(Chudnovsky)級数の円周率を高速に計算する多倍長プログラムをAIに書いてとお願いしていますが、書いてくれません。しばっちさまに2016年10月円周率のプログラムを書いて頂き、今も時々計算結果出力してます。40万桁で約1時間かかります。

Re: Linux版basicとBASファイルの関連付け - SHIRAISHI Kazuo

2026/02/09 (Mon) 15:34:56

-OR open and run
-NR no run
です。
-NRを付加して起動すると，ファイルを開くだけになります。
-ORを付加して起動すると，ファイルを開いて，実行します。

Re: Linux版basicとBASファイルの関連付け - 島村1243

2026/02/10 (Tue) 09:41:51

白石先生、早速にご教示賜り有難うございます。
ご教示頂いたオプション「-NR」を使用したら、目的を達成出来ました！！

ネット検索しても「-OR オプションを使用する。」とのAI回答表示しかなく、途方にくれていました。