Ex 29 : méthode de Héron

un

t/un

un+1=1/2*(un+t/un)

t/un+1

Valeur approchée de

Commencer par saisir la partie entière de racine(t),
 puis modifier le curseur en saisissant la valeur 
suivante  un+1                                    et ainsi de suite ... 
Le rectangle devient un peu plus carré et on 
obtient une valeur approchée de  racine(t)

un = 3.741 (slider from 0 to 20, step 0.001)
A = Point(xAxis)
B = Point(Circle[A, un])
a = Segment(A, B)
t = 14 (slider from 0 to 25, step 1)
c = Circle(A, t / un)
b = OrthogonalLine(A, xAxis)
D = Intersect(c, b, 2)
d = OrthogonalLine(D, b)
e = OrthogonalLine(B, a)
C = Intersect(d, e)
f = Segment(A, D)
g = Segment(D, C)
h = Segment(B, C)
E = (-2.4179285704915507, -5.27328298293869)
F = Point(Circle[E, (1 / 2 * (un + t / un))])
i = Segment(E, F)
k = Circle(E, t / i)
j = OrthogonalLine(E, i)
H = Intersect(k, j, 2)
l = OrthogonalLine(H, j)
m = OrthogonalLine(F, i)
G = Intersect(l, m)
n = Segment(E, H)
p = Segment(H, G)
q = Segment(F, G)
poly1, a_1, b_1, d_1, c_1 = Polygon(A, B, C, D)
poly2, e_1, f_1, h_1, g_1 = Polygon(E, F, G, H)
texte3 = LaTeX(i, true, true)