FUNCTION RNDN(q,p,a,b)     ! a: 優秀度（１〜１０）  b: 難易度（１〜１０） →  正解１  負正解０
   LET X=RND
   IF X=<(a-1)/10+0.05+(1.05-a/10)*((5-b)/5) THEN
      LET RNDN=1
   ELSE
      LET RNDN=0
   END IF
END FUNCTION

RANDOMIZE
CALL SetColor(0,0,0)
PRINT
PRINT "  比較的滑らかに分布した能力の解答者と、"
PRINT "  比較的滑らかに分布した難易度の問題における、"
PRINT "ＳＰ表による滑らか分布度の判定 のシミュレーション"
PRINT
PRINT
PRINT "（１） １回目のシミュレーション"
PRINT
DIM T(10,10),TT(10,10),Q(10),QQ(10),R(10),RR(10),WS(10),W(10),WW(10)
LET L1=0
LET L2=0
FOR XX=1 TO 1000
   FOR i=1 TO 10
      FOR j=1 TO 10
         LET T(11-i,j)=RNDN(i,j,i,j)
      NEXT j
   NEXT i
   MAT TT=TRN(T)
   IF XX=1 THEN
      PRINT "元の行列： 行は事前評価で能力の高い順の解答者、列は事前評価で容易な順の問題。"
      PRINT
   END IF
   FOR i=1 TO 10
      LET KEI1=O
      LET KEI2=O
      FOR j=1 TO 10
         IF XX=1 THEN
            PRINT USING "#####":T(i,j);
         END IF
         LET KEI1=KEI1+T(i,j)
         LET KEI2=KEI2+T(i,j)*10^(10-j)
         IF j=10 THEN
            LET Q(i)=KEI1
            LET R(i)=KEI2
            IF XX=1 THEN
               CALL SetColor(0,0,250)
               PRINT USING "#######":Q(i)
               PRINT
            END IF
         END IF
      NEXT j
   NEXT i
   IF XX=1 THEN
      FOR i=1 TO 10
         LET WS(i)=0
         FOR j=1 TO 10
            LET WS(i)=WS(i)+T(j,i)
         NEXT j
      NEXT i
      CALL SetColor(0,0,250)
      FOR i=1 TO 10
         PRINT USING "#####":WS(i);
      NEXT i
      PRINT
      PRINT
   END IF
   CALL SORT(Q,R,T)
   FOR i=1 TO 9
      FOR j=i+1 TO 10
         IF Q(i)=Q(j) THEN
            IF R(i)<R(j) THEN
               LET z1=R(j)
               LET R(j)=R(i)
               LET R(i)=z1
               FOR k=1 TO 10
                  LET z2=T(j,k)
                  LET T(j,k)=T(i,k)
                  LET T(i,k)=z2
               NEXT k
            END IF
         END IF
      NEXT j
   NEXT i
   IF XX=1 THEN
      CALL SetColor(0,0,0)
      PRINT
      PRINT "ＳＰ分析表： 行は結果的能力順解答者、列は事前評価で容易な順の問題。"
      PRINT
   END IF
   FOR i=1 TO 10
      FOR j=1 TO 10
         IF XX=1 THEN
            PRINT USING "#####":T(i,j);
         END IF
      NEXT j
      IF XX=1 THEN
         CALL SetColor(0,0,255)
         PRINT USING "#######":Q(i)
         PRINT
      END IF
   NEXT i
   FOR i=1 TO 10
      LET W(i)=0
      FOR j=1 TO 10
         LET W(i)=W(i)+T(j,i)
      NEXT j
   NEXT i
   IF XX=1 THEN
      CALL SetColor(0,0,255)
      FOR i=1 TO 10
         PRINT USING "#####":W(i);
      NEXT i
      PRINT
      PRINT
      PRINT
   END IF
   LET G11=0
   LET G12=0
   FOR i=1 TO 10
      LET G11=G11+SQR((Q(i)-(10-i))^2)
      LET G12=G12+SQR((Q(i)-(11-i))^2)
   NEXT i
   IF XX=1 THEN
      CALL SetColor(0,0,0)
      PRINT "      �@  解答者の能力の滑らか分布度： ";
      PRINT USING "###.#":ROUND(100-(G11+G21)/2,2);
      PRINT;" %"
      PRINT
      PRINT
      PRINT "元の行列の転置行列： 行は事前評価で容易な順の問題、行は事前評価で能力の高い順の解答者。"
      PRINT
   END IF
    
   FOR i=1 TO 10
      LET EI1=O
      LET EI2=O
      FOR j=1 TO 10
         IF XX=1 THEN
            PRINT USING "#####":TT(i,j);
         END IF
         LET EI1=EI1+TT(i,j)
         LET EI2=EI2+TT(i,j)*10^(10-j)
         IF j=10 THEN
            LET QQ(i)=EI1
            LET RR(i)=EI2
            IF XX=1 THEN
               CALL SetColor(0,0,250)
               PRINT USING "#######":QQ(i)
               PRINT
            END IF
         END IF
      NEXT j
   NEXT i
   FOR i=1 TO 10
      LET WW(i)=0
      FOR j=1 TO 10
         LET WW(i)=WW(i)+TT(j,i)
      NEXT j
   NEXT i
   IF XX=1 THEN
      CALL SetColor(0,0,255)
      FOR i=1 TO 10
         PRINT USING "#####":WW(i);
      NEXT i
      PRINT
      PRINT
   END IF
    
   CALL SORT(QQ,RR,TT)
   FOR i=1 TO 9
      FOR j=i+1 TO 10
         IF QQ(i)=QQ(j) THEN
            IF RR(i)<RR(j) THEN
               LET z1=RR(j)
               LET RR(j)=RR(i)
               LET RR(i)=z1
               FOR k=1 TO 10
                  LET z2=TT(j,k)
                  LET TT(j,k)=TT(i,k)
                  LET TT(i,k)=z2
               NEXT k
            END IF
         END IF
      NEXT j
   NEXT i
   IF XX=1 THEN
      CALL SetColor(0,0,0)
      PRINT
      PRINT "Ｓ'Ｐ'分析表： 行は結果的容易さ順問題、列は事前評価で能力の高い順の解答者。"
      PRINT
   END IF
   FOR i=1 TO 10
      FOR j=1 TO 10
         IF XX=1 THEN
            PRINT USING "#####":TT(i,j);
         END IF
      NEXT j
      IF XX=1 THEN
         CALL SetColor(0,0,250)
         PRINT USING "#######":QQ(i)
         PRINT
      END IF
   NEXT i
    
   IF XX=1 THEN
      CALL SetColor(0,0,250)
      FOR i=1 TO 10
         PRINT USING "#####":WW(i);
      NEXT i
      PRINT
      PRINT
      PRINT
   END IF
   LET G21=0
   LET G22=0
   FOR i=1 TO 10
      LET G21=G21+SQR((QQ(i)-(10-i))^2)
      LET G22=G22+SQR((QQ(i)-(11-i))^2)
   NEXT i
   IF XX=1 THEN
      CALL SetColor(0,0,0) 
      PRINT "      �A  問題の難易度の滑らか分布度： ";
      PRINT USING "###.#":ROUND(100-(G21+G22)/2,2);
      PRINT;" %"
      PRINT 
      PRINT "      �B  総合的な滑らか分布度      ： ";
      PRINT USING "###.#":ROUND(100-(G11+G12+G21+G22)/4,2);
      PRINT " %"
      PRINT
      PRINT
   END IF
   LET L1=L1+(100-(G11+G12)/2)
   LET L2=L2+(100-(G21+G22)/2)
