在我的示例中,我想模拟椭圆轴:
a <- x # semi-major axis
e <- 1/sqrt(2) # eccentricity
b <- a * sqrt(1 - e^2) # semi-minor axis
c <- 6.3 # ellipse area
当我具有平均偏心率(e)和不同的短轴(b)和椭圆区域(c)时,我的代码需要估算999个重复的平均主轴(a)。任何想法? 谢谢
答案 0 :(得分:1)
a = x
b = x * sqrt(1 - e^2)
area = pi * x * x * sqrt(1 - e^2) = pi * x^2 * sqrt(1 - e^2)
x = sqrt( area / (pi * sqrt(1 - e^2)) );