function res=myjinc(x)
%res=jinc(x)
x=x+eps; % avoids errors at x=0

[tmp,err] = besselj(1,pi*x);
res = real(tmp./(2*x));
if sum(err(:))
    disp('error');
end


