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
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
|
from manimlib.imports import*
class cuber(ThreeDScene):
def construct(self):
axes=ThreeDAxes()
cube=Cube()
# cube.scale(1)
cube.shift(RIGHT+DOWN+OUT)
sq3=Square(color=RED, fill_opacity=0.85)
sq3.rotate(PI/2, axis=UP)
sq3.shift(DOWN+OUT+2*RIGHT)
x=TextMobject("x")
y=TextMobject("y")
z=TextMobject("z")
x.rotate(PI/2, axis=RIGHT)
x.rotate(PI/4,axis=OUT)
x.shift(5.8*DOWN)
y.rotate(PI/2, axis=RIGHT)
y.rotate(PI/8,axis=OUT)
y.shift(5.8*RIGHT)
z.rotate(PI/2, axis=RIGHT)
z.rotate(PI/5,axis=OUT)
z.shift(3.2*OUT+0.4*LEFT)
axis_label=VGroup(x,y,z)
v1=Vector(color=YELLOW,buff=15)
v1.rotate(PI/4,axis=RIGHT)
v1.shift(2*RIGHT+1*DOWN+1*OUT)
n1=TextMobject(r"$\vec{n}$",color=YELLOW)
n1.scale(0.8)
n1.rotate(PI/2,axis=RIGHT)
n1.rotate(PI,axis=OUT)
n1.shift(3*RIGHT+1.3*OUT+DOWN)
spaceloc = [[0,0,2],[1,0,2],[-1,0,2],[2,0,2],[-2,0,2],[3,0,2],[-3,0,2],
[0,1,2],[1,1,2],[-1,1,2],[2,1,2],[-2,1,2],[3,1,2],[-3,1,2],
[0,-1,2],[1,-1,2],[-1,-1,2],[2,-1,2],[-2,-1,2],[3,-1,2],[-3,-1,2],
[0,2,2],[1,2,2],[-1,2,2],[2,2,2],[-2,2,2],[3,2,2],[-3,2,2],
[0,-2,2],[1,-2,2],[-1,-2,2],[2,-2,2],[-2,-2,2],[3,-2,2],[-3,-2,2],
[0,3,2],[1,3,2],[-1,3,2],[2,3,2],[-2,3,2],[3,3,2],[-3,3,2],
[0,3,2],[1,3,2],[-1,3,2],[2,3,2],[-2,3,2],[3,3,2],[-3,3,2],
[0,4,2],[1,4,2],[-1,4,2],[2,4,2],[-2,4,2],[3,4,2],[-3,4,2],
[0,4,2],[1,4,2],[-1,4,2],[2,4,2],[-2,4,2],[3,4,2],[-3,4,2],
[0,5,2],[1,5,2],[-1,5,2],[2,5,2],[-2,5,2],[3,5,2],[-3,5,2],
[0,5,2],[1,5,2],[-1,5,2],[2,5,2],[-2,5,2],[3,5,2],[-3,5,2],
[0,6,2],[1,6,2],[-1,6,2],[2,6,2],[-2,6,2],[3,6,2],[-3,6,2],
[0,1.5,2],[1,1.5,2],[-1,1.5,2],[2,1.5,2],[-2,1.5,2],[3,1.5,2],[-3,1.5,2],
[0,3,2],[1,3,2],[-1,3,2],[2,3,2],[-2,3,2],[3,3,2],[-3,3,2]]
veclist1=[Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY)]
[veclist1[i].rotate(PI/4,axis=RIGHT) for i in range(42)]
[veclist1[i].rotate(PI/6,axis=OUT) for i in range(42)]
[veclist1[i].rotate(PI/8,axis=DOWN) for i in range(42)]
vectorfield1=VGroup(*veclist1)
[veclist1[i].shift(spaceloc[i]) for i in range(42)]
veclist2=[Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY)]
[veclist2[i].rotate(PI/4,axis=RIGHT) for i in range(42)]
[veclist2[i].rotate(PI/4,axis=RIGHT) for i in range(42)]
[veclist2[i].rotate(PI/6,axis=OUT) for i in range(42)]
[veclist2[i].rotate(PI/8,axis=DOWN) for i in range(42)]
vectorfield2=VGroup(*veclist2)
[veclist2[i].shift(spaceloc[i]) for i in range(42)]
veclist3=[Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY)]
[veclist3[i].rotate(PI/4,axis=RIGHT) for i in range(42)]
[veclist3[i].rotate(PI/4,axis=RIGHT) for i in range(42)]
[veclist3[i].rotate(PI/6,axis=OUT) for i in range(42)]
[veclist3[i].rotate(PI/8,axis=DOWN) for i in range(42)]
vectorfield3=VGroup(*veclist3)
[veclist3[i].shift(spaceloc[i]) for i in range(42)]
veclist4=[Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY)]
[veclist4[i].rotate(PI/4,axis=RIGHT) for i in range(42)]
[veclist4[i].rotate(PI/4,axis=RIGHT) for i in range(42)]
[veclist4[i].rotate(PI/6,axis=OUT) for i in range(42)]
[veclist4[i].rotate(PI/8,axis=DOWN) for i in range(42)]
vectorfield4=VGroup(*veclist4)
[veclist4[i].shift(spaceloc[i]) for i in range(42)]
vectorfield1.shift(1.5*DOWN)
vectorfield2.shift(IN+1.5*DOWN)
vectorfield3.shift(2*IN+1.5*DOWN)
vectorfield4.shift(3*IN+1.5*DOWN)
vectors=[vectorfield1,vectorfield2,vectorfield3,vectorfield4]
vectorfield=VGroup(*vectors)
vectorfield.scale(1.3)
fv=[Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),Vector(color=GREY),
Vector(color=GREY)
]
spaceloc2 = [[1.5,0.5,0.5],[1.5,1,0.5],[1.5,1.5,0.5],
[1.5,0.5,1],[1.5,1,1],[1.5,1.5,1],
[1.5,0.5,1.5],[1.5,1,1.5],[1.5,1.5,1.5]
]
[fv[i].rotate(PI/4,axis=RIGHT) for i in range(9)]
[fv[i].rotate(PI/6,axis=OUT) for i in range(9)]
[fv[i].rotate(PI/8,axis=DOWN) for i in range(9)]
[fv[i].shift(spaceloc2[i]) for i in range(9)]
fvfield=VGroup(*fv)
fvfield.shift(2*DOWN)
fvfield.scale(1.3)
flux=TextMobject("Flux through one side of the cube").set_color(GOLD)
flux.shift(3.5*UP+0.5*LEFT)
alll=[vectorfield1,vectorfield2,vectorfield3,vectorfield4,fvfield]
alllvectors=VGroup(*alll)
self.set_camera_orientation(phi=70 * DEGREES,theta=-75*DEGREES)
self.play(ShowCreation(axes),ShowCreation(axis_label))
self.begin_ambient_camera_rotation(rate=0.2)
self.play(ShowCreation(alllvectors))
# self.add(fvfield)
self.play(ShowCreation(cube, run_time=1))
self.wait(1)
self.play(ShowCreation(sq3))
self.wait(1)
self.play(FadeOut(cube))
self.play(FadeOut(vectorfield))
self.add_fixed_in_frame_mobjects(flux)
# self.play(ShowCreation(flux))
self.wait(1)
self.play(ShowCreation(v1),ShowCreation(n1))
self.wait(5)
self.stop_ambient_camera_rotation()
self.wait(1)
|
