shavpi - Half-Cycle Haversine hav(πx)

#include <msca.h> 
 
double shavpi (double x); 
float  shavpif(float x);

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