3

Facebookでこの画像を見つけて、何か新しいことを学ぶ良い機会になると思ったので、それをプロットしたいと思いました。

愛の公式
(出典:akamaihd.net

コード:

import numpy as np
import scipy.special
import pylab

x = np.linspace( -1.0, 1.0, 200 )
y = np.sqrt(  1-np.square(np.complex64(x) ) ) + np.complex64( scipy.special.cbrt( np.square( x ) ) )
pylab.plot( np.complex64( x ), y )
pylab.show()

出力:

IMG

図のようなプロットが得られません。プロットするために、数式をフォームに並べ替えて、y = f(x)上記のコードに変換しました。

4

1 に答える 1

7
import matplotlib.pyplot as plt
import numpy as np

y, x = np.ogrid[-1:2:100j, -1:1:100j]
plt.contour(x.ravel(), y.ravel(), x**2 + (y-((x**2)**(1.0/3)))**2, [1])
plt.axis('equal')
plt.show()

ここに画像の説明を入力


アスペクト比をいじって、曲線をもう少しハートのようにすることができます。

fig = plt.figure()
ax = fig.add_subplot(1, 1, 1)
y, x = np.ogrid[-1.5:2:1000j, -2:2:1000j]
plt.contour(x.ravel(), y.ravel(), x**2 + (y-((x**2)**(1.0/3)))**2, [1])
ax.set_aspect(0.75) 
plt.show()

ここに画像の説明を入力

于 2013-02-15T19:09:19.373 に答える