Write out the Gaussian equation.
Implement gaussian(x, mean, stddev) to return the Gaussian probability density at x: 1 / (stddev * sqrt(2π)) * exp(-0.5 * ((x - mean) / stddev)²). Inputs are numbers with stddev > 0; return a floating-point value. Results are graded within a small numerical tolerance.
def gaussian(x, mean, stddev):