目标：已有要在yy1,yy2两个函数之间的区域填上灰色（Confidence band ）

请问如何做？

ps：我的问题代码
plot(x,yy1,col="blue",type = "p")  #,xlim = c(0,2)
points(x,yy2)
#polygon(c(x,rev(x)),c(x^3+q,rev(x^3-q)),col = "gray")
#polygon(c(x),c(x^3+q,rev(x)^3-q),col = "gray")
polygon(c(xx),c(yy1,yy2),col = "gray")