【曲線-3】フェルマー螺旋を描いてみよう

フェルマー螺旋の数式は下記ですが、
x=a*sqr(th)*cos(th)
y=a*sqr(th)*sin(th)
パラメータがあまりないので、いじりようがないかもしれません。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
%!PS-Adobe-3.0 EPSF-3.0

%%BoundingBox: 0 0 240 240

%%Title:Felmat spiral(フェルマー螺旋)を描く

/a 1 def % 大きさ

/sw 0 def % 始点を置くためのスイッチ

0 0 240 240 rectfill % 黒塗りつぶし背景

120 120 translate % 座標の原点を中央に移動

newpath % パスの初期化

0 1 12000 { % 0から始めて１づつ増分し12000まで{ }内を繰り返す

    /th exch def % 制御変数をth(角度)に入れる

    % フェルマー螺旋の数式

    /x a th sqrt mul th cos mul def 
    /y a th sqrt mul th sin mul def
    % swが０なら始点を置き、さもなくば線を引く

    sw 0 eq { x y moveto /sw 1 def }{ x y lineto } ifelse

} for

.5 .7 .8 sethsbcolor % 色設定

stroke % 線を描画する