x = var('x')
outside = (sqrt(cos(x))*cos(100*x)+sqrt(abs(x))-0.4)*(4-x*x)^0.1
inside = (sqrt(cos(x))*cos(100*x)+sqrt(abs(x))-0.4)*(4-x*x)^0.1
heart = [plot(outside, xmin=-1.6, xmax=1.6, ymin = -2, ymax = 2, gridlines = True ),
plot(inside, xmin=-1.6, xmax=1.6, ymin = -2, ymax = 2, gridlines = True )]
p = sum(heart)
p.show()