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" );