<p>My problem is some complex, because some condition is not supplied as a given value<br>please have a look at<br><a href="http://www.lusiya.com/bbs/attachment/thumb/Mon_0906/80_3201_d2fa7bd75d28675.jpg">http://www.lusiya.com/bbs/attachment/thumb/Mon_0906/80_3201_d2fa7bd75d28675.jpg</a><br>
because I have to do futher calculation based on previous steps as <a href="http://www.lusiya.com/bbs/attachment/thumb/Mon_0906/80_3201_d342e52235049e0.jpg">http://www.lusiya.com/bbs/attachment/thumb/Mon_0906/80_3201_d342e52235049e0.jpg</a> shows, I want to get a precise value if possible, that is to say, (sqrt(2), 0), but not (1.414.0)</p>

<p>[description]<br>(clear condition)<br>A=point(0,1)<br>B=point(1,1)<br>C=point(0,0)<br>D=point(1,0)<br>AC=segment(A, C)<br>BD=segment(B, D)<br>CD=segment(C, D)<br>AB=segment(A, B)</p>
<p>I=midpoint(B, D)</p>
<p>(vague condition)<br>E=pointOn(AC)<br>F=pointOn(BD)<br>EF=segment(E, F)</p>
<p>G=mirror(A, EF)<br>G=pointOn(CD)</p>
<p>H=mirror(I, EF)<br>H=pointOn(AB)</p>
<p><br>(solve)<br>E<br>F</p>
<p>[/description]</p>
<p>Above is the description of my question, but I also hope I can write a problem like that to solve it. <br>The only lib I know is sympy, but it seems that I must write different and long code accroding to different problem</p>

<p>So is there any advanced (and lazy of cause) lib to handle my question?</p>
<p>thank you</p>