function dy = besselde(x,y)
%  dy = besselde(x,y)
%  Bessel equation of order 1/2 for 2-point BVP

dy = [ y(2);
       -(1-1/(4*x^2))*y(1)-y(2)/x ];