Nothing
Kr_theta <- function(theta, thetaS, thetaR, n, Ks, f=0.5) {
m <- 1-(1/n)
Se <- (theta-thetaR)/(thetaS-thetaR)
a <- (1-(1-Se^(1/m))^m)^2
out <- Ks*(Se^f)*a
return(out)
}
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.