有，你看这个可以么：
xrange=c(-15,15)<br />
yrange=c(0,16)<br />
plot(0,xlim=xrange,ylim=yrange,type='n')<br />
<br />
yr=seq(yrange[1],yrange[2],len=50)<br />
offsetFn=function(y){2*sin(0+y/3)}<br />
offset=offsetFn(yr)<br />
leftE = function(y){-10-offsetFn(y)}<br />
rightE = function(y){10+offsetFn(y)}<br />
<br />
<br />
xp=c(leftE(yr),rev(rightE(yr)))<br />
yp=c(yr,rev(yr))<br />
polygon(xp,yp,col="#ffeecc",border=NA)<br />
<br />
h=9<br />
<br />
xt=seq(0,rightE(h),len=100)<br />
yt=log(1+log(1+log(xt+1)))<br />
yt=yt-min(yt)<br />
yt=h*yt/max(yt)<br />
<br />
x=c(leftE(h),rightE(h),rev(xt),-xt)<br />
y=c(h,h,rev(yt),yt)<br />
polygon(x,y,col="red",border=NA)
