Y=VL_GAUSSIAN(X) computes the standard (zero mean, unit variance) Gaussian density.

To obtain the Gaussian density of standard deviation S do

  Y = 1/S * VL_GAUSSIAN(X/S).