besselJ(a_,x_):= if(a<0, -besselJ(-a,x) );
besselJ(1/2,x_):=sqrt(2/pi/x)*sin(x);