平面代数曲線による一筆書き

以下のネタは https://mathtod.online/@satie/263761 より

In [1]:
# 以下は必ず実行しておく. PyPlot も宣言しておかないと、plot_implict で kernel が死ぬ.

using SymPy
using PyPlot
In [2]:
# 数式処理用の変数を宣言

@vars x y
Out[2]:
(x,y)
In [3]:
f=x*(x^2-1)*(x^2-4)*(x^2-9)*(x^2-16)-(1+0/10)*y*(y^2-1)*(y^2-4)*(y^2-9)*(y^2-16)
Out[3]:
$$x \left(x^{2} - 16\right) \left(x^{2} - 9\right) \left(x^{2} - 4\right) \left(x^{2} - 1\right) - 1.0 y \left(y^{2} - 16\right) \left(y^{2} - 9\right) \left(y^{2} - 4\right) \left(y^{2} - 1\right)$$
In [4]:
plot_implicit(f,(x,-10,10),(y,-10,10))
Out[4]:
PyObject <sympy.plotting.plot.Plot object at 0x000000001F630DD8>
In [5]:
g=x*(x^2-1)*(x^2-4)*(x^2-9)*(x^2-16)-(1+5/10)*y*(y^2-1)*(y^2-4)*(y^2-9)*(y^2-16)
Out[5]:
$$x \left(x^{2} - 16\right) \left(x^{2} - 9\right) \left(x^{2} - 4\right) \left(x^{2} - 1\right) - 1.5 y \left(y^{2} - 16\right) \left(y^{2} - 9\right) \left(y^{2} - 4\right) \left(y^{2} - 1\right)$$
In [6]:
plot_implicit(g,(x,-10,10),(y,-10,10))
Out[6]:
PyObject <sympy.plotting.plot.Plot object at 0x000000001F656208>