sahav - Archaversine hav-1(x)

function sahav (x) 
function sahavf(x)

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