>>590
面積を S とする
S = 4πa∫[0,π/2]√{a^2*cos^2(θ)+c^2*sin^2(θ)}sin(θ)dθ

r = a/c, s = √|r^2-1| とすると
S = 2π{1 + (log(r+s)/(rs))}　　(r＞1 のとき)
S = 2π{1 + (arcsin(s)/(rs))}　　(r＜1 のとき)