sj1 - Bessel function J1(x)

#include <msca.h> 
 
double sj1 (double x); 
float  sj1f(float x);

Given a value $x$ this function computes the Bessel function $J (x) 1$ .