dat<- data.frame(t=seq(0, 2*pi, by=0.1) )
xhrt <- function(t) 16*sin(t)^3
yhrt <- function(t) 5*(13*cos(t)-5*cos(2*t)-2*cos(3*t)-cos(4*t))
dat$y=yhrt(dat$t)
dat$x=xhrt(dat$t)
with(dat, plot(x,y, type="l"))
with(dat, polygon(x,y,col="#E74A4A"))

x <- seq(-2.5, 2.5, by=0.000001)
plot(x, (sqrt(cos(x))*cos(200*x)+sqrt(abs(x))-0.7)*(4-x*x)^0.01, ylab="y", col="#E74A4A", "l", lwd=2)
## Warning in sqrt(cos(x)): NaNs produced