mu = 1; sd = 0.5;
    x = randn(1000,1) * sd + mu;