２次曲線について ５. 2次曲線の非線形写像によるイメージ（リメイク）

　math さんの問題について、「2次曲線について１．2次曲線の非線形写像によるイメージ」を Mathematica と手計算で、解きました。

　先日、その問題 (1) に関して、math さんよりMathematica のプログラムをいただきました。とても参考になったので記事にさせていただきました。（メモリーの関係でグラフィックは省略させていただきました）

(1) において、xを消去するのに Eliminateという関数が用いられています。オリジナルバージョンでは、式の形をみながら、手動でxを消去しました。自動で消去してくれるのですね。素晴らしい！

(2) の整数解に関しては、「x^2+8xy+4y^2+2x+4y+1=0の整数解！」で、WolframAlphaで、解きました。双曲型の不定方程式といわれ、この場合は無限個の解を持ちます。（どうやって解くかはわかりません）

オリジナル

before =(x-1)*y-(2*x-3)
X0=16+2(*-7*) ;(*<------可変です*)
F[x_,y_]:={(2-y)/(2x*y-2x-y), (1-x)/(2x*y-2x-y)}//Evaluate
F[x,(2*x-3)/(x-1)];
%//Together;
{X,Y}==%
Eliminate[%,x];
%[[2]]-%[[1]]//Expand
EliDE=%/.{X->x,Y->y}(*//TraditionalForm*)
pts={x,y}/.FindInstance[%==0&&-X0<x<X0&&-X0<y<X0,{x,y},Integers,170]
Pair=ListPlot[%,PlotStyle->{Magenta,PointSize[0.025]},PlotRange->All];
Fa=ContourPlot[{EliDE},{x,-X0,X0},{y,-X0,X0},PlotPoints->200,AspectRatio->Automatic,Axes->True,Contours->69];
AS=ContourPlot[{1/9 (-3 x+(6 Sqrt[3]-12) y+Sqrt[3]-3) (3 x+(12+6 Sqrt[3]) y+Sqrt[3]+3)==0,EliDE==0,before==0},{x,-X0,X0},{y,-X0,X0},AspectRatio->Automatic,Axes->True]
Show[Fa,AS,Pair]
%//Timing