Find the gaussian value of a position

gauss.pde

// loc = distance from center of distribution
float gaussian( float loc, float variance )
{ 
  variance *= variance;
  float a = loc * loc / ( 2.0 * variance );
  float b = exp( -a );
  float c = sqrt( PI * 2.0 * variance );
  return b / c;
}