20/12/10 21:50:59.93 Whv7aYu3.net
>>918（訂正）
URLﾘﾝｸ(i.imgur.com)
> green = function(x) sqrt(r^2-x^2) # 円Oの方程式
> blue = function(x) sqrt(s^2-(x-c)^2) # 円Cの方程式
> red = function(x) sqrt(R^2-(x-a)^2) # 円Aの方程式
>
> VG = integral(function(x) pi*green(x)^2, -r,b0) # volume green
> VB = integral(function(x) pi*blue(x)^2, b0,c+s) # volume green
> VR = integral(function(x) pi*red(x)^2, b0,a+R) - VB # volume red
> c(VG,VB,VR)
[1] 0.86504771150 0.04346965773 2.39616784810
>
> AG = integral(function(x) 2*pi*green(x), -r,b0) # area green
> AB = integral(function(x) 2*pi*blue(x),b0,c+s) # area blue
> AR = integral(function(x) 2*pi*red(x)^2, c+s,a+R) # area red
> c(AG,AB,AR)
[1] 2.957868072 0.183726052 4.689248957