El Partenón. Atenas.

Ictinio y Calícrates (arquitectos) y Fidias (escultor). El Partenón. Del 447 al 432 a.c. Atenas.
El Partenón es uno de los templos de la Acrópolis dedicado a la diosa Atenea.

Pulsa el botón Play para iniciar una secuencia de animación que muestra elementos geométricos de la composición de la obra:

(1) En primer lugar recomponemos la cubierta del frontón a partir de los restos para tener una idea de la fachada completa del templo.
(2) Dibujamos una cuadrado que utilizaremos para construir un rectángulo áureo que enmarque la fachada del Partenón. Dentro de ese rectángulo vamos a construir la espiral. Veremos cómo los sucesivos pasos nos llevan a distintos elementos arquitectónicos prefijados por la construcción.
(3) Para cada rectángulo eliminaremos un cuadrado a la vez que construimos el arco correspondiente a ese cuadrado y su centro. En la primera división encontramos que, de las 8 columnas de la fachada, separamos 5 a la izquierda y 3 a la derecha, los tres números: 8, 5 y 3 son elementos de la sucesión de Fibonacci que está íntimamente ligada a la proporción áurea (Φ es el límite de los cocientes entre los términos consecutivos de la sucesión).
(4) Al eliminar el segundo cuadrado, llegamos a una línea horizontal que indica el comienzo del arquitrabe.
(5) El tercero nos aisla una columna -la tercera por la izquierda- y el cuarto es la línea de la cornisa. El nuevo rectángulo áureo abarca ahora la franja que contiene el entablamento (franja situada entre las columnas y el frontón).
(6) Continuamos el proceso de construcción de la espiral. 
(7) Obtenemos el punto de corte de las diagonales (en verde) de los dos primeros rectángulos es el punto de convergencia de la espiral.

En esta dirección http://jmora7.com/Arte/arte.htm se encuentra una colección de cuadros analizados siguiendo esta forma de trabajar.

