R画图

 

> plot(x,z, col="green")
> lines(x, x*0+10, col="red", lwd=5)
> title("limit=10")
> tmp1 <- locator(1)
> text(tmp1, "actual flow")