どのように暗黙的に定義された関数のグラフを、Pythonでプロットすることができますか？

Question

私の関数fは（1、x、y）という形式のベクトルをとり、出力は実数です。

x-y平面でf（1、x、y）== 0となるような線を描きたいと思います。

私の試み：= Y1（...）は、fここ

delta = 0.025 
xrange = arange(0, 12, delta) 
yrange = arange(0, 12, delta) 
p, q = meshgrid(xrange, yrange) 

mesh_point = [] 
for i in range(len(p[0])): 
    for j in range (len(q[0])): 
     mesh_point.append([p[0][i], q[0][j]]) 

for i in range (0, len(p[0])): 
plt.contour(p, q, (y1([1, p[0][i], q[0][i]])-y2([1, p[0][i], q[0][i]])) , [0]) 
plt.show() 

- Y2（...）

誰でも助けることができますか？

Answer 1

f(1,x,y) = 0には、contourに左側を入力し、1つのレベル（0）を選択して等高線を描くことができます。

import matplotlib.pyplot as plt 
import numpy as np 

delta = 0.025 
x = np.arange(0, 12, delta) 
y = np.arange(0, 12, delta) 
p, q = np.meshgrid(x, y) 
# define some function f(n,x,y) 
f = lambda n, x, y: np.sin(np.sqrt(np.abs(x-3-0.5*y)**1.5))+np.exp(-(y-6)**2)-1 
z=f(1, p,q) 

# plot contour line of f(1,x,y)==0 
plt.contour(p, q, z , [0], colors=["k"]) 

#make legend 
proxy, = plt.plot([], color="k") 
plt.legend(handles=[proxy], labels=["f(1,x,y) = 0"]) 
plt.show()