function vals=testorder(nd)
h(1)=0.5;
for i=2:20
    h(i) = h(i-1)/2;
end

for i=1:20
    er(i)=abs(cos(1)-nd(1,h(i), @sin));
end

for i=1:19
    vals(i)=log(er(i))/log(h(i));
end

