uniformball.jl

using Plots,Random, Distributions
function uniformball(n)
# 1-säteisessä pallossa n tasajakautunutta pistettä
# Mullerin menetelmä
# Juhani Kaukoranta 9.1.2023
    r = cbrt.(rand(n)) # tasajakauman kuutiojuuresta säde r
    d = Normal(0,1); # normaalijakama μ=0,σ=1.0
    u = rand(d,n);     # koordinaatit normaalijakautuneita
    v = rand(d,n);
    w = rand(d,n);
    normi = @. sqrt(u^2+ v^2+ w^2); #norm = (u*u + v*v + w*w)**(0.5)
    # lasketaan tasajakautuneet x,y,z-koordinaatit
    x =   r .* u ./ normi ; # n kpl x:n arvoa x[1],x[2],x[3],...
    y =   r .* v ./ normi ; # n kpl y:n arvoa
    z =   r .* w ./ normi ; # n kpl z:n arvoa
    plt3d= Plots.plot(x,y,z,seriestype=:scatter,markersize=1,legend=false)
    display(plt3d) # plottaa jakauman
end