[B][COLOR=DarkSlateGray]=(ACOS(1-2*H/D)-SIN(2*ACOS(1-2*H/D))/2)*D^2*L/4[/COLOR][/B]
[B][COLOR=DarkSlateGray]=(ACOS(1-2*H/D)-SIN(2*ACOS(1-2*H/D))/2)/PI()[/COLOR][/B]
[B][COLOR=DarkSlateGray]Function vol(r#)
Application.Volatile [COLOR=DarkOrange]'Facultatif[/COLOR]
Dim u0#, u1#, u2#, y0#, y1#, y2#
r = r * 3.14159265358979
u1 = 1.6: y1 = r - u1 + Sin(2 * u1) / 2
u2 = 1.5: y2 = r - u2 + Sin(2 * u2) / 2
Do Until Abs(y2) < 0.0000000000001 Or y2 = y1
u0 = u1: u1 = u2: y0 = y1: y1 = y2
u2 = u0 - y0 * (u1 - u0) / (y1 - y0): y2 = r - u2 + Sin(2 * u2) / 2
Loop
vol = (1 - Cos(u2)) / 2
End Function[/COLOR][/B]
[B][COLOR=DarkSlateGray]=vol(p)*D[/COLOR][/B]
[COLOR=DarkSlateGray][B]Function boule(r#)
Application.Volatile [COLOR=DarkOrange]'Facultatif[/COLOR]
Dim u0#, u1#, u2#, y0#, y1#, y2#
u1 = 0.6: y1 = r - u1 ^ 2 * (3 - 2 * u1)
u2 = 0.5: y2 = r - u2 ^ 2 * (3 - 2 * u2)
Do Until Abs(y2) < 0.0000000000001 Or y2 = y1
u0 = u1: u1 = u2: y0 = y1: y1 = y2
u2 = u0 - y0 * (u1 - u0) / (y1 - y0): y2 = r - u2 ^ 2 * (3 - 2 * u2)
Loop
boule = u2
End Function[/B][/COLOR]