#define RND_NBIT 32
#define RAND_MAX 2147483647
#define RND_P 9689
#define RND_Q 4187
void ranset(int raninit, int rseed[]);
void rnd(int r[], int n, int rseed[]);
void drnd(double r[], int n, int rseed[]);