sahavd - Archaversine hav-1(x)

#include <msca.h> 
 
double sahavd (double x); 
float  sahavdf(float x);

Given a value $x$ this function computes the archaversine of $x$ with output in degrees, $hav-1(x)$ .