bonyolult egyenletrendszer egyszerűsítése, computer algebra system

Sziasztok,

hobbiprojektem magját ma fejeztem be, egy bonyolult egyenletrendszer. a projektet egyelőre yacas-ban írom, és amint elkészültem papíron az egyenletekkel, ráengedtem a Solve() függvényt, ami már órák óra megy, és kezdek kétségbe esni, mert annyira nem lehet bonyolult.

a fenti művelettől azt vártam, hogy gyorsan dob egy eredményt. nem tette, úgyhogy felmerültek bennem a következő kérdések:

egyik kérdésem az, hogy szerintetek az alábbi egyenleteket hogyan lehetne átírni valahogyan egy Java metódusba? mondjuk, a metódus megkapná az A*, B*, cu*, ra* változók értékeit, és visszatérne a C*, D* változók értékeivel. azt tudom, hogy egyenletrendszermegoldó képessége nincs annak a nyelvnek, de talán át lehetne írni valahogyan..

szerintetek hogyan lehetne a fentieket egyszerűbb alakra hozni?

egy másik kérdésem, hogy ti melyik cas rendszereket használtátok, mik a tapasztalataitok. melyik a legnépszerűbb? érdemes cas rendszerre költeni (maxima, maple), vagy az opensource alternatívák vannak legalább olyan jók? én a yacas-sal és a maxima-val ismerkedtem. a yacas volt szimpatikus, de az oo lehetősége nekem mindkettőből nagyon hiányzik.

legvégül a kérdéses borzadmány:

get2Ddistance(x1,y1,x2,y2):=Sqrt((x1-x2)^2+(y1-y2)^2);

(Cx-Ax)^2+(Cy-Ay)^2+(Cz-Az)^2=(currentWidth*ratioWtoAC)^2 , és
(Cx-Bx)^2+(Cy-By)^2+(Cz-Bz)^2=(currentWidth*ratioWtoBC)^2 , továbbá,
(Dx-Cx)^2+(Dy-Cy)^2+(Dz-Cz)^2=(currentWidth*ratioWtoCD)^2 , ahol, 
(Dx-Bx)^2+(Dy-By)^2+(Dz-Bz)^2=(currentWidth*ratioWtoBD)^2 , ráadásul
get2Ddistance(Cx,Cy,0,0)/get2Ddistance(Dx,Dy,0,0)=get2Ddistance(Ax,Ay,0,0)/get2Ddistance(Bx,By,0,0) , de
Cz<0 , és
Dz>0

válaszaitokat előre is köszönöm,
zp