e1:x0^2+y0^2=t^2*((x1-x0)^2+(y1-y0)^2)+d2;

e2:x1^2+y1^2=(1-t)^2*((x1-x0)^2+(y1-y0)^2)+d2;

e3:solve([e1,e2],[d2,t]);

e4:factor(e3);

e5:expand((y1-y0)^2+(x1-x0)^2)=L2;

e6:subst(e5,e4);