2020/03/18

放物螺旋1(Parabolic Spiral)の書き方

放物螺旋とは


 放物螺旋(ほうぶつらせん、Parabolic Spiral)は極座標の方程式 (1) によって表される曲線である。渦は外側にいくほど(t が大きくなるほど)間隔が狭くなっていく。

Wikipediaより引用


数式表現


(1)r(t)=at

書き方


 R言語での書き方は以下のとおり。

 条件:a=1 のとき

1
2
3
4
5
6
7
8
9
10
t=seq(0, 12*pi, length=300)
a=1
r=a*sqrt(t)
x=r*cos(t)
y=r*sin(t)
svg("parabolic1.svg", width=5, height=5)
plot(x, y, type="l", asp=1, lwd=2, col=2)
abline(h=0)
abline(v=0)
dev.off()

0 件のコメント:

コメントを投稿

マンデルブロ集合の彩色方法(5)

06.knit 1 発散判定式を変更する mandelbrot() 内の発散判定式 |zn|>2 を変更する...