Computations

with(linalg):

>

> A:=[5*cos(theta),5*sin(theta)];

> B:=[10+5*cos(theta),5*sin(theta)];

> C:=[5*cos(theta)+2,5*sin(theta)];

> DD:=[6*cos(theta),6*sin(theta)];

> E:=[10,0];

> P:=[9*cos(theta),sin(theta)];

> Q:=C+[alpha,beta];

> eq1:=dotprod(Q-C,Q-C,orthogonal)=1;

> eq2:=dotprod(Q-DD,Q-DD,orthogonal)=4;

> Q:=subs(solve({eq1,eq2},{alpha,beta})[2],Q);

> R:=subs(solve({eq1,eq2},{alpha,beta})[2],C+[4*alpha,4*beta]);

> theta:='theta';

> with(plottools): with(plots):

> OA:=scalarmul(A,t):

> LOA:=[OA[1],OA[2],t=0..1]:

> OE:=scalarmul(E,t):

> LOE:=[OE[1],OE[2],t=0..1]:

> EB:=matadd(scalarmul(E,t),scalarmul(B,1-t)):

> LEB:=[EB[1],EB[2],t=0..1]:

> AB:=matadd(scalarmul(A,t),scalarmul(B,1-t)):

> LAB:=[AB[1],AB[2],t=0..1]:

> AD:=matadd(scalarmul(A,t),scalarmul(DD,1-t)):

> LAD:=[AD[1],AD[2],t=0..1]:

> AC:=matadd(scalarmul(A,t),scalarmul(C,1-t)):

> LAC:=[AC[1],AC[2],t=0..1]:

> DQ:=matadd(scalarmul(DD,t),scalarmul(Q,1-t)):

> LDQ:=[DQ[1],DQ[2],t=0..1]:

> CQ:=matadd(scalarmul(C,t),scalarmul(Q,1-t)):

> LCQ:=[CQ[1],CQ[2],t=0..1]:

> CR:=matadd(scalarmul(C,t),scalarmul(R,1-t)):

> LCR:=[CR[1],CR[2],t=0..1]:

> AP:=matadd(scalarmul(A,t),scalarmul(P,1-t)):

> LAP:=[AP[1],AP[2],t=0..1]:

> PR:=matadd(scalarmul(P,t),scalarmul(R,1-t)):

> LPR:=[PR[1],PR[2],t=0..1]: