из pylab импорт * из numpy импорт NaN из matplotlib.colors импорт LogNormс = - .123 , .745res = .001 # разрешение сетки max_it = 1000 # максимальное количество итерацийфигура ( размер_фиг. = ( 10 , 10 ))def jul ( a ): z = a для n в диапазоне ( 1 , max_it + 1 ): z = z ** 2 + s [ 0 ] + 1 j * s [ 1 ] если abs ( z ) > 2 : вернуть n вернуть NaNX = диапазон ( -2 , 2 + разрешение , разрешение ) Y = диапазон ( -1,5 , 1,5 + разрешение , разрешение ) Z = нули ( ( len ( Y ), len ( X ) )) для iy , y в перечислении ( Y ): print ( iy + 1 , "of" , len ( Y )) для ix , x в перечислении ( X ): Z [ - iy - 1 , ix ] = jul ( x + 1 j * y )сохранить ( "douady" , Z ) # сохранить массив в файлimshow ( Z , cmap = plt.cm.viridis , interpolation = 'none' , norm = LogNorm ( ), extends = ( X.min ( ) , X.max ( ), Y.min ( ) , Y.max ())) xlabel ( "Re(z)" ) ylabel ( " Im (z)" ) axis ( ( - 1.5 , 1.5 , -1.5 , 1.5 ) ) savefig ( " douady_rabbit_python.svg " ) show ( )