function [k]=optstruna(l,n,f,U_pred)
%[k]=optstruna(l,n,f,U_pred)
%
% Funkce najde tuhost struny k tak, aby U se rovnalo priblizne U_pred, 
% tj. resi ulohu min ||U - U_pred||^2, kde U je deformace struny
%
% l ... delka struny
% n ... pocet deleni
% f ... sila pusobici na strunu
% U ... vysledna deformace struny
%
% k ... tuhost pruziny

[A,F]=struna_data(l,n,f);

global U_pred A F;

k=newton('opt','dopt','ddopt',1,1e-5);

clear global

%end

