function h = interptable(X, Y, u)
h = interp1(X, Y, u, 'linear', 'extrap');