A method for generating and managing random numbers, comprising steps of: an input unit setting an N bits binary series R, where N is an integer and N>2, and multi-dimensional coordinate information i for K dimensions (i1,..., ik,..., iK) where ik is an integer and ik ≧0, k:1, 2,..., K, wherein R is a decimal portion of an initial value x 0 for a nonlinear function, x t +i = 4 x t (1 - x t) (hereinafter referred to as LMAP, 0 < x t <1, t=0,1,2,...), and said coordinate information i is then transformed into coordinate value of each dimension (i1,..., ik,..., iK) by a transforming circuit; and x 0, i1,..., ik,..., iK are stored in a register; 1) according to the initial value x 0 and ik, an LMAP circuit implementing LMAP calculation having the accuracy of N bits through a fixed decimal point calculation, which is repeatedly done by chaos computing to generate N bits binary series B ik, and bits of B ik are individually constituted by b k,0, b k, 1,... b k,N -1, wherein b k,0 = [2 × x N × ik], b k,1 = [2 × x N × ik +1],...and b k,N -1 = [2 × x N × ik + N -1], and [] represents a calculation of removing decimal portion; 2) a chaos computing unit implementing the chaos computing on B ik as a decimal portion of the initial value x 0 after LMAP and repeatedly conducting calculation of the LMAP to generate and store an N bits binary series R i 1,... ik, where R i 1,..., ik is R i 1 when k=1 and R i 1,.., ik is R i 1,..., ik when k=K, wherein the bits of R i 1,..., ik are individually constituted by r k,0, r k,i,... and r k,N -1, and r k,0 = [2 × x N], r k,1 =[2 × x N + i],..., r k, N -1 = [2 × x 2 N -1], which are stored in a random number register; 3) using said R i 1,..., ik stored in the random number register as a decimal portion of initial value x 0 of the LMAP and storing it to said register, except when k=K; and a random number generation control unit implementing the above operations 1), 2) and 3) for each dimension k=1, 2,..., K, but when k=K, operation 3) is stopped, and operation 2) is used to generate R i 1,..., ik as random numbers. an input unit setting an N bits binary series R, where N is an integer and N>2, and multi-dimensional coordinate information i for K dimensions (i1,..., ik,..., iK) where ik is an integer and ik ≧0, k:1, 2,..., K, wherein R is a decimal portion of an initial value x 0 for a nonlinear function, x t +i = 4 x t (1 - x t) (hereinafter referred to as LMAP, 0 < x t <1, t=0,1,2,...), and said coordinate information i is then transformed into coordinate value of each dimension (i1,..., ik,..., iK) by a transforming circuit; and x 0, i1,..., ik,..., iK are stored in a register; 1) according to the initial value x 0 and ik, an LMAP circuit implementing LMAP calculation having the accuracy of N bits through a fixed decimal point calculation, which is repeatedly done by chaos computing to generate N bits binary series B ik, and bits of B ik are individually constituted by b k,0, b k, 1,... b k,N -1, wherein b k,0 = [2 × x N × ik], b k,1 = [2 × x N × ik +1],...and b k,N -1 = [2 × x N × ik + N -1], and [] represents a calculation of removing decimal portion; 2) a chaos computing unit implementing the chaos computing on B ik as a decimal portion of the initial value x 0 after LMAP and repeatedly conducting calculation of the LMAP to generate and store an N bits binary series R i 1,... ik, where R i 1,..., ik is R i 1 when k=1 and R i 1,.., ik is R i 1,..., ik when k=K, wherein the bits of R i 1,..., ik are individually constituted by r k,0, r k,i,... and r k,N -1, and r k,0 = [2 × x N], r k,1 =[2 × x N + i],..., r k, N -1 = [2 × x 2 N -1], which are stored in a random number register; 3) using said R i 1,..., ik stored in the random number register as a decimal portion of initial value x 0 of the LMAP and storing it to said register, except when k=K; and a random number generation control unit implementing the above operations 1), 2) and 3) for each dimension k=1, 2,..., K, but when k=K, operation 3) is stopped, and operation 2) is used to generate R i 1,..., ik as random numbers. 1) according to the initial value x 0 and ik, an LMAP circuit implementing LMAP calculation having the accuracy of N bits through a fixed decimal point calculation, which is repeatedly done by chaos computing to generate N bits binary series B ik, and bits of B ik are individually constituted by b k,0, b k, 1,... b k,N -1, wherein b k,0 = [2 × x N × ik], b k,1 = [2 × x N × ik +1],...and b k,N -1 = [2 × x N × ik + N -1], and [] represents a calculation of removing decimal portion; 2) a chaos computing unit implementing the chaos computing on B ik as a decimal portion of the initial value x 0 after LMAP and repeatedly conducting calculation of the LMAP to generate and store an N bits binary series R i 1,... ik, where R i 1,..., ik is R i 1 when k=1 and R i 1,.., ik is R i 1,..., ik when k=K, wherein the bits of R i 1,..., ik are individually constituted by r k,0, r k,i,... and r k,N -1, and r k,0 = [2 × x N], r k,1 =[2 × x N + i],..., r k, N -1 = [2 × x 2 N -1], which are stored in a random number register; 3) using said R i 1,..., ik stored in the random number register as a decimal portion of initial value x 0 of the LMAP and storing it to said register, except when k=K; and a random number generation control unit implementing the above operations 1), 2) and 3) for each dimension k=1, 2,..., K, but when k=K, operation 3) is stopped, and operation 2) is used to generate R i 1,..., ik as random numbers.