#include <stdlib.h>
#include "distrib.h"

float rand_normal(float mean, float sd)
{
  int i;
  float rn = 0;
  for (i=0; i<12; i++)
    rn +=(float)random()/(float)RAND_MAX;
  return sd*(rn-6.0)+mean;
}

float rand_uniform(float low, float high)
{
  return ((float)random()/(float)RAND_MAX)*(high-low)+low;
}