3839.sagews

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()

x = var('x')

outside = (sqrt(cos(x))*cos(200*x)+sqrt(abs(x))-0.4)*(4-x*x)^0.1
inside =  (sqrt(cos(x))*cos(200*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()

x = var('x')

outside = (sqrt(cos(x))*cos(300*x)+sqrt(abs(x))-0.4)*(4-x*x)^0.1
inside =  (sqrt(cos(x))*cos(300*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()