Início

A = (2, 1.68)
campo, f, g, h, i, j, l, m, n, p, q, r, s = Polygon((-0.8, 1), (-0.8, 100), (-100, 100), (-100, -100), (-0.8, -100), (-0.8, -1), (1, -1), (1, -100), (100, -100), (100, 100), (1, 100), (1, 1))
A': (x(A), (-y(A)))
anima = 81.1 (slider from 0 to 100, step 0.1)
B1 = PointIn(campo)
B2 = PointIn(campo)
B3 = PointIn(campo)
B4 = PointIn(campo)
B5 = PointIn(campo)
B6 = PointIn(campo)
B7 = PointIn(campo)
B8 = PointIn(campo)
B9 = PointIn(campo)
B10 = PointIn(campo)
B11 = PointIn(campo)
B12 = PointIn(campo)
k = 0.1
r1: 1
puntos: CellRange(B1,B12,"point")
puntosC = Join({puntos, {B1}})
b = Spline(puntosC, 3)
B13: B1
vectoresReg = Sequence(Vector[Element[puntosC, k], Element[puntosC, k + 1]] - (r1 * UnitVector[Element[puntosC, k + 1] - Element[puntosC, k]]), k, 1, 12)
A1: 1
A2: 2
A3 = 3
A4 = 4
A5 = 5
A6 = 6
A7 = 7
A8 = 8
A9 = 9
A10 = 10
A11 = 11
A12 = 12
C1 = Element(vectoresReg, A1)
C2 = Element(vectoresReg, A2)
C3 = Element(vectoresReg, A3)
C4 = Element(vectoresReg, A4)
C5 = Element(vectoresReg, A5)
C6 = Element(vectoresReg, A6)
C7 = Element(vectoresReg, A7)
C8 = Element(vectoresReg, A8)
C9 = Element(vectoresReg, A9)
C10 = Element(vectoresReg, A10)
C11 = Element(vectoresReg, A11)
C12 = Element(vectoresReg, A12)
vectoresAA' = Sequence(Vector[(0, -1)] / abs(Element[puntos, k] - A) + Vector[(0, 1)] / abs(Element[puntos, k] - A'), k, 1, 12)
D1 = Element(vectoresAA', A1)
D2 = Element(vectoresAA', A2)
D3 = Element(vectoresAA', A3)
D4 = Element(vectoresAA', A4)
D5 = Element(vectoresAA', A5)
D6 = Element(vectoresAA', A6)
D7 = Element(vectoresAA', A7)
D8 = Element(vectoresAA', A8)
D9 = Element(vectoresAA', A9)
D10 = Element(vectoresAA', A10)
D11 = Element(vectoresAA', A11)
D12 = Element(vectoresAA', A12)
Centro: Sum[puntos] / 12
vectoresFin = Sequence(Vector[Element[puntos, k] - (-5, 0)] / ((abs(Element[puntos, k] - (-5, 0)))^(2) + 1), k, 1, 12)
F1 = Element(vectoresFin, A1)
F2 = Element(vectoresFin, A2)
F3 = Element(vectoresFin, A3)
F4 = Element(vectoresFin, A4)
F5 = Element(vectoresFin, A5)
F6 = Element(vectoresFin, A6)
F7 = Element(vectoresFin, A7)
F8 = Element(vectoresFin, A8)
F9 = Element(vectoresFin, A9)
F10 = Element(vectoresFin, A10)
F11 = Element(vectoresFin, A11)
F12 = Element(vectoresFin, A12)
v: UnitVector[(-5, 0) - Centro] / 6
vectoresReg2 = Sequence((abs(Element[puntos, k]) / 5 * (Vector[Element[puntos, k], Centro] - ((1.93 * r1) * UnitVector[Centro - Element[puntos, k]]))), k, 1, 12)
G1 = Element(vectoresReg2, A1)
G2 = Element(vectoresReg2, A2)
G3 = Element(vectoresReg2, A3)
G4 = Element(vectoresReg2, A4)
G5 = Element(vectoresReg2, A5)
G6 = Element(vectoresReg2, A6)
G7 = Element(vectoresReg2, A7)
G8 = Element(vectoresReg2, A8)
G9 = Element(vectoresReg2, A9)
G10 = Element(vectoresReg2, A10)
G11 = Element(vectoresReg2, A11)
G12 = Element(vectoresReg2, A12)
B = (1.46, 3.98)
C = (1.46, 0.98)
u = Vector(B, C)
C' = Mirror(C, xAxis)
B' = Mirror(B, xAxis)
w = Vector(B', C')
D = (7, 3)
E = (3, 3)
c = Vector(D, E)
E1: (k * (G1 + D1 + C1 + F1)) + v
E2: (k * (G2 + D2 + C2 + F2)) + v
E3: (k * (G3 + D3 + C3 + F3)) + v
E4: (k * (G4 + D4 + C4 + F4)) + v
E5: (k * (G5 + D5 + C5 + F5)) + v
E6: (k * (G6 + D6 + C6 + F6)) + v
E7: (k * (G7 + D7 + C7 + F7)) + v
E8: (k * (G8 + D8 + C8 + F8)) + v
E9: (k * (G9 + D9 + C9 + F9)) + v
E10: (k * (G10 + D10 + C10 + F10)) + v
E11: (k * (G11 + D11 + C11 + F11)) + v
E12: (k * (G12 + D12 + C12 + F12)) + v