指定された平均値と標準偏差を持つ正規乱数ベクトルを返す

平均値 mu，標準偏差 sigma を与え，サンプルサイズ n の 1　変数データを返す

function gendat(n::Int64=1000; mu::Float64=0.0, sigma::Float64=1.0)
        x = randn(n)
        (x .- mean(x)) ./ std(x) .* sigma .+ mu
end

x = gendat(10);
using Statistics
println("mean = $(mean(x)) sd = $(std(x))")

mean = -1.3322676295501878e-16 sd = 1.0