sagd - Inverse Gudermannian function gd-1(x)

#include <msca.h> 
 
double sagd (double x); 
float  sagdf(float x);

Given a value $x$ this function computes the inverse Gudermannian function, $gd-1(x)$ .