廢話不多，需要安裝python，和matplotlib、numpy等庫，安裝不全的，自行安裝一下

在windows下測試通過

# -*- coding: utf-8 -*-

#需要安裝matplotlib、numpy等庫才能執行
import numpy as np
import pylab as pl
import time
from matplotlib import cm

def iter_point2(x,y):
    #初始化公式的引數，需要改的話，可以改引數
    a = 1
    b = 1.75
    iter = 0 #初始化迭代次數
    for i in xrange(1, 10):
        iter = i
        dist = (x*x + y*y)  #計算摸長，看是否超過最大限度
        dist *= dist
        #print "dist:" + str(dist)
        if  dist > 200:  #如果超出了最大長度，就跳出迴圈，返回這個迭代次數
            break
        #臨時儲存一下x和y
        tempx = x
        tempy = y
        #這裡是公式 X = a-b|x| + y ; Y = 0.3x;
        x = a - b*abs(tempx) + tempy
        y = 0.3*tempx
    return iter
    
#繪製圖形時，以cx，xy為中心，距離為d
def draw_lozi(cx, cy, d):
    size = 400
    x0, x1, y0, y1 = cx-d, cx+d, cy-d, cy+d 
    y, x = np.ogrid[y1:y0:size*1j, x0:x1:size*1j] #使用範圍生成陣列，後面用這個進行迭代
    c = x + y*1j
    x.shape = -1 #轉化成線性陣列
    y.shape = -1

    start = time.clock()
    lozi = np.ones(c.shape)
    #遍歷每一個點，計算迭代次數，賦值給陣列lozi
    for j in range(0,size):
        for i in range(0,size):
            lozi[j][i] = iter_point2(x[i],y[j])
            pass

    print "time=",time.clock() - start
    #使用陣列lozi，繪圖， 使用藍色調色盤，顯示到圖上的座標範圍是x0,x1,y0,y1
    pl.imshow(lozi, cmap=cm.Blues_r, extent=[x0,x1,y0,y1])
    #不顯示橫縱座標
    pl.gca().set_axis_off()
    #重新整理畫布
    pl.show()
    
#滑鼠點選觸發執行的函式
def on_press(event):
    global g_size
    print event
    print dir(event)
    newx = event.xdata
    newy = event.ydata
    print newx
    print newy
    
    #不合理的滑鼠點選，直接返回，不繪製
    if newx == None or newy == None  or event.dblclick == True:
        return None
    #不合理的滑鼠點選，直接返回，不繪製
    if event.button == 1:  #button ==1 代表滑鼠左鍵按下， 是放大影象
        g_size /= 2
    elif event.button == 3: #button == 3 代表滑鼠右鍵按下， 是縮小影象
        g_size *= 2
    else:
        return None
    print g_size
    
    draw_lozi(newx,newy,g_size)
fig, ax = pl.subplots(1)

g_size = 4.5

#註冊滑鼠事件
fig.canvas.mpl_connect('button_press_event', on_press)

#初始繪製一個圖
draw_lozi(0,0,g_size)