(1)求一个随机数的平均数:
#include<stdio.h>
#define INITIAL_SEED 17
#define MULTIPLIER 25173
#define INCREMENT 13849
#define MODULUS 65536
#define FLOATING_MODULUS 65536.0
static unsigned seed=INITIAL_SEED;
double probability(void)
{
seed=(MULTIPLIER*seed+INCREMENT)%MODULUS;
return(seed/FLOATING_MODULUS);
}
int main(void)
{
double a=0,b,v;
int i;
for(i=1;i<=100;i++)
{
v=probability();
printf("v=%lf/n",v);
a=a+v;
}
b=a/100;
printf("b=%lf",b);
printf("/n");
}
(2)求1000个随机数的平均数
#include<stdio.h>
#define INITIAL_SEED 17
#define MULTIPLIER 25173
#define INCREMENT 13849
#define MODULUS 65536
#define FLOATING_MODULUS 65536.0
static unsigned seed=INITIAL_SEED;
double probability(void)
{
seed=(MULTIPLIER*seed+INCREMENT)%MODULUS;
return(seed/FLOATING_MODULUS);
}
int main(void)
{
double a=0,b,v;
int i;
for(i=1;i<=1000;i++)
{
v=probability();
printf("v=%lf/n",v);
a=a+v;
}
b=a/1000;
printf("b=%lf",b);
printf("/n");
}