function [U]=struna2(l,n,f,k)
%[U]=struna2(l,n,f,k)
%
% Funkce vypocte a vykresli prohnuti struny zatizene silou s konstantni hustotou f,
% tj. vyresi rovnici -k*u''(x) = f s okrajovou podminkou u(0)=u(l)=0.
%
% l ... delka struny
% n ... pocet deleni
% f ... konstantni hustota sily pusobici na strunu
% k ... tuhost struny
%
% U ... vysledna deformace struny

% Bez Dirichletovych podminek

% PREPROCESSING

% Diskretizace
h=l/n;

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

% SOLVER

% Reseni soustavy
U(1) = 0;
U(2:n) = (k*A)\F;
U(n+1) = 0;

% POSTPROCESSING

% Vykresleni reseni
x=0:h:l;
plot(x,U)

%end

