t <- seq(1,60,.01) y0 <- 100 + 0*t y1 <- y0 + 5 * cos(2*pi*1*t/12) - 3 * sin(2*pi*1*t/12) y2 <- y0 + 5 * cos(2*pi*2*t/12) - 3 * sin(2*pi*2*t/12) y3 <- y0 - 5 * cos(2*pi*3*t/12) + 3 * sin(2*pi*3*t/12) y4 <- y0 + 5 * cos(2*pi*4*t/12) - 3 * sin(2*pi*4*t/12) y5 <- y0 - 5 * cos(2*pi*5*t/12) + 3 * sin(2*pi*5*t/12) win.graph(height=5.5, width=7.0) par(mfrow=c(2,3)) plot(y0 ~ t, type="l", col="purple", ylim=c(50,150)) plot(y1 ~ t, type="l", col="red", ylim=c(50,150)) plot(y2 ~ t, type="l", col="green", ylim=c(50,150)) plot(y3 ~ t, type="l", col="blue", ylim=c(50,150)) plot(y4 ~ t, type="l", col="brown", ylim=c(50,150)) plot(y5 ~ t, type="l", col="orange", ylim=c(50,150)) y0 <- 100 + 0*t y1 <- y0 + 5 * cos(2*pi*1*t/12) - 3 * sin(2*pi*1*t/12) y2 <- y1 + 5 * cos(2*pi*2*t/12) - 3 * sin(2*pi*2*t/12) y3 <- y2 - 5 * cos(2*pi*3*t/12) + 3 * sin(2*pi*3*t/12) y4 <- y3 + 5 * cos(2*pi*4*t/12) - 3 * sin(2*pi*4*t/12) y5 <- y4 - 5 * cos(2*pi*5*t/12) + 3 * sin(2*pi*5*t/12) win.graph(height=5.5, width=7.0) par(mfrow=c(2,3)) plot(y0 ~ t, type="l", col="purple", ylim=c(50,150)) plot(y1 ~ t, type="l", col="red", ylim=c(50,150)) plot(y2 ~ t, type="l", col="green", ylim=c(50,150)) plot(y3 ~ t, type="l", col="blue", ylim=c(50,150)) plot(y4 ~ t, type="l", col="brown", ylim=c(50,150)) plot(y5 ~ t, type="l", col="orange", ylim=c(50,150))