Bonjour,

 

Je ne crois pas que le problème vienne de %pi ou %phi. La commande solve échoue 
généralement à résoudre des équations ou système d’équations
dont les inconnues ce trouvent sous des exposant rationnel.

 

Il existe en revanche un module en cours de développement « to_poly_solve ».qui 
doit être en mesure de résoudre votre équation.

Pour l’utiliser faite:

 

load(to_poly_solve);

%solve([c(x) = f2(x)], x);

 

Pour les méthodes numériques il y a entre autre:

find_root, bf_find_root, allroots, newton et le module mnewton.

 

 

Laurent.

 

-----Message d'origine-----
De : Patrick Gelin [mailto:[email protected]] 
Envoyé : jeudi 19 novembre 2015 17:31
À : [email protected]
Objet : [Maxima-lang-fr] Problème avec la définition des constantes %pi et %phi 
...

 

Bonjour,

C'est une question de débutant sous Maxima. J'ai deux fonctions, définies 
ci-dessous :
f2(x):=sqrt(%pi^(3/2)*sqrt(16*%phi^(3/2)*x+%phi*%pi+%pi)−32*x^2+8*sqrt(%phi)*%pi*x+%phi*%pi^2)/2^(5/2);
c(x):=sqrt(−16*x^2+8*sqrt(%phi)*%pi*x−%phi*%pi^2+16*%phi)/4;

Voici le graphique, j'espère que la mailing liste le laisse en place... On peut 
voir un point d'intersection autour de x=1 :





Je cherche la valeur précise de x, absice du point d'intersection c(x)=f2(x). 
J'ai essayé de trouver la solution algébrique avec la commande solve :
solve([c(x)=f2(x)],[x]);
[sqrt(%pi^(3/2)*sqrt(16*%phi^(3/2)*x+(%phi+1)*%pi)−32*x^2+8*sqrt(%phi)*%pi*x+%phi*%pi^2)
 = sqrt(2)*sqrt(−16*x^2+8*sqrt(%phi)*%pi*x−%phi*%pi^2+16*%phi)]

On voit que la commande pose le problème sans le résoudre... J'ai l'impression 
que maxima considère les constantes %pi et %phi non pas comme des constantes 
mais comme des variables supplémentaires de l'équation... Un message d'erreur 
sur un autre test me l'a suggéré. Alors 3 variables pour une équation, il ne 
peut donc pas trouver la solution... Comment lui faire reconnaitre que %pi et 
%phi ne sont pas des variables mais bien des constantes définies ? J'ai essayé 
de supprimer le % devant pi et phi, cela ne change pas le problème.

Sinon existerait t-il une autre commande avec un autre algorithme qui pourrait 
trouver la valeur numérique de la solution,  et non chercher son expression 
algébrique ?

Merci pur votre aide.
Patrick




 
------------------------------------------------------------------------------
_______________________________________________
Maxima-lang-fr mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/maxima-lang-fr

Répondre à