close all;
clear;
clc;

%---------------------------------------
%   A
%---------------------------------------

xa=linspace(15,55);
fxa=sin(xa).*xa.^(3/2) - 5.*(xa-10);

plot(xa,fxa,'-k','LineWidth',2), grid on, hold on
xlabel('X'), ylabel('F(x)')

%---------------------------------------
%  B
%---------------------------------------

xb=linspace(20,50,21);
k=length(xb);
F=inline('sin(xb).*xb.^(3/2) - 5.*(xb-10)');

%szukamy miejsc zerowych
for i=1:k
    x0(i)=fzero(F,xb(i));  % x0 - wektor miejsc zerowych
end

%----------------------------------------
%  C
%----------------------------------------

x0=sort(x0);

%usuwanie miejsc zerowych nie mieszczących się w przedziale 20-50

n=1;
indeksy = find(20<x0 & x0<50); % zwraca indeksy wektora x0 które spełniają
                               % warunek 20<x0 & x0<50
for k=indeksy
    tab(n)=x0(k);              % prezniesienie wartości spelniających warunek zad.
    n=n+1;
end

tab

% usuwanie powtarzających się elementów

tab2=unique(tab)


%---------------------------
%  D
%---------------------------

plot(tab2,0,'-ro')