Résolution système d'équations avec Maxima

Tout ce qui concerne par exemple les outils de calcul numérique, de calcul formel ou de géométrie.
[participation réservée aux membres inscrits]
Règles du forum
Merci de soigner la rédaction de vos messages et de consulter ce sujet avant de poster. Pensez également à utiliser la fonction recherche du forum.
polio

Résolution système d'équations avec Maxima

Message non lu par polio »

Bonjour,

Nouvel utilisateur de Maxima, je cherche à résoudre un système d'équations. Mon code est le suivant mais ne fonctionne pas :

Code : Tout sélectionner

/* 1er birapport */
PAB:matrix([xp,xa,xb],[yp,ya,yb],[1,1,1]);
PAC:matrix([xp,xa,xc],[yp,ya,yc],[1,1,1]);
PCB:matrix([xp,xc,xb],[yp,yc,yb],[1,1,1]);
PCD:matrix([xp,xc,xd],[yp,yc,yd],[1,1,1]);
br1:(determinant(PAB)/determinant(PAC))/(determinant(PCB)/determinant(PCD));

/* 2nd birapport */
PDC:matrix([xp,xd,xc],[yp,yd,yc],[1,1,1]);
PDB:matrix([xp,xd,xb],[yp,yd,yb],[1,1,1]);
PAB:matrix([xp,xa,xb],[yp,ya,yb],[1,1,1]);
PAC:matrix([xp,xa,xc],[yp,ya,yc],[1,1,1]);
br2:(determinant(PDC)/determinant(PDB))/(determinant(PAB)/determinant(PAC));

/* Résolution du système d'équations */
solve([br1=i,br2=j],[xp,yp]);
J'ai "simplifié" l'équation sous cette forme :

Code : Tout sélectionner

solve([((- xa * (yp - yb) + xb * (yp - ya) + xp * (ya - yb)) * (- xc * (yp - yd) + xd * (yp - yc) + xp * (yc - yd))) / ((- xa * (yp - yc) + xc * (yp - ya) + xp * (ya - yc)) * (xb * (yp - yc) - xc * (yp - yb) + xp * (yc - yb))) = i,((- xa * (yp - yc) + xc * (yp - ya) + xp * (ya - yc)) * (xc * (yp - yd) - xd * (yp - yc) + xp * (yd - yc))) / ((- xa * (yp - yb) + xb * (yp - ya) + xp * (ya - yb)) * (xb * (yp - yd) - xd * (yp - yb) + xp * (yd - yb))) = j], [xp, yp]);
... Mais elle ne fonctionne pas non plus.

Merci pour l'aide.