function [A,F]=struna_data(l,n,f)
%[A,F]=struna_data(l,n,f,k)
%
% Funkce sestavi matici A a vektor F struny zatizene silou s konstantni hustotou f
%
% l ... delka struny
% n ... pocet deleni
% f ... konstantni hustota sily pusobici na strunu

% Diskretizace
h=l/n;

% Sestaveni soustavy po diskretizaci
A=(1/h^2).*(diag(2*ones(n+1,1)) + diag(-ones(n,1),1) + diag(-ones(n,1),-1));
F=f*ones(n+1,1);

% Dirichletova podminka
A([1,n+1],:)=zeros(2,n+1);
A(:,[1,n+1])=zeros(n+1,2);
A([1,n+1],[1,n+1])=eye(2,2);

F(1)=0;
F(n+1)=0;

%end

