int et precision

[frederic han]

salut,

j'ai l'impression qu'a partir de 16 chiffres il y a un probleme avec l'integration qui a l'ai de tourner sans fin:

evalf(int(1/sqrt(x*(x+4)*(x+1)),x,0,+infinity),15) est OK

evalf(int(1/sqrt(x*(x+4)*(x+1)),x,0,+infinity),16) a l'air de boucler sans message d'erreur.

Fred

[parisse]

en effet, ca vient de la borne infinie, xcas calcule alors
1/sqrt(tan(x)^3+5*tan(x)^2+4*tan(x))*(1+tan(x)^2)
pour x=0..pi/2, avec une divergence en pi/2, donc appel de l'acceleration de Romberg en utilisant le point milieu (et non les trapezes) ... qui n'est implemente qu'en double precision. Je vais commencer par convertir correctement les bornes d'integration en double et ajouter un warning.

[parisse]

voila, c'est corrige, je viens de voir qu'il y avait deja un warning, maintenant ca ne bouclera plus mais ce ne sera pas pour autant plus precis (en tout cas pour le moment).

[frederic han]

OK merci.

Y a t'il des options pour imposer une methode de calcul du genre trapeze, romberg...

a+

[parisse]

oui, partiellement, avec la commande plotarea d'une part et romberg d'autre part (sachant que evalf(int()) appelle romberg).