もともとは，
次の方程式のグラフの概形をかけ。
というものだったけど，描いてみたら，片思いハートだったので，両思いにして，さらにたくさん描いてみた。
fxy = function(x, y) x^2+(y-abs(x)^(2/3))^2-1
x = seq(-2.5, 2.5, length=1000)
y = seq(-2.5, 3.5, length=1000)
z = matrix(0, 1000, 1000)
for (i in 1:1000) {
for (j in 1:1000) {
z[i,j] = fxy(x[i], y[j])
}
}
contour(x, y, z, drawlabels=FALSE, levels=c(-9:1/10, 0:5), asp=1, col="red")
