仮想計算機構

IT業界と無縁な派遣社員のブログ

【Sage】Trefoil knot の可視化

Parametric surfacewikipedia に載っている CG を Sage で再現した。
参考:https://en.wikipedia.org/wiki/Parametric_surface

環境

OS:Windows 10
Python 3.7.7
Sage version 9.2

プログラム

u,v=var('u,v')
r = 5
parametric_plot3d((r * sin(3 * u) / (2 + cos(v)),
                   r * (sin(u) + 2 * sin(2 * u)) / (2 + cos(v + pi * 2 / 3)),
                   r / 2 * (cos(u) - 2 * cos(2 * u)) * (2 + cos(v)) * (2 + cos(v + pi * 2 / 3)) / 4),
                   (u,-pi,3*pi),
                   (v,-pi,3*pi),
                  plot_points=[100,100],
                  opacity=0.3)

実行結果