Fast Random Number for Sin^HCos^M

P(x) = sin^H(Pi.y).cos^M(Pi.y).dy
= sqrt^H-1(1-x²).x^M.dx

x-> [0<xini,xfin<1[
y -> [-1<yini:yfin<1[
H=real>=1, M=even integer

DOWNLOAD:
To use Sin_Cos.c you must have SU2.res= table. You can download the executable create_Sin_Cos_res.out to create it.

Sin_Cos.c : random number in C
create_Sin_Cos.out : (Linux, 20 Ko): can choose H, M, xmin, xmax, rate
create_Sin_Cos.exe : (Windows, 20 Ko): can choose H, M, xmin, xmax, rate

PROBABILITIES PERFORMANCES :

Performance: 2.3 times faster than the uniform method, 2 times faster than the metropolis.
Figure: The probability P(x) for M=0, H= from 2 to 20.
Figure: The probability P(x) for M=from 2 to 20, H= from 2 to 20.

• mail to D. Loison

Fast Random Number for SinHCosM

Fast Random Number for Sin^HCos^M