-
Notifications
You must be signed in to change notification settings - Fork 1
/
fig5.py
91 lines (71 loc) · 2.04 KB
/
fig5.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
import numpy as np
import matplotlib.pyplot as plt
from polygon import polygon
from diskmap import diskmap
plt.figure(figsize=(6.4, 9.6))
r = np.arange(0.1, 1, 0.07)
theta = np.arange(0, 2*np.pi, np.pi/12)
p = polygon([1-0.5j, np.inf, 1j, np.inf, -1, np.inf, 1-1j],
[3/2, -1/4, 3/2, 0, 3/2, -1/4, 1])
m = diskmap(p)
plt.subplot(3,2,1)
m.plot(r, theta, 'b', lw=0.5)
p.plot('r')
plt.axis('equal')
plt.axis('off')
plt.axis([-2.5, 2.5, -3, 2])
############################################################
p = polygon([1-1j, np.inf, 1+1j, np.inf, -0.3+1.3j,
np.inf, -1+0.5j, np.inf, -1-1j],
[1, -1/6, 2, -1/3, 2, -1/4, 2, -1/4, 1])
m = diskmap(p)
plt.subplot(3,2,2)
m.plot(r, theta, 'b', lw=0.5)
p.plot('r')
plt.axis('equal')
plt.axis('off')
plt.axis([-2, 2.3, -1.02, 3.28])
############################################################
p = polygon([-1j, np.inf, 0], [1/2, -1, 3/2])
m = diskmap(p, 1)
plt.subplot(3,2,3)
m.plot(r, theta, 'b', lw=0.5)
p.plot('r')
plt.axis('equal')
plt.axis('off')
plt.axis([-2, 2, -1.02, 2.98])
############################################################
p = polygon([-1j, 1-1j, np.inf, 0],
[1/2, 2, -2, 3/2])
m = diskmap(p, 1)
plt.subplot(3,2,4)
m.plot(r, theta, 'b', lw=0.5)
p.plot('r')
plt.axis('equal')
plt.axis('off')
plt.axis([-2,3,-2,3])
############################################################
p = polygon([-2j, np.inf, 1.5-0.5j, np.inf, 2j,
np.inf, -1+0.5j, np.inf, -3j],
[2, -1/6, 2, -5/6, 2, -1/3, 2, -2/3, 1])
m = diskmap(p)
plt.subplot(3,2,5)
m.plot(r, theta, 'b', lw=0.5)
p.plot('r')
plt.axis('equal')
plt.axis('off')
plt.axis([-3,3,-3,3])
############################################################
p = polygon([-3-1j, -3-3j, np.inf, -3+3j, -3+1j,
-1+1j, -3+1j, -3-1j, -2-1j],
[1/2, 1/2, 0, 1/2, 1/2, 2, 1/2, 1/2, 2])
m = diskmap(p)
plt.subplot(3,2,6)
m.plot(r, theta, 'b', lw=0.5)
p.plot('r')
plt.axis('equal')
plt.axis('off')
plt.axis([-3.02, 3.02, -3.02, 3.02])
plt.tight_layout()
plt.savefig('fig5.eps')
plt.show()