カッシーニ曲線 x<-y<-seq(from=-2,to=2,by=0.01) xy<-expand.grid(x,y) a<-4 b<-4 z<-(xy[,1]^2+xy[,2]^2)^2-2*a^2*(xy[,1]^2-xy[,2]^2)+a^4 image(x,y,matrix(z,length(x),length(y)))

こちら(ikuro先生のページは本当に素晴らしい) ２定点からの距離の和が一定となる点の軌跡は楕円(２定点が一致すれば、円)、差が一定の点の軌跡は双曲線(２定点が一致すれば、バツ印)、商が一定の点は円（アポロニウスの円）（２定点が一致すれば、商はつね…