Geheugensteuntje - Wiki

Daarna stijgt de totale fout weer omdat de afrondingsfouten nu sterker stijgen.
x = 1/16 is klein; t.e.m. x = 1/512 geen probleem:
[5*ex - 5]/x = 5*[ex - 1]*(1/x)
met ex voldoende verschillend van 1.
x = 1/1024 (= 1/(210)) te klein:
gevaarlijke aftrekking
=> groot verlies aan j.b.c.
=> grote fout.
Oefening 4:
Mogelijke commando's:
K = 5
x_waarden = 0:0.01:1
functie_waarden = atan(K*x_waarden)
figure
plot(x_waarden, functie_waarden, 'b')
x_1 = 1
aantal = 50
perturbatie = 1*10^(-15)
[abskond, relkond, f_x, f_pert_x, x, x_pert] = oefz2oef4(x_1, aantal, K, perturbatie)
figure
semilogy(1:aantal,abskond,'b',1:aantal,relkond,'r')
figure
semilogy(x,abskond,'b',x,relkond,'r')
abs_verschil = abs(f_x - f_pert_x)
rel_verschil = abs_verschil ./ abs(f_x)
figure
semilogy(x, abs_verschil, 'b', x, rel_verschil, 'r')
figure
plot(x, f_x, 'b');
figure
plot(x_pert, f_pert_x, 'r');
neg_x = zeros(size(x)) - x;
figure
plot(neg_x, f_x, 'g');
neg_x_pert = zeros(size(x_pert)) - x_pert;
figure
plot(neg_x_pert, f_pert_x, 'm')
figure
plot(x, f_x, 'b');
hold
plot(neg_x, f_x, 'g');
figure
plot(x_pert, f_pert_x, 'r');
hold
plot(neg_x_pert, f_pert_x, 'm');
figure
semilogy(x, abskond, 'b');