帮我找一段绘制Julia分形集的python代码

以下是一个绘制Julia分形集的Python代码示例：

import numpy as np
import matplotlib.pyplot as plt

def julia_set(c, n=100, m=1000, xmin=-1.5, xmax=1.5, ymin=-1.5, ymax=1.5):
    x, y = np.meshgrid(np.linspace(xmin, xmax, m), np.linspace(ymin, ymax, n))
    z = x + y * 1j
    for i in range(100):
        z = z**2 + c
    return z

c = -0.8 + 0.156j
julia = julia_set(c)

plt.imshow(np.abs(julia) < 2, cmap='gray')
plt.xticks([])
plt.yticks([])
plt.show()

在这个例子中，我们使用了numpy和matplotlib库来绘制Julia分形集。函数julia_set()使用了复数运算，循环100次来计算Julia分形集的像素值。最后，我们使用imshow()函数将像素值转换为黑白图像并显示出来