球面調和関数（２）

球面調和関数（１）からの続き

ぜんぜん，球面調和関数にはなっていないけれど。JuliaのPlotsの散布図で色をつける方法がなんとなくわかった。各点の色に対応する配列 c を用意して，カラースキームを指定して，おまじないの marker_z = c，をオプションにすればいいようだ。

とりあえず，マーカーをcricleにして半径を関数値にあわせて変化させる。色についてはこれとは独立な関数を用意すればよい。α値なども同様に変えられるかもしれない。

using Plots
using Random
using ColorSchemes

Random.seed!(0)
N=10000

x, y, z, c, t, s = rand(N).-0.5, rand(N).-0.5, rand(N).-0.5, rand(N), rand(N), rand(N)
r = (x.*x .+ y.*y .+ z.*z).^0.5
t = acos.(z./r)
s = atan.(y./x)
c = 100 .* z.^2 .* exp.((-5).*r) .+ 0.5

#println(x[10],y[10],z[10],c[10])

scatter(x, y, z,
marker_z = c,
markershape = :circle,
markersize = c,
markeralpha = 0.5,
markercolor = :jet,
markerstrokewidth = 0.2,
markerstrokealpha = 0.3,
markestrokecolor = :white,markerstrokestyle = :dot,
legend = :right,
camera = (60,60)
)

図 Julia の Plots の散布図の例