shavd - Haversine hav(x)

#include <msca.h> 
 
double shavd (double x); 
float  shavdf(float x);

Given a value $x$ this function computes the haversine of $x$ given in degres, $hav(x)$ .