#==========
Julia の修行をするときに,いろいろなプログラムを書き換えるのは有効な方法だ。
以下のプログラムを Julia に翻訳してみる。
MATLAB や Python(matplotlib) の例はよくあるが。
ファイル名: torus.jl 関数名: torus
翻訳するときに書いたメモ
meshgrid って,結局は単純な物だ
==========#
using Plots
function torus(n=100, a=2, b=9)
pyplot()
u = reshape(repeat(range(0, 2π, length=n), n), n, n);
v = u';
x = (b .+ a .* cos.(u)) .* cos.(v);
y = (b .+ a .* cos.(u)) .* sin.(v);
z = a .* sin.(u);
surface(x, y, z, lims=(-10, 10))
end
torus(30)
torus(50, 5, 10)
savefig("fig.png")