> bAB := x - 1/2; bAB := x - 1/2 -------------------------------------------------------------------------------- > bAC := a*x+b*y-a^2/2-b^2/2 ; 2 2 bAC := a x + b y - 1/2 a - 1/2 b -------------------------------------------------------------------------------- > solve({bAB, bAC}, {x, y}); 2 2 a + b - a {y = 1/2 -----------, x = 1/2} b -------------------------------------------------------------------------------- > bBC := (1-a)*x - b*y +b^2/2 -(1-a)*(a+1)/2; 2 bBC := (1 - a) x - b y + 1/2 b - 1/2 (1 - a) (a + 1) -------------------------------------------------------------------------------- > expand(bBC); 2 2 x - a x - b y + 1/2 b - 1/2 + 1/2 a -------------------------------------------------------------------------------- > solve({bAB, bBC}, {x, y}); 2 2 a + b - a {y = 1/2 -----------, x = 1/2} b -------------------------------------------------------------------------------- > ?plot -------------------------------------------------------------------------------- > plot([0,0,1,0,1,1,0,0], style=line, color=green, axes=framed); -------------------------------------------------------------------------------- > op(1,""); CURVES([[0, 0], [1., 0], [1., 1.]]) -------------------------------------------------------------------------------- > op(1, "); 2 2 a + b - a y = 1/2 ----------- b -------------------------------------------------------------------------------- > op(2, "); 2 2 a + b - a 1/2 ----------- b -------------------------------------------------------------------------------- > lprint("); 1/2*(a^2+b^2-a)/b -------------------------------------------------------------------------------- > read(`public/Maple/thingy.mpl`); poe := proc(a,b) local Mx,My,triang,circ,rad; triang := plot([0,0,1,0,a,b,0,0],style = line,color = green); Mx := 1/2; My := 1/2*(a^2+b^2-a)/b; rad := sqrt(Mx^2+My^2); circ := plot([rad*cos(t)+Mx,rad*sin(t)+My,t = 0 .. 2*Pi],color = red) ; plot({triang,circ},axes = framed,projection = constrained) end -------------------------------------------------------------------------------- > poe(1,3); -------------------------------------------------------------------------------- > ?rand -------------------------------------------------------------------------------- > R := rand(1..100); R := proc() local t; global _seed; _seed := irem(427419669081*_seed,999999999989); t := _seed; irem(t,100)+1 end -------------------------------------------------------------------------------- > poe(R()/100, R()/100); Error, (in plot/options2d) unknown or bad argument, projection = constrained -------------------------------------------------------------------------------- > R(); 77 -------------------------------------------------------------------------------- > R(); 39 -------------------------------------------------------------------------------- > R(); 86 -------------------------------------------------------------------------------- > ?plot3d -------------------------------------------------------------------------------- > ?constrained Try one of the following topics: {convert, content, constant, conjugate, Content, const, contourplot, convergs, conformal, cond, confrac, conic, conj, contrib, connect, concur, connectivity, contribution, convexhull, concyclic, condition, confracform, contract} -------------------------------------------------------------------------------- > ?projection Try one of the following topics: {geom3d[projection], geometry[projection], plot3d[projection]} -------------------------------------------------------------------------------- > read(`public/Maple/thingy.mpl`); poe := proc(a,b) local Mx,My,triang,circ,rad; triang := plot([0,0,1,0,a,b,0,0],style = line,color = green); Mx := 1/2; My := 1/2*(a^2+b^2-a)/b; rad := sqrt(Mx^2+My^2); circ := plot([rad*cos(t)+Mx,rad*sin(t)+My,t = 0 .. 2*Pi],color = red) ; plot({triang,circ},axes = framed) end bisector := proc(a,b,c,d) (d-b)*(y-1/2*b-1/2*d)+(c-a)*(x-1/2*a-1/2*c) end -------------------------------------------------------------------------------- > bAB; x - 1/2 -------------------------------------------------------------------------------- > bisector(0,0,1,0); x - 1/2 -------------------------------------------------------------------------------- > bAC; 2 2 a x + b y - 1/2 a - 1/2 b -------------------------------------------------------------------------------- > bisector(0,0,a,b); b (y - 1/2 b) + a (x - 1/2 a) -------------------------------------------------------------------------------- > expand("); 2 2 a x + b y - 1/2 a - 1/2 b -------------------------------------------------------------------------------- > solve({bAB, bAC}, {x, y}); 2 2 a + b - a {x = 1/2, y = 1/2 -----------} b -------------------------------------------------------------------------------- > subs(",x); 1/2 -------------------------------------------------------------------------------- > subs("",y); 2 2 a + b - a 1/2 ----------- b -------------------------------------------------------------------------------- > read(`public/Maple/thingy.mpl`); poe := proc(a,b) local bAB,bAC,M,Mx,My,triang,circ,rad; triang := plot([0,0,1,0,a,b,0,0],style = line,color = green); bAB := bisector(0,0,1,0); bAC := bisector(0,0,a,b); M := solve({bAB,bAC},{y,x}); Mx := subs(M,x); My := subs(M,y); rad := sqrt(Mx^2+My^2); circ := plot([rad*cos(t)+Mx,rad*sin(t)+My,t = 0 .. 2*Pi],color = red) ; plot({triang,circ},axes = framed) end bisector := proc(a,b,c,d) (d-b)*(y-1/2*b-1/2*d)+(c-a)*(x-1/2*a-1/2*c) end -------------------------------------------------------------------------------- > poe(2,1);\ -------------------------------------------------------------------------------- >