membuat program(fisika komputasi)

kesalahan pemotongan (round off)

% rndoff - demo program untuk round-off error
clear; help rndoff;
h=1 % intial value for h
for j=1:21
    temp = (10 + h )-10; % bila tidak ada round off temp=h
    hplot(j)= h;
    eplot (j)= abs (temp-h)/h;
    h=h/10;
end
loglog (hplot, eplot,'*');
xlabel('h');
ylabel('fractional error');
title('round off error')

hasilnya