A = (0.030909090909087417, 2.28)
l = 14 (slider from 0 to 15)
A1: A + (l, 0)
a = 20 (slider from 0 to 20, step 0.001)
a2 = Max(0, Min[1, a - 2])
a1 = Max(0, Min[1, a - 1])
a4 = Max(0, Min[1, a - 5])
a7 = Max(0, Min[1, a - 7])
a8 = Max(0, Min[1, a - 8])
a9 = Max(0, Min[1, a - 9])
a6 = Max(0, Min[1, a - 6])
a11 = Max(0, Min[1, a - 11])
a12 = Max(0, Min[1, a - 12])
a15 = Max(0, Min[1, a - 15])
a16 = Max(0, Min[1, a - 16])
α = 0.4363323129985824 rad (slider from 0 to 6.283185307179586, step 1°)
a18 = Max(0, Min[1, a - 18])
a19 = Max(0, Min[1, a - 19])
a3 = Max(0, Min[1, a - 3])
a5 = Max(0, Min[1.11, a - 5])
a10 = Max(0, Min[1, a - 10])
a13 = Max(0, Min[1, a - 13])
a14 = Max(0, Min[1, a - 14])
b = 11.5 (slider from 0 to 15, step 0.1)
B = (1.1709090909090862, 3.44)
B1: B + (b, 0)
BC = Midpoint(B, B1)
f1 = Line(A, A1)
s1 = OrthogonalLine(BC, f1)
D1 = (-0.14909090909091238, 8.22)
D2 = (11.010909090909074, 10.98)
s2 = Line(D1, D2)
s3 = Mirror(s2, s1)
D3 = (-0.3490909090909121, 8.72)
s4 = Line(D3, s2)
s5 = Mirror(s4, s1)
C1 = Point(s4)
C3 = Intersect(s4, s1)
k1 = Segment(C1, C3)
C2 = Mirror(C1, s1)
D4 = Mirror(C3, s1)
k2 = Segment(C2, D4)
t1 = OrthogonalLine(C3, s1)
t2 = OrthogonalLine(C1, t1)
P1 = Intersect(t1, t2)
P4 = Mirror(P1, s1)
d = Distance(P1, P4)
dd: d / 1.618
P2: P1 - (0, dd)
D5 = Mirror(P2, t2)
t3 = OrthogonalLine(P2, t2)
t4 = OrthogonalLine(P4, t1)
P3 = Intersect(t3, t4)
D6 = Intersect(k1, k2)
h4, i2, i1, i3, i4 = Polygon(P1, P2, P3, P4)
c1 = Circle(P2, P1)
Q2 = Intersect(c1, i1, 1)
r1 = OrthogonalLine(Q2, i1)
c2 = Circle(P3, Q2)
Q3 = Intersect(c2, i3, 1)
r22 = OrthogonalLine(Q3, i3)
c3_1 = Circle(P4, Q3)
Q4 = Intersect(c3_1, i4, 1)
r3 = OrthogonalLine(Q4, i4)
QQ1 = Intersect(r1, i4)
c4 = Circle(Q4, QQ1)
QQ4 = Intersect(c4, r3, 2)
r4 = OrthogonalLine(QQ4, r3)
Q5 = Intersect(r1, r4)
QQ2 = Intersect(r1, r22)
c5 = Circle(Q5, QQ2)
QQ5 = Intersect(c5, r4, 2)
r5 = OrthogonalLine(QQ5, r4)
Q6 = Intersect(r22, r5)
QQ3 = Intersect(r3, r22)
d1 = CircleArc(QQ1, P1, Q2)
d2 = CircleArc(QQ2, Q2, Q3)
d3 = CircleArc(QQ3, Q3, Q4)
d4 = CircleArc(QQ4, Q4, Q5)
d5 = CircleArc(QQ5, Q5, Q6)
c6 = Circle(QQ3, Q6)
Q7 = Intersect(c6, r3, 1)
r2 = OrthogonalLine(QQ3, r3)
r6 = OrthogonalLine(Q7, r3)
N = Intersect(r4, r6)
QQ6 = Intersect(r6, r5)
d6 = CircleArc(QQ6, Q6, Q7)
c7 = Circle(QQ4, Q7)
QQ91 = Intersect(c7, r4, 1)
r7 = OrthogonalLine(QQ91, r4)
QQ7 = Intersect(r6, r7)
polígono1, j1, j2, j3, j4, P6, P8 = Polygon(P1, P2, 4)
C4: C1 + (a1 * (C3 - C1))
g1 = Segment(C1, C4)
C5: C3 + (a2 * (C2 - C3))
g2 = Segment(C3, C5)
P91 = Midpoint(P6, P2)
P92 = Midpoint(j4)
h1 = Segment(P91, P92)
P91' = Rotate(P91, a5, P91)
P92' = Rotate(P92, a5, P91)
a55: (a5 * π / 2)
h2 = Segment(P91, P8)
h3 = Segment(Rotate[P91, (-a5), P91], Rotate[P8, (-a5), P91])
P81 = Rotate(P8, (-a5), P91)
R1 = Rotate(P1, (a7 * π / 2), P8)
e1 = CircleArc(P8, P1, R1)
m1 = Segment(P8, P6)
R = Rotate(P1, (a7 * π / 2), P8)
R2 = Rotate(P6, (a8 * π / 2), QQ2)
m2 = Segment(QQ2, Q3)
e2 = CircleArc(QQ2, P6, R2)
Q36 = NaN
R3 = Rotate(Q3, (a9 * π / 2), QQ3)
e3 = CircleArc(QQ3, Q3, R3)
R4 = Rotate(Q4, (a10 * π / 2), QQ4)
e4 = CircleArc(QQ4, Q4, R4)
R5 = Rotate(Q5, (a11 * π / 2), QQ5)
R6 = Rotate(Q6, (a12 * π / 2), QQ6)
e5 = CircleArc(QQ5, Q5, R5)
e6 = CircleArc(QQ6, Q6, R6)
R7 = Rotate(Q7, (a13 * π / 2), QQ7)
m3 = Segment(Q4, QQ3)
m4 = Segment(Q5, QQ4)
m5 = Segment(QQ5, Q6)
m6 = Segment(QQ6, Q7)
m7 = Segment(QQ7, QQ91)
c = Circle(QQ91, QQ5)
QQ8 = Intersect(c, m7, 1)
R8 = Rotate(QQ91, (a14 * π / 2), QQ8)
S1: P2 + (a15 * (P4 - P2))
l1 = Segment(P2, S1)
S2: P8 + (a16 * (P3 - P8))
l2 = Segment(P8, S2)
e7 = CircleArc(QQ7, Q7, R7)
f = OrthogonalLine(QQ8, i3)
QQ92 = Intersect(f, m5)
e8 = CircleArc(QQ8, QQ91, R8)
e = CircleArc(P91, P81, P8)