function [f]=priklad3(l,k,n,pred)
%[f]=priklad3(l,k,n,pred)
%
% Funkce najde maximalni hustotu sily pusobici na strunu tak, aby se pri urcenych datech
% neprohla vice, nez je predepsano.
%
% l    ... delka struny
% k    ... tuhost struny
% n    ... pocet deleni
% pred    ... maximalni predepsane prohnuti struny (pred je zaporne)
%
% f ... maximalni hustoty sily, ktera muze pusobit na strunu, aby se
% prohnula nejvyse o hodnotu pred (f je zaporne)

h=l/n;
x=0:h:l;

f=0;

U=f/(2*k).*x.*(l-x);

while min(U)>=pred
    f=f-1e-2*k;
    U=f/(2*k).*x.*(l-x);
end
f=f+1e-2*k;
U=f/(2*k).*x.*(l-x);

plot(x,U)
fprintf('Hustota sily pusobici na strunu muze byt nejvyse %f.\n',f)

%konec