i2pi · June 22, 2021 06:27
diff --git a/gistfile1.md b/gistfile1.md
diff --git a/love_simple.R b/love_simple.R
 heart <- function (w) {
  x <- 16*sin(w)^3  
  y <- 13*cos(w)-5*cos(2*w)-2*cos(3*w)-cos(4*w)
  return (cbind (x,y)/32)
 }

 wave <- function (a, R, p, w) {
  h <- heart(a)
  x <- (0.5 + 0.5*R*sin(p+w*a))*h[,1]
  y <- (0.5 + 0.5*R*sin(p+w*a))*h[,2]
  return (cbind(x,y))
 }

 par(bg=rgb(148/255,131/255,171/255))
 par(fg='white')

 t <- seq(0, 2*pi, length.out=5000)
 M<-0.35
 plot(c(-1, 1)*M, c(-1, 1)*M, pch='', bty='n', axes='F', xlab='', ylab='')
 for (i in 0:7) lines(wave(t,0.3,2*pi*i/8,2), lwd=3)
	heart <- function (w) {
	x <- 16*sin(w)^3
	y <- 13cos(w)-5cos(2w)-2cos(3w)-cos(4w)
	return (cbind (x,y)/32)
	}

	wave <- function (a, R, p, w) {
	h <- heart(a)
	x <- (0.5 + 0.5Rsin(p+wa))h[,1]
	y <- (0.5 + 0.5Rsin(p+wa))h[,2]
	return (cbind(x,y))
	}

	par(bg=rgb(148/255,131/255,171/255))
	par(fg='white')

	t <- seq(0, 2*pi, length.out=5000)
	M<-0.35
	plot(c(-1, 1)M, c(-1, 1)M, pch='', bty='n', axes='F', xlab='', ylab='')
	for (i in 0:7) lines(wave(t,0.3,2pii/8,2), lwd=3)