R代码唯美 — 一鸣惊人

代码如下：
xrange = c(-15, 15) ;
yrange = c(0, 16) ;
plot(0, xlim = xrange, ylim = yrange, type = “n”) ;
yr = seq(yrange[1], yrange[2], len = 50) ;
offsetFn = function(y) {     2 * sin(0 + y/3) }  ;
offset = offsetFn(yr)  ;
  leftE = function(y) {     -10 – offsetFn(y) }  ;
  rightE = function(y) {     10 + offsetFn(y) } ;
  xp = c(leftE(yr), rev(rightE(yr)))   ;
   yp = c(yr, rev(yr)) ;
    polygon(xp, yp, col = “#ffeecc”, border = NA) ;
     h = 9 ;
     xt = seq(0, rightE(h), len = 100);
     yt = log(1 + log(1 + log(xt + 1))) ;
     yt = yt – min(yt) ;
     yt = h * yt/max(yt)  ;
     x = c(leftE(h), rightE(h), rev(xt), -xt) ;
     y = c(h, h, rev(yt), yt) ;
     polygon(x, y, col = “red”, border = NA) ;
 
去掉坐标后的制品：
xrange = c(-15, 15) ;
yrange = c(0, 16) ;
plot(xrange, yrange, type=”n”, ann=FALSE, axes=FALSE) ;
yr = seq(yrange[1], yrange[2], len = 50) ;
offsetFn = function(y) {     2 * sin(0 + y/3) }  ;
offset = offsetFn(yr)  ;
  leftE = function(y) {     -10 – offsetFn(y) }  ;
  rightE = function(y) {     10 + offsetFn(y) } ;
  xp = c(leftE(yr), rev(rightE(yr)))   ;
   yp = c(yr, rev(yr)) ;
    polygon(xp, yp, col = “#ffeecc”, border = NA) ;
     h = 9 ;
     xt = seq(0, rightE(h), len = 100);
     yt = log(1 + log(1 + log(xt + 1))) ;
     yt = yt – min(yt) ;
     yt = h * yt/max(yt)  ;
     x = c(leftE(h), rightE(h), rev(xt), -xt) ;
     y = c(h, h, rev(yt), yt) ;
     polygon(x, y, col = “red”, border = NA) ;