#include <math.h>
#define M_SQRT1_2 1.0/sqrt(2.0)
/*
XOPEN(4)
*/