
function demo_regula_falsi

% nota: le operazioni sono scalari, non serve il "." ma lo metto perche'
% impariate a usarlo.
f=@(x) x.^3-4*x.^2+5*x-2;
sol=2;
a=1.5; b=2.5; toll=1e-10; maxit=10000;

% applicazione regula falsi
[aa,bb,c,wr,flag]= regula_falsi(f,a,b,toll,maxit);

% estremi intervalli
fprintf('\n \n ------ ESTREMI INTERVALLI ------ \n \n');
for k=1:length(aa)
    fprintf('\n \t k: %5.0f aa: %1.15e bb: %1.15e',k,aa(k),bb(k));
end

% lunghezze intervalli
ampiezze=abs(aa-bb);
fprintf('\n \n ------ AMPIEZZE INTERVALLI ------ \n \n');
for k=1:length(ampiezze)
    fprintf('\n \t k: %5.0f ampiezze: %1.15e',k,ampiezze(k));
end

% residui pesati
fprintf('\n \n ------ RESIDUI PESATI ------ \n \n');
for k=1:length(wr)
    fprintf('\n \t k: %5.0f wr: %1.5e',k,wr(k));
end

% errori assoluti
abserr=abs(c-sol);
fprintf('\n \n ------ ERRORI ASSOLUTI ------ \n \n');
for k=1:length(abserr)
    fprintf('\n \t k: %5.0f abserr: %1.5e',k,abserr(k));
end

fprintf('\n \n');

% grafico errori assoluti
semilogy(1:length(abserr),abserr);








% % COMMENTO
% % 1. IL PLOT MOSTRA UNA CONVERGENZA LINEARE.
% % 2. LE AMPIEZZE SONO TALI PERCHE' c(k) SOSTITUISCE SEMPRE a(k).

% % demo_regula_falsi
% %
% %  	 k:     1 aa: 1.500000000000000e+00 bb: 2.500000000000000e+00
% %  	 k:     2 aa: 1.600000000000000e+00 bb: 2.500000000000000e+00
% %  	 k:     3 aa: 1.702127659574469e+00 bb: 2.500000000000000e+00
% %  	 k:     4 aa: 1.794249201277955e+00 bb: 2.500000000000000e+00
% %  	 k:     5 aa: 1.867251073444433e+00 bb: 2.500000000000000e+00
% %  	 k:     6 aa: 1.918829899212133e+00 bb: 2.500000000000000e+00
% %  	 k:     7 aa: 1.952198429137602e+00 bb: 2.500000000000000e+00
% %  	 k:     8 aa: 1.972519697635290e+00 bb: 2.500000000000000e+00
% %  	 k:     9 aa: 1.984430821580826e+00 bb: 2.500000000000000e+00
% %  	 k:    10 aa: 1.991253966794777e+00 bb: 2.500000000000000e+00
% %  	 k:    11 aa: 1.995110744470290e+00 bb: 2.500000000000000e+00
% %  	 k:    12 aa: 1.997274280581629e+00 bb: 2.500000000000000e+00
% %  	 k:    13 aa: 1.998482772413684e+00 bb: 2.500000000000000e+00
% %  	 k:    14 aa: 1.999156185690423e+00 bb: 2.500000000000000e+00
% %  	 k:    15 aa: 1.999530932864180e+00 bb: 2.500000000000000e+00
% %  	 k:    16 aa: 1.999739320203891e+00 bb: 2.500000000000000e+00
% %  	 k:    17 aa: 1.999855151041641e+00 bb: 2.500000000000000e+00
% %  	 k:    18 aa: 1.999919520067006e+00 bb: 2.500000000000000e+00
% %  	 k:    19 aa: 1.999955286367192e+00 bb: 2.500000000000000e+00
% %  	 k:    20 aa: 1.999975158303017e+00 bb: 2.500000000000000e+00
% %  	 k:    21 aa: 1.999986198813432e+00 bb: 2.500000000000000e+00
% %  	 k:    22 aa: 1.999992332598880e+00 bb: 2.500000000000000e+00
% %  	 k:    23 aa: 1.999995740309486e+00 bb: 2.500000000000000e+00
% %  	 k:    24 aa: 1.999997633498102e+00 bb: 2.500000000000000e+00
% %  	 k:    25 aa: 1.999998685274511e+00 bb: 2.500000000000000e+00
% %  	 k:    26 aa: 1.999999269596267e+00 bb: 2.500000000000000e+00
% %  	 k:    27 aa: 1.999999594219938e+00 bb: 2.500000000000000e+00
% %  	 k:    28 aa: 1.999999774566567e+00 bb: 2.500000000000000e+00
% %  	 k:    29 aa: 1.999999874759183e+00 bb: 2.500000000000000e+00
% %  	 k:    30 aa: 1.999999930421762e+00 bb: 2.500000000000000e+00
% %  	 k:    31 aa: 1.999999961345421e+00 bb: 2.500000000000000e+00
% %  	 k:    32 aa: 1.999999978525233e+00 bb: 2.500000000000000e+00
% %  	 k:    33 aa: 1.999999988069574e+00 bb: 2.500000000000000e+00
% %  	 k:    34 aa: 1.999999993371985e+00 bb: 2.500000000000000e+00
% %  	 k:    35 aa: 1.999999996317769e+00 bb: 2.500000000000000e+00
% %  	 k:    36 aa: 1.999999997954317e+00 bb: 2.500000000000000e+00
% %  	 k:    37 aa: 1.999999998863509e+00 bb: 2.500000000000000e+00
% %  	 k:    38 aa: 1.999999999368616e+00 bb: 2.500000000000000e+00
% %  	 k:    39 aa: 1.999999999649231e+00 bb: 2.500000000000000e+00
% %  	 k:    40 aa: 1.999999999805129e+00 bb: 1.999999999805129e+00
% %
% %  ------ AMPIEZZE INTERVALLI ------
% %
% %
% %  	 k:     1 ampiezze: 1.000000000000000e+00
% %  	 k:     2 ampiezze: 8.999999999999999e-01
% %  	 k:     3 ampiezze: 7.978723404255310e-01
% %  	 k:     4 ampiezze: 7.057507987220448e-01
% %  	 k:     5 ampiezze: 6.327489265555670e-01
% %  	 k:     6 ampiezze: 5.811701007878673e-01
% %  	 k:     7 ampiezze: 5.478015708623980e-01
% %  	 k:     8 ampiezze: 5.274803023647097e-01
% %  	 k:     9 ampiezze: 5.155691784191745e-01
% %  	 k:    10 ampiezze: 5.087460332052227e-01
% %  	 k:    11 ampiezze: 5.048892555297100e-01
% %  	 k:    12 ampiezze: 5.027257194183710e-01
% %  	 k:    13 ampiezze: 5.015172275863164e-01
% %  	 k:    14 ampiezze: 5.008438143095768e-01
% %  	 k:    15 ampiezze: 5.004690671358201e-01
% %  	 k:    16 ampiezze: 5.002606797961087e-01
% %  	 k:    17 ampiezze: 5.001448489583589e-01
% %  	 k:    18 ampiezze: 5.000804799329939e-01
% %  	 k:    19 ampiezze: 5.000447136328081e-01
% %  	 k:    20 ampiezze: 5.000248416969830e-01
% %  	 k:    21 ampiezze: 5.000138011865676e-01
% %  	 k:    22 ampiezze: 5.000076674011196e-01
% %  	 k:    23 ampiezze: 5.000042596905139e-01
% %  	 k:    24 ampiezze: 5.000023665018980e-01
% %  	 k:    25 ampiezze: 5.000013147254891e-01
% %  	 k:    26 ampiezze: 5.000007304037331e-01
% %  	 k:    27 ampiezze: 5.000004057800624e-01
% %  	 k:    28 ampiezze: 5.000002254334330e-01
% %  	 k:    29 ampiezze: 5.000001252408166e-01
% %  	 k:    30 ampiezze: 5.000000695782381e-01
% %  	 k:    31 ampiezze: 5.000000386545789e-01
% %  	 k:    32 ampiezze: 5.000000214747666e-01
% %  	 k:    33 ampiezze: 5.000000119304264e-01
% %  	 k:    34 ampiezze: 5.000000066280148e-01
% %  	 k:    35 ampiezze: 5.000000036822305e-01
% %  	 k:    36 ampiezze: 5.000000020456834e-01
% %  	 k:    37 ampiezze: 5.000000011364907e-01
% %  	 k:    38 ampiezze: 5.000000006313836e-01
% %  	 k:    39 ampiezze: 5.000000003507685e-01
% %  	 k:    40 ampiezze: 0.000000000000000e+00
% %
% %  ------ RESIDUI PESATI ------
% %
% %
% %  	 k:     1 wr: 1.15200e-01
% %  	 k:     2 wr: 1.04146e-01
% %  	 k:     3 wr: 8.14243e-02
% %  	 k:     4 wr: 5.61564e-02
% %  	 k:     5 wr: 3.54011e-02
% %  	 k:     6 wr: 2.11041e-02
% %  	 k:     7 wr: 1.21863e-02
% %  	 k:     8 wr: 6.91466e-03
% %  	 k:     9 wr: 3.88624e-03
% %  	 k:    10 wr: 2.17285e-03
% %  	 k:    11 wr: 1.21140e-03
% %  	 k:    12 wr: 6.74319e-04
% %  	 k:    13 wr: 3.75028e-04
% %  	 k:    14 wr: 2.08474e-04
% %  	 k:    15 wr: 1.15858e-04
% %  	 k:    16 wr: 6.43773e-05
% %  	 k:    17 wr: 3.57689e-05
% %  	 k:    18 wr: 1.98727e-05
% %  	 k:    19 wr: 1.10408e-05
% %  	 k:    20 wr: 6.13386e-06
% %  	 k:    21 wr: 3.40773e-06
% %  	 k:    22 wr: 1.89320e-06
% %  	 k:    23 wr: 1.05178e-06
% %  	 k:    24 wr: 5.84322e-07
% %  	 k:    25 wr: 3.24624e-07
% %  	 k:    26 wr: 1.80347e-07
% %  	 k:    27 wr: 1.00193e-07
% %  	 k:    28 wr: 5.56626e-08
% %  	 k:    29 wr: 3.09237e-08
% %  	 k:    30 wr: 1.71798e-08
% %  	 k:    31 wr: 9.54434e-09
% %  	 k:    32 wr: 5.30241e-09
% %  	 k:    33 wr: 2.94578e-09
% %  	 k:    34 wr: 1.63655e-09
% %  	 k:    35 wr: 9.09193e-10
% %  	 k:    36 wr: 5.05107e-10
% %  	 k:    37 wr: 2.80615e-10
% %  	 k:    38 wr: 1.55897e-10
% %  	 k:    39 wr: 8.66092e-11
% %
% %  ------ ERRORI ASSOLUTI ------
% %
% %
% %  	 k:     1 abserr: 4.00000e-01
% %  	 k:     2 abserr: 2.97872e-01
% %  	 k:     3 abserr: 2.05751e-01
% %  	 k:     4 abserr: 1.32749e-01
% %  	 k:     5 abserr: 8.11701e-02
% %  	 k:     6 abserr: 4.78016e-02
% %  	 k:     7 abserr: 2.74803e-02
% %  	 k:     8 abserr: 1.55692e-02
% %  	 k:     9 abserr: 8.74603e-03
% %  	 k:    10 abserr: 4.88926e-03
% %  	 k:    11 abserr: 2.72572e-03
% %  	 k:    12 abserr: 1.51723e-03
% %  	 k:    13 abserr: 8.43814e-04
% %  	 k:    14 abserr: 4.69067e-04
% %  	 k:    15 abserr: 2.60680e-04
% %  	 k:    16 abserr: 1.44849e-04
% %  	 k:    17 abserr: 8.04799e-05
% %  	 k:    18 abserr: 4.47136e-05
% %  	 k:    19 abserr: 2.48417e-05
% %  	 k:    20 abserr: 1.38012e-05
% %  	 k:    21 abserr: 7.66740e-06
% %  	 k:    22 abserr: 4.25969e-06
% %  	 k:    23 abserr: 2.36650e-06
% %  	 k:    24 abserr: 1.31473e-06
% %  	 k:    25 abserr: 7.30404e-07
% %  	 k:    26 abserr: 4.05780e-07
% %  	 k:    27 abserr: 2.25433e-07
% %  	 k:    28 abserr: 1.25241e-07
% %  	 k:    29 abserr: 6.95782e-08
% %  	 k:    30 abserr: 3.86546e-08
% %  	 k:    31 abserr: 2.14748e-08
% %  	 k:    32 abserr: 1.19304e-08
% %  	 k:    33 abserr: 6.62801e-09
% %  	 k:    34 abserr: 3.68223e-09
% %  	 k:    35 abserr: 2.04568e-09
% %  	 k:    36 abserr: 1.13649e-09
% %  	 k:    37 abserr: 6.31384e-10
% %  	 k:    38 abserr: 3.50769e-10
% %  	 k:    39 abserr: 1.94871e-10
% %


