
function demo_interpolazione_inversa

x=2.1:0.1:2.9;
y=besselj(0,x);
xbar=interpolazione_inversa(x,y);

fprintf('\n \t xbar: %1.15e',xbar);
sol=2.404825557695780;
err=abs(sol-xbar);
fprintf('\n \t err: %1.3e \n \n',err);

%--------------------------------------------------------------------------
% >> demo_interpolazione_inversa
% 
%  	 xbar: 2.404825571062216e+00
%  	 err: 1.337e-08 
%  
% >> 
%--------------------------------------------------------------------------