Pour ce genre de choses, je ne vois pas pourquoi tu n'approximes pas ton ellipse pas des segments de droites. Avec le truc super basique ci dessous, j'obtiens 3.14142 comme approximation de pi, ce qui n'est pas trop mal. En faisant un truc plus évolué avec un pas adaptatif, on peut très bien obtenir quelque chose de suffisamment précis et de rapide.
u=5cm;
beginfig(1)
path elipse;
pair a;
elipse=(u,0)..
for i=1 upto 99:
((u,0) rotated (3.6*i))..
endfor cycle;
s=0;
l=.01;
for i=l step l until length(elipse):
a:=(point i of elipse)-(point (i-l) of elipse);
s:=s+((xpart a)++(ypart a));
draw point i of elipse;
endfor;
show s/(2u);
endfig;
end
[^] # Re: Et par rapport à pgf/Tikz et autres ?
Posté par fmaz fmaz . En réponse au journal METAPOST 1.0 est sorti. Évalué à 3.
u=5cm;
beginfig(1)
path elipse;
pair a;
elipse=(u,0)..
for i=1 upto 99:
((u,0) rotated (3.6*i))..
endfor cycle;
s=0;
l=.01;
for i=l step l until length(elipse):
a:=(point i of elipse)-(point (i-l) of elipse);
s:=s+((xpart a)++(ypart a));
draw point i of elipse;
endfor;
show s/(2u);
endfig;
end