もうしらん
7c7
< }yumiya[1000];
---
> }yumiya[N];
15,18c15,16
< {
< if((z->x)*(z->x)+(z->y)*(z->y)<=1){
< return 0;
< }else{
---
> {
> if((z->x)*(z->x)+(z->y)*(z->y)<=1.){ // sqrt(x*x + y*y) <= 1 <=> x*x + y*y <= 1
19a18,19
> }else{
> return 0;
30c30
< return double(a/(a+b));
---
> return a/double(a+b);
40c40,41
< n+=det_arrow(yumiya[i]);
---
> set_arrow(&yumiya[i]);
> n+=det_arrow(&yumiya[i]);
45c46
< printf("%lf",probabi(n,h));
---
> printf("%f\n",4*probabi(n,h)); // quadrant sphere