Nothing
#############################################################################
##
## Auxiliary Functions.
##
#############################################################################
## --------------------------------------------------------------------------
arc.mean <- function(x,y) {
## ----------------------------------------------------------------------
## Compute "arc-mean" of two points.
## ----------------------------------------------------------------------
## x, y ... two points
## ----------------------------------------------------------------------
## return: arc-mean
## ----------------------------------------------------------------------
## We need x < y.
p <- sort(c(x,y))
if (p[2] < -1.e3 || p[1] > 1.e3)
## Use harmonic mean.
return (2./(1./p[1] + 1./p[2]))
a = atan(p)
if (abs(a[1]-a[2]) < 1.e-6)
## Use arithmetic mean.
return (0.5*p[1] + 0.5*p[2])
## Use "arc mean",
return (tan(0.5*a[1] + 0.5*a[2]))
} ## -- end of arc.mean() -- ##
## --------------------------------------------------------------------------
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.