sahavpi - Half-Cycle Archaversine hav-1(x)/π

#include <msca.h> 
 
double sahavpi (double x); 
float  sahavpif(float x);

Given a value $x$ this function computes the half-cycle archaversine of $x$ , $hav-1(x)/π$ .