w=@(x)((cos(4.*x)).*x.^(3./2) -4*(x-10));
px=linspace(5,20);
plot(px,w(px));
hold on
X=fzero(@jestes,8);

for n=5:20
    X(n-4)=fzero(@jestes,n);
end
Xkoniec=unique(X);
plot(Xkoniec,0,'ro');


A jestes to taka funkcja
function G = jestes(x)
    G=[((cos(4.*x)).*x.^(3./2) -4*(x-10))];
end