Nothing
sdfarma <- function(n, phi=numeric(0), theta=numeric(0)) {
lams <- 2*pi*seq(from=1/n, to=1/2, by=1/n)
nf <- length(lams)
if (length(theta)>0) {
a <- outer(lams, 1:length(theta))
C <- cbind(1, cos(a))%*% c(1, -theta)
S <- sin(a) %*% theta
} else {
C <- 1
S <- 0
}
num <- as.vector(C*C + S*S)/(2*pi)
if (length(phi)>0) {
a <- outer(lams, 1:length(phi))
C <- cbind(1, cos(a))%*% c(1, -phi)
S <- sin(a) %*% phi
} else {
C <- 1
S <- 0
}
den <- as.vector(C*C+S*S)
num/den
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.