days=seq(1,5,length=100)
y=numeric(100)
y[1:50]=2*days[1:50]
y[51:100]=rep(2*days[51],50)
z=rnorm(100,0,.15)
y=y+z
plot(days,y)