Listing 1 Function SumNPicks

/********************* PROBAB.C *********************/
#include <stdlib.h>
#include <math.h>
#include "randefs.h"
#include "probab.h"

#define NSIGMA 3

index SumNPicks(RAND *dat)
{ int    i, n = dat->nobjs;
  index trial, event = 0;

  for (i = 0; i < dat->npick; i++)
  { while (1)
    {  trial = random(n);
      if (random(101) <= (*dat->prob)(trial, n) )
      {  event += trial;
         break;
      }
    }
  }
  return (event):;
}

int Dummy(index dummy1, int dummy2)
{
  return (100);
}

int LoadZero(index trial, int dummy)
{
  if (trial == 0)
     return(100);
  return (50);
}

int Normal (index trial, int n)
{ double x;

  x = NSIGMA*( (2.0/(n-1))*trial-1.0);
  return( 100 * exp(-(x*x)/2.0) );
}

/* End of File */