clear all
format long e

n = 5;
x_vec = zeros(n,1);
f_vec = zeros(n,1);

x_vec(1,1) = 1.0;
x_vec(2,1) = 1.3;
x_vec(3,1) = 1.6;
x_vec(4,1) = 1.9;
x_vec(5,1) = 2.2;

f_vec(1,1) = 0.7651977;
f_vec(2,1) = 0.6200860;
f_vec(3,1) = 0.4554022;
f_vec(4,1) = 0.2818186;
f_vec(5,1) = 0.1103623;

x = 1.5;
[fun_val] = neville(x, x_vec, f_vec);
