function [xmin,fmin]=MinimumBrutal(jmenof,a,b,nmax)
% Nalezeni minima funkce f(x) metodou hrube sily
% [xmin,fmin]=MinimumBrutal(jmenof,a,b,nmax);
% xmin ... souradnice minima
% fmin ... nalezena minimalni hodnota f(x)
% jmenof ... nazev funkce f(x) jako retezec znaku
% a ... dolni mez definicniho oboru f(x)
% b ... horni mez definicniho oboru f(x)
% nmax ... maximalni pocet vycisleni funkce f(x)
h=(b-a)/(nmax-1); 	% krok vycisleni funkce
fmin=inf; 	% doposud nejlepsi reseni
for k=1:nmax 	% bezpecny cykl
   x=a+h*(k-1); 	% novy bod
   f=feval(jmenof,x); 	% vycislit v novem bode
   if f<fmin 	% je to lepsi?
      fmin=f; 	% aktualizace hodnoty minima
      xmin=x; 	% aktualizace polohy minima
   end
end