Merci de votre réponse rapide

laurent couraud wrote:
> Je suppose que certain doublons ne sont pas identifié par le listify(setify())
> Si les expressions ne sont pas simplifiées.
>   
J'ai en effet peut-être des doublons, mais ce qui est étrange, c'est 
plus loin que le problème se situe. En fait, dans la boucle for ensuite, 
je prends les droites une à une, et je supprime les points qui sont du 
mauvais côté de la droite. L'algo sans le ratsimp me donne le bon 
résultat, avec le bon nombre de points, et aucun doublons, alors que 
lorsque je mets le ratsimp, j'ai des points qui sont supprimés alors 
qu'ils ne le devraient pas.
> Vous pouvez certainement le vérifié en faisant la différence des solutions
>  avec et sans ratsimp()
>   
J'essaierai, merci du conseil
> Quand à une méthode plus simple pour obtenir le même résultat, j'avoue que je 
> n'ai
>  pas très bien compris ce que vous cherchiez à faire.
>   
En entrée, j'ai un domaine défini par des inéquations du type 0 <= 
ax+by+c <= 1, je cherche à calculer le polygone correspondant à la zone 
où x et y vérifient l'ensemble de ces inéquations.
> En tous les cas vous pourriez utiliser la fonction solve pour rechercher 
> l'intersection de 2
> droite.
>
> intersectDroite(droite1,droite2):=solve([equationDroite(droite1),equationDroite(droite2)],[x,y]
> )$
>   
C'est ce que je faisais au début, mais il me donnait en résultat une 
liste de la forme [x=truc,y=machin], et il me fallait [truc,machin], 
donc je devait rajouter map(rhs,la_solution), et de plus, pour gérer les 
droites parallèles, je devais faire en sorte de supprimer l'erreur pour 
les équations inconsistantes. Au final, le programme se révélait plus 
court de cette façon...
Merci beaucoup pour vos conseils, en tout cas

Vincent

-------------------------------------------------------------------------
This SF.net email is sponsored by: Microsoft
Defy all challenges. Microsoft(R) Visual Studio 2008.
http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/
_______________________________________________
Maxima-lang-fr mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/maxima-lang-fr

Répondre à