反演变换

OO={1/6,1/6};
r[t_]:={5/3 Cos[t], Sin[t]}
sol=Solve[(r[t]-OO).(x r[t]-OO)==1,x]
M=x r[t]/.sol//Flatten
ParametricPlot[{r[t],M}//Evaluate,{t,0,2 Pi}]

以OO为反演中心，以1为反演半径，r[t]是一个椭圆，M就是r[t]的反演图形。