21/02/14 20:58:48.87 68y2TfCD.net
>>83
こっちのほうが楽だった。pracmaはベクトルの外積計算でしか使ってなかった。
密度z = 1/((1+sqrt(x^2+y^2)) * 1/(1+sqrt((x-1)^2+y^2)))を
x=r*cos(t)
y=r*sin(t)
としてヤコビアンをかけて
A＝∫[0,2π]∫[0,2] r/(1 + r) * 1/(1 + sqrt((r*cos(t)-1)^2 + (r*sin(t))^2)) dr dt
B＝∫[0,2π]∫[1,2] r/(1 + r) * 1/(1 + sqrt((r*cos(t)-1)^2 + (r*sin(t))^2)) dr dt
を求める
"
library(pracma)
f <- function(r,t) r/(1+r) * 1/(1 + sqrt((r*cos(t)-1)^2 + (r*sin(t))^2))
(A=integral2(f, 0,2,-pi,pi)$Q)
(B=integral2(f, 1,2,-pi,pi)$Q)
B/A*4000