double angle = ( double )( pos[e] - pos[b] ) * pi * 2.0 / ( double )n;

	return ( double ) d * ( double ) d / 4.0 * sin ( angle ) / 2.0;