2021，吾以笛卡尔之心开始

“

R version 4.03

谨以此献给我的爱人，感谢她对我的支持与帮助######

Date:2021.1.1

Author:杨博士聊生信

”

笛卡尔是法国著名的哲学家，数学家和科学家。笛卡尔心形函数来自于一个美丽的爱情传说。相传笛卡尔落魄在瑞典之时，被瑞典王国一位喜欢数学的公主所欣赏，把他请进了王宫当老师，他们在研究数学的世界中日久生情，国王知道后，大发雷霆要处死笛卡尔。公主百般求情，笛卡尔被驱逐回法国，不久得了黑死病，他给公主写了13封信，公主收到了，信上面有一个方程：r=a（1-sinθ），画出来的就是一个心形线。他用一生至爱的方程告诉他爱她，这也是被人们所认为的最美的表达爱情的公式，被世人称为笛卡尔心形线。
在这借用笛卡尔的心形函数，感谢我媳妇一直以来对我的支持。

以下是代码，喜欢的小伙伴可以试一下

x<-seq(-1.1,1.1,length = 3000)
rep<-30
y<-matrix(,3000,rep) #构建矩阵
yx<-matrix(,3000,rep)
xx<-matrix(,3000,rep)
c<-seq(0,1,length=rep) ##设置立体效果
for (i in 1:rep)
for (p in 1:3000){
  y[p,i]<-(1-(x[p])^2)^(1/2)+((x[p])^2)^(1/4)+c[i]
  yx[p,i]<--(1-(x[p])^2)^(1/2)+((x[p])^2)^(1/4)+c[i]
  xx[,i]<-x}
for (q in 1:100){
par(bg="black")
plot(main = "You are My heart", col.main = q, c(rev(xx),xx),c(rev(y),yx), type = 'l' ,xlim = c(-1,1), xlab = "I love three things in the world, 
     Sun for day, moon for night and you forever",col.lab = q,  ylab = "",  axes = F,col=q)   #我们可根据需要改动main和xlab后面对应的文字
Sys.sleep(1)} #Sys.sleep(1)动态效果

ref:
1.r=a(1−sinθ)'>r=a(1−sinθ) r=a(1−sinθ) 笛卡尔心形函数
2.https://blog.csdn.net/hhtnan/article/details/79282218：绘制动态心形图案：：R语言绘制心形图
3.https://www.bilibili.com/read/cv7166339/