% skrypt 32
disp('skrypt oblicza pierwiastki, maksima i minima funkcji bessla pierwszego rodzaju')
z=zeros(1,5);
%% liczymy mielsca zerowe
z(1)=fzero('be',0);
z(2)=z(1);
for ji=2:5,
    while z(ji)==z(ji-1) 
        x=3;
        z(ji)=fzero('be',z(ji-1)+x);
        if ji<5,
            z(ji+1)=z(ji);
        end
        x=x+0.5;
    end
end
disp(sprintf('miejsce zerowe to: %d \n', z));
%% liczymy minima
for jk=(1:2),
    mi(jk)=fminbnd('be',z(jk+1),z(jk+3));
end

disp(sprintf('minimum to: %d \n', mi));
li
for jk=(1:2),
    mx(jk)=fminbnd('bem',z(jk),z(jk+2));
end
mx(3)=fminbnd('bem',z(4),15);

disp(sprintf('maximum to: %d \n', mx));

cla;
hold on;
fplot('be',[0 15],'-k');
plot(z,[0 0 0 0 0],'xr');
plot(mi,be(mi),'+b');
plot(mx,be(mx),'+g');
hold off;

