% Skript zur Musterlösung
%
% Numerische Mathematik, SS 2013
%
% Visualisierung Rundungsfehler
%
% Ingolf Busch, Juni 2013


k = 0:360;
x = 1.606+k*eps;
a = [4 -59 324 -751 622];
b = [1 -14 72 -151 112];

figure;
y = 10^12*(my_horner(x,a,b)-8.752376580778);
plot(k,y,'+');