#include <math.h>
#define M_PI_2 3.14159/2.0
/*
XOPEN(4)
*/