Thanks to visit codestin.com
Credit goes to github.com

Skip to content

Commit df9eda9

Browse files
Revert "use np to solve equation"
This reverts commit 115772a
1 parent 115772a commit df9eda9

File tree

1 file changed

+2
-43
lines changed

1 file changed

+2
-43
lines changed

31.py

+2-43
Original file line numberDiff line numberDiff line change
@@ -1,16 +1,10 @@
11
from PIL import Image
22
from PIL import ImageSequence
3-
import numpy as np
4-
import matplotlib.pyplot as plt
53

64

75
def main():
86
pic_path = "ufo/mandelbrot.gif"
9-
# first_step(pic_path)
10-
# second_step()
11-
# plt.imshow(mandelbrot(1000, 1000))
12-
# plt.show()
13-
numPyVersion()
7+
first_step(pic_path)
148

159

1610
def first_step(input):
@@ -24,44 +18,9 @@ def first_step(input):
2418
height = 480
2519
for x in range(width):
2620
for y in range(height):
27-
print(im.getpixel((x, y)), end=" ")
21+
print(im.getpixel((x, y)), end = " ")
2822
print()
2923

3024

31-
def second_step():
32-
a = np.arange(15).reshape(3, 5)
33-
print(a)
34-
print(a.shape)
35-
print(a.ndim)
36-
print(a.dtype.name)
37-
print(a.dtype)
38-
print(a.itemsize)
39-
print(type(a))
40-
41-
b = np.array([6, 7, 8])
42-
print(b)
43-
print(b.shape)
44-
pass
45-
46-
47-
def mandelbrot(h, w, maxit=20):
48-
y, x = np.ogrid[-1.4:1.4:h * 1j, -2:0.8:w * 1j]
49-
c = x + y * 1j
50-
z = c
51-
divtime = maxit + np.zeros(z.shape, dtype=int)
52-
53-
for i in range(maxit):
54-
z = z ** 2 + c
55-
diverge = z * np.conj(z) > 2 ** 2
56-
div_now = diverge & (divtime == maxit)
57-
divtime[div_now] = i
58-
z[diverge] = 2
59-
return divtime
60-
61-
62-
def numPyVersion():
63-
print(np.version.version)
64-
65-
6625
if __name__ == "__main__":
6726
main()

0 commit comments

Comments
 (0)