The bifolium curve plotted using Sage software using the following code.


var('x y')
f = (x^2+y^2)^2-x^2*y
pl = implicit_plot(f, (x,-.4,.4), (y,-.05,.3), aspect_ratio=1.1,cmap='hsv')
pl.axes_labels(['$x$','$y$'])
pl.show()
Own work
Krishnavedala