Comment intégrer (numériquement) ?

[inviteaea18ef8]

Bonjour, j'ai hésité à poster dans le forum informatique mais je crois que ma question est un peu plus générale.

Je cherche à réaliser une intégrale. Je ne cracherais pas sur la solution littérale, mais ça semble illusoire donc une approximation numérique m'irait très bien.

Mon intégrale c'est ça :


J'utilise le logiciel de calcul formel Maxima, mais quand j'entre ceci :

%lambda : 2.2;
ri : 2.1;
hi : 3.056;
he : 17.778;
re(h) := -0.033*h + 3.652;
Y(h) := (2*%pi)/(1/(re(h)*he)+log(re(h)/ri)/%lambda+1/(ri*hi));
YI : integrate(Y(h),h,2,3);

Il me renvoit l'intégrale non calculée. Vu que ça marche avec des intégrale plus simples, j'en déduis qu'il y a une simplification que je ne sais pas faire, ou qu'il n'y arrive pas.

Quelqu'un sait utiliser mieux Maxima?
Quelqu'un peut me dire si sur son programme de prédilection cette intégrale est calculée? (ça ce serait trop sympa).

Merci
-----

[invite4793db90]

Salut,

essaie l'une des commandes suivantes :

romberg(Y(h), h, 2, 3);
quad_qags(Y(h), h, 2, 3);

Cordialement.

[invite1e1a1a86]

Mapple arrive a calculer cette intégrale après un "evalf" donc Maxima aussi

sauf erreur: 22.74042801 (si c'est un log à base dix) et 15.22436263 (si c'est ln népérien)

[inviteaea18ef8]

Merci pour vos réponses.

Mais j'ai déjà essayé, ça ne semble pas marcher :
ça me renvoit romberg(...tout le tintoin...) sans afficher le résultat.

J'ai pensé aussi qu'il y a avait un moyen de le forcer à évaluer (et j'espère toujours un peu), mais j'essaie en vain "ev(%)", "float(%)"...

[A voir en vidéo sur Futura]

[invite4793db90]

Salut,

c'est surprenant, car maxima donne chez moi une sortie numérique (15.2..). Disposes-tu de la dernière version ?

Cordialement.

[inviteaea18ef8]

Vous n'allez pas me croire, mais après réinstallation ça marche! Exactement en lançant le même script. Je ne l'avais pas téléchargée depuis le sourceforge officiel, mais c'était en apparence la même version! (pour preuve : http://forum.framasoft.org/viewtopic.php?f=17&t=34188 )

Bon merci à vous en tous cas! Je suis bien content de pouvoir continuer avec un logiciel libre, même s'il paie pas de mine.