wi=0 wii=NA H=2 L=30 v=1 wi=rep(0,83) plot(wi,type="l",xlim=c(0,90),ylim=c(-5,5)) abline(v=83,lwd=3) title("Calm conditions") for(i in 1:10) { date_time=Sys.time() while((as.numeric(Sys.time()) - as.numeric(date_time))<0.1){} } for(i in 1:83) { wi[1:i]=H/2*cos(2*pi*(seq(1,i)/L-i/v/T)-i/7) plot(wi,type="l",xlim=c(0,90),ylim=c(-5,5)) title("Incoming wave") abline(v=83,lwd=3) date_time=Sys.time() while((as.numeric(Sys.time()) - as.numeric(date_time))<0.1){} } for(i in 1:83) { wi=H/2*cos(2*pi*(seq(7.5,90)/L-i/v/T)-i/7) wii[(84-i):83]=H/2*cos(2*pi*(seq((90-i),89)/L-i/v/T)+i/7) plot(wi,type="l",xlim=c(0,90),ylim=c(-5,5)) title("Reflected wave") lines(wii,col="red") lines(wi+wii,col="blue",lwd=2) abline(v=83,lwd=3) date_time=Sys.time() while((as.numeric(Sys.time()) - as.numeric(date_time))<0.1){} } #for(i in 1:198) for(i in 1:176) #for(i in 1:99) { wi=H/2*cos(2*pi*(seq(7.5,90)/L-i/v/T)-i/7) wii=H/2*cos(2*pi*(seq(7.5,90)/L-i/v/T)+i/7) plot(wi,type="l",xlim=c(0,90),ylim=c(-5,5)) title("Standing wave") lines(wii,col="red") lines(wi+wii,col="blue",lwd=2) abline(v=83,lwd=3) date_time=Sys.time() while((as.numeric(Sys.time()) - as.numeric(date_time))<0.1){} }