scahav - Complex Archaversine hav-1(x)

#include <msca.h> 
 
double complex scahav (double complex x); 
float complex  scahavf(float complex x);

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