h = 1;

hlist = [];
flist = [];

for n = 1:16
    h = h / 10;
    f = sin(h) / h;

    hlist = [hlist h];
    flist = [flist f];
end

% Print values as a table (matrix)
format long
[hlist' flist']
