>>109
モンテカルロ法での判定関数が思いつかないのでガウスグリーンの定理
  https://manabitimes.jp/math/900
に数値微分と数値積分を使って算出。

> fx1=\(p) genD(fx,p)$D[1]
> fy1=\(p) genD(fy,p)$D[1]
> func=\(p) fx(p)*fy1(p)/2-fy(p)*fx1(p)/2
> func=Vectorize(func)
> integrate(func,0,Inf,rel.tol = 1e-12)
0.25 with absolute error < 5.5e-13

答は0.25であるみたいだな。
なんとかモンテカルロ法で検証したいのだが、どうしたらよいものやら。