Homemade leaf

I wrote a little program in R (a programming language for statistics) to produce a homemade leaf (for nature enthusiasts in quarantine):

 

postscript("Homemade_Leaf.ps", width = 10, height = 10)

 

phi=seq(-pi,pi,pi/200)

 

### ???

r15=0.25*(1+cos(phi))*(1+cos(5*phi))+0.1*(1+cos(phi))*(1+cos(15*phi))+(1-0.2*cos(phi))*(1+cos(phi))

phi15=phi+0.025*((1+cos(phi))*(-sin(5*phi)-2*sin(15*phi)/5)+(1-cos(phi))*(-sin(phi)))

 

r15=r15/max(r15)

plot(r15*sin(phi15),r15*cos(phi15),type='l',xlim=c(-1,1),ylim=c(-0.3,1))

MAX=8

for(i in 1:MAX){

ri=r15*i/MAX

lines(ri*sin(phi15),ri*cos(phi15),col='gray')

}

for(j in 1:length(r15)){

vx=rep(0,MAX+1)

vy=rep(0,MAX+1)

for(i in 1:MAX){

rji=r15[j]*i/MAX

vx[i]=rji*sin(phi15[j])

vy[i]=rji*cos(phi15[j])

}

lines(vx,vy,col='darkgreen')

}

 

dev.off()

 

 

Done