function [x0,iteracje]=fals( f,a,b)

p=0.000001; z=0.000001;

ilosc=0; f=inline(f);

fa=f(a); fb=f(b);

x1=a; x0=b;

if(fa*fb > 0)

x0='Funkcja nie spełnia założeń'

iteracje=0 else

while(abs(x1-x0) > p) x1=x0;

x0=a-(fa*((b-a)/(fb-fa)));

f0=f(x0);

if(abs(f0) < z) break; end

if(fa*f0 < 0) b=x0; fb=f0;

else

a=x0; fa=f0; end

ilosc=ilosc+1; end

iteracje=ilosc

x0=x0

end end

function [x0,iter]=sampol( f,a,b)

prec=0.000000001; zer=0.0000000001;

count=0; f=inline(f);

fa=f(a); fb=f(b);

if(fa*fb > 0)

x0='Funkcja nie spełnia założeń'

iter=0 else

while(abs(a-b) > prec)

x0=(a+b)/2;

f0=f(x0);

if(abs(f0) < zer) break; end

if(fa*f0 < 0)

b=x0; else a=x0;

fa=f0; end

count=count+1; end

iter=count

x0=x0 end end

function[w,k]=samsiecz( f,a,b )

EPSX=0.000000001; EPS0=0.000000001;

count=0; f=inline(f);

f1=f(a); f2=f(b); i=500;

while(i>0 && (abs(a-b) > EPSX))

if(abs(f1-f2) < EPS0)

w='Złe punkty startowe'

k=0 break; end

x0=a-f1*((b-a)/(f2-f1));

f0=f(x0);

if(abs(f0) < EPS0) break; end

b=a; f2=f1; a=x0; f1=f0;

count=count+1; i=i-1; end

if(i==0) x0='Za dużo obiegów'; end

w=count k=x0 end