sj1 - Bessel function J1(x)

function sj1 (x) 
function sj1f(x)

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