R/get.d.R

In egg: Distances and areas on the surface of an egg

get.d = function(z.hat, egg) {
  # Get the value of d for a given value of z.hat
  d = (egg$a*exp(((-z.hat^2)/(2*egg$b^2)) + ((egg$c*z.hat)/(egg$b^2)) - ((egg$c^2)/(2*egg$b^2)))*sqrt(1 - z.hat)*sqrt(z.hat))/(pi*egg$b)
  return(d)
}