$B!|!|(B $BM=A*LdBj$NFbMF!J5vMF;D:9(B : 1e-11$B!K(B $B!|!|(B
$BM=A*$N9TNs$O!"0J2<$N<0$G@8@.$9$k5?;wMp?t$H$9$k!#(B

int ir;

double random(){
    ir = ( 1229 * ir + 351750 ) % 1664501 ;
    return (double)ir / 1664501.0 ;
}

$BM=A*LdBj$N9TNs$NBg$-$5(B:n $B5Z$S!"Mp?t$N<o(B:ir $B$O(B
key1   n = 1023 , ir = 2001
key2   n = 1024 , ir =   95
key3   n = 1025 , ir = 1000

$B5?;wMp?t$r0J2<$N=gHV$G!"9TNs$*$h$S%Y%/%H%k$K3JG<$9$k!#(B

                                       $B!J9TNs#A$r@8@.$7=*$($F$+$i!K(B
    $B9TNs#A(B[][] $B$KMp?t$r3JG<$9$k=gHV(B     $B%Y%/%H%k#B(B[]$B$K3JG<$9$k=gHV(B

                $BNs"*(B
       0   1   2   3   4   5   6 ...              1
      10  11  12  13  14  15  16 ...              2
$B9T(B    20  21  22  23  24  25  26 ...              3
$B"-(B              $B!'(B                               $B!'(B          	
                $B!'(B                               $B!'(B          	
     110 111 112 113 114 115 116 ...              10
     120 121 122 123 124 125 126 ...              11
                $B!'(B                                $B!'(B          	



$B!|!|(B $BK\A*LdBj$NFbMF!J5vMF;D:9(B : 1e-10$B!K(B $B!|!|(B
$BK\A*$N%G!<%?$O(B key1 key3 $B$OM=A*$HF1$8Mp?t$r;H$$!"(Bkey2 $B$O!"@QJ,JQ49$N5UJQ(B
$B49$rN%;62=$7$?$b$N$H$9$k(B

key1   n = 1000 , ir = 1012
key3   n = 2001 , ir = 1013

key2$B$O!"@QJ,JQ49(B $B"i(BK(x,y)f(y)dy = g(x) $B$N5UJQ49$rN%;62=$7$F5a$a$k9TNs$G(B
$B$9!#K\%3%s%F%9%H$G$O!"(BK(x,y)$B$H$7$F!"(BLorenz$B7?$NJ,2rG=4X?t$N(Bkernel 

	                a
	K(x,y) = -----------------
	         ( x - y )^2 + a^2

$B$^$?!"1&JU$H$7$F!"(B

	           sin 10x
	g(x) = -----------------
	       (x - 0.5)^2 + c^2

$B$r9M$($^$7$?!#=>$C$F!"9TNs(B A[i][j] $B5Z$S!"%Y%/%H%k(B b[] $B$O(B

	           1            1
	A[i][j] = --- * -------------------
	           n    (x[i]-y[i])^2 + a^2

	          sin(10 * x[i])
	b[i] = ---------------------
	       (x[i] - 0.5)^2 + c^2

$B$H$J$j$^$9!#$?$@$7!"9TNs%5%$%:$r(B n $B$KBP$7$F(B

	        j                (1 - 2a)
	y[j] = ---- , x[i] = a + -------- i
	        n                   n

	a = 0.1 , b = 0.2
$B$H$7$^$9!#Nc$($P!"(Bn = 10 $B$N;~(B $B$N(B A[][]$B$NCM$H(B b[]$B$NCM$r0J2<$K<($7$^$9!#(B

A[0][0]  A[0][1]  A[0][2]  ... b[0]
A[1][0]  A[1][1]  A[1][2]  ... b[1]
A[2][0]  A[2][1]  A[2][2]  ... b[2]
$B$H$7$F!"(B

 500.00 12.20  2.76  1.19  0.66  0.42  0.29  0.21  0.16  0.13 0.36
  8.50 609.76 11.68  2.71  1.17  0.65  0.41  0.29  0.21  0.16 4.55
  2.35  8.82 735.29 11.19  2.65  1.16  0.64  0.41  0.28  0.21 6.98
  1.08  2.40  9.16 862.07 10.73  2.60  1.14  0.64  0.41  0.28 1.29
  0.62  1.10  2.44  9.52 961.54 10.31  2.54  1.12  0.63  0.40 -15.52
  0.40  0.62  1.11  2.49  9.90 1000.00  9.90  2.49  1.11  0.62 -23.97
  0.28  0.40  0.63  1.12  2.54 10.31 961.54  9.52  2.44  1.10 -6.02
  0.21  0.28  0.41  0.64  1.14  2.60 10.73 862.07  9.16  2.40 7.99
  0.16  0.21  0.28  0.41  0.64  1.16  2.65 11.19 735.29  8.82 7.88
  0.13  0.16  0.21  0.29  0.41  0.65  1.17  2.71 11.68 609.76 2.50
