rm(list=ls())
require(animation)
require(rv)

par(bg = "white")
plot(0:40,seq(-20,20,by=1),type="n",xlab="Ship's Y-Axis",ylab="Ship's Z-Axis")

q<-rnorm(30,mean=0,sd=7)
qWave<-c(mean(q))
r<-0

for (n in seq(1,40,by=1)) {

r<-c(r,n)
v<-rnorm(1,mean=0,sd=7)
q<-c(q,v)
qWave<-c(qWave,mean(q[seq((1+n),(29+n),by=1)]))

}

plot(0:40,seq(-20,20,by=1),type="n",xlab="Ship's Y-Axis",ylab="Ship's Z-Axis")
lines(r,qWave,col="blue")


## set animation options
oopt <- ani.options(interval = 0.1)

for (m in seq(41,200,by=1)) {

plot((max(0,m-40):max(40,m)),seq(-20,20,by=1),type="n",xlab="Ship's Y-Axis",ylab="Ship's Z-Axis")

r<-c(r,m)
v<-rnorm(1,mean=0,sd=7)
q<-c(q,v)
qWave<-c(qWave,mean(q[seq((1+m),(29+m),by=1)]))
lines(r,qWave,col="blue")
##points(m,v)

    ani.pause()

}