Write a function to simulate a specific probability distribution.
Implement simulate_normal(n, mean, stddev, uniforms) to generate n normally distributed samples using the Box-Muller transform. Use consecutive pairs from uniforms, round every result to 6 decimal places, and return the samples in order. For an odd n, discard the second value from the final generated pair.
Signature: def simulate_normal(n, mean, stddev, uniforms):
Each uniforms value is in (0, 1], and its length is at least 2 * ceil(n / 2). Return a list of n numbers.
def simulate_normal(n, mean, stddev, uniforms):