a=0;
b=pi/2;
eps = 1e-4;
%   xvect - ciag kolejnych przyblizen  x(i)
%    xdif - ciag roznic |x(i+1) - x(i)|
%      fx - ciag wartosci funkcji f(x(i))
%  it_cnt - liczba wykonanych iteracji
[xvectb, xdifb, fxb, it_cntb] = bisect(a, b, 1000, eps, 'funckcjaf3');
qb = [1:size(xvectb,1)];

[xvects, xdifs, fxs, it_cnts] = secant(a, b, 1000, eps, 'funckcjaf3');
qs =[1:size(xvects,1)];

[xvectn, xdifn, fxn, it_cntn] = newton(b/2, 1000, eps, 'funckcjaf3', 'funkcjaf3dif');
qn = [1:size(xvectn,1)];

plot(qb, xvectb, 'rx-');
hold on
plot(qs, xvects, 'gx-');
hold on
plot(qn, xvectn, 'bx-');

qb = [1:size(xdifb,1)];

qs =[1:size(xdifs,1)];

qn = [1:size(xdifn,1)];
figure
plot(qb, xdifb, 'rx-');
hold on
plot(qs, xdifs, 'gx-');
hold on
plot(qn, xdifn, 'bx-');