NEXT XX
PRINT
PRINT "（２） １０００回のシミュレーションの結果の平均："
PRINT
PRINT "      解答者の能力の滑らか分布度： ";
PRINT USING "###.#":ROUND(L1/1000,2);
PRINT;" %"
PRINT "      問題の難易度の滑らか分布度： ";
PRINT USING "###.#":ROUND(L2/1000,2);
PRINT" %"
PRINT "      総合的な滑らか分布度      ： ";
PRINT USING "###.#":ROUND((L1+L2)/2000,2);
PRINT " %"
PRINT
END

EXTERNAL SUB SORT(A(),B(),C(,))
FOR i=1 TO 9
   FOR j=i+1 TO 10
      IF A(i)<A(j) THEN
         LET s1=A(j)
         LET A(j)=A(i)
         LET A(i)=s1
         LET s2=B(j)
         LET B(j)=B(i)
         LET B(i)=s2
         FOR k=1 TO 10
            LET s3=C(j,k)
            LET C(j,k)=C(i,k)
            LET C(i,k)=s3
         NEXT k
      END IF
   NEXT j
NEXT i
END SUB

EXTERNAL SUB SetColor(R,G,B)
OPTION CHARACTER Byte
SUB SendMessage(hwnd,msg,wparam,lparam$)
   ASSIGN "USER32.DLL","SendMessageA"
END SUB
LET EM_SETCHARFORMAT=BVAL("0444",16)
LET CHARFORMAT$=CHR$(60) & REPEAT$(CHR$(0),59)
LET CHARFORMAT$(8:8)=CHR$(64)
LET CHARFORMAT$(21:24)=CHR$(R) & CHR$(G) & CHR$(B) & CHR$(0)
CALL SendMessage(WinHandle("RICHEDIT"),EM_SETCHARFORMAT,1, CHARFORMAT$)
END SUB