with(plots):

heart_polar := (t) -> [t, 2 - 2*sin(t) + sin(t)*sqrt(abs(cos(t))) / (sin(t) + 1.4)]:


display(
    plot(heart_polar(t), t=0..2*Pi, coords=polar, color=red, thickness=2),  # Heart shape
    scaling=constrained,  # Ensure equal scaling
    title="I Love Quantum Computing"
);