2017-11-04 Sa Oct 消参
设切点 $ A(a, b) $，圆心 $M(x, y)$，则有 $R_M = MA = MP = R_C-MC$：
$$ \left{ \begin{aligned}
(a-3)^2 + b^2 &= 8^2 \
\sqrt{(x-a)2+(y-b)2} &= \sqrt{(x+3)2+y2} \
&= 8-\sqrt{(x-3)2+y2}
\end{aligned} \right. $$
Maxima:
solve([(a-3)^2 + b^2 = 8^2, sqrt((x-a)^2+(y-b)^2) = sqrt((x+3)^2+y^2), sqrt((x+3)^2+y^2) = 8-sqrt((x-3)^2+y^2)],[x]);
solve([(a-3)^2 + b^2 = 8^2], [a])
(%i4) solve([(a-3)^2 + b^2 = 8^2], [a]);
2 2
(%o4) [a = 3 - sqrt(64 - b ), a = sqrt(64 - b ) + 3]
solve([sqrt((x+3)^2+y^2) = 8-sqrt((x-3)^2+y^2)], [x])
解不出来……
UPD 2017-11-10 Fr 10:39PM
周一在学校的时候想了一下，由 $ (a-3)^2 + b^2 = 8^2 $ 可知$a$ $b$关系，由 $ \sqrt{(x+3)2+y2} = 8-\sqrt{(x-3)2+y2} $ 可知 $x$ $y$ 关系，这样就只剩下两个未知数，然后再带入最后一个方程就行了。
I simplfied in hand but the last step involved a polynomial with too much terms. So I decided to go to Maxima at weekend.
Maxima:
Maxima 5.25.0 http://maxima.sourceforge.net
using Lisp Clozure Common Lisp Version 1.7-r14925M (WindowsX8632)
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) solve([(a-3)^2 + b^2 &= 8^2], [a])
Maxima 5.25.0 http://maxima.sourceforge.net
using Lisp Clozure Common Lisp Version 1.7-r14925M (WindowsX8632)
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) solve([(a-3)^2 + b^2 = 8^2], [a])
;
2 2
(%o1) [a = 3 - sqrt(64 - b ), a = sqrt(64 - b ) + 3]
(%i2) a
;
(%o2) a
(%i3) a;
(%o3) a
(%i4) solve([sqrt((x-a)^2+(y-b)^2) = 8-sqrt((x-3)^2+y^2)], [y])
;
2 2 2 2 2 2
(%o4) [sqrt(y - 2 b y + x - 2 a x + b + a ) = 8 - sqrt(y + x - 6 x + 9)]
(%i5) solve([sqrt((x-a)^2+(y-b)^2) = 8-sqrt((x-3)^2+y^2)], [y]);
2 2 2 2 2 2
(%o5) [sqrt(y - 2 b y + x - 2 a x + b + a ) = 8 - sqrt(y + x - 6 x + 9)]
(%i6) solve([(x-a)^2+(y-b)^2 = 64 + (x-3)^2+y^2 - 16*sqrt((x-3)^2+y^2)], [y])
;
2 2 2 2
16 sqrt(y + x - 6 x + 9) - 2 a x + 6 x + b + a - 73
(%o6) [y = -------------------------------------------------------]
2 b
(%i7)
Seems I don't know how to use Maxima to solve equation correctly...
I gave it the wrong equation in %i4.. 'a' should not be involved. Try again:
Maxima 5.25.0 http://maxima.sourceforge.net
using Lisp Clozure Common Lisp Version 1.7-r14925M (WindowsX8632)
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) solve([sqrt((x+3)^2+y^2) = 8-sqrt((x-3)^2+y^2)], [y]);
2 2 2 2
(%o1) [sqrt(y + x + 6 x + 9) = 8 - sqrt(y + x - 6 x + 9)]
(%i2) solve([(x+3)^2+y^2 = 64 + (x-3)^2+y^2 - 16*sqrt((x-3)^2+y^2)], [y]);
Is 3 x - 16 positive, negative, or zero?
negative
;
2 2
sqrt(7) sqrt(16 - x ) sqrt(7) sqrt(16 - x )
(%o2) [y = - ---------------------, y = ---------------------]
4 4
(%i3)
Well.. let's do it step by step.
$$ \left{ \begin{aligned}
(a-3)^2 + b^2 &= 8^2 \
\sqrt{(x-a)2+(y-b)2} &= \sqrt{(x+3)2+y2} \
8-\sqrt{(x-3)2+y2} &= \sqrt{(x+3)2+y2}
\end{aligned} \right. $$
