function [y,dy] = f2(x)
    y  = +besselj(0,x);
    dy = -besselj(1,x);
    [x,y]  % displayed for each evaluation of the function
end%function 