function [U]=animace(l,n,g,f,k,v)
%[U]=animace(l,n,g,f,k,v)
%
% Funkce vypocte a vykresli animaci prohnuti mostu zatizeneho gravitacni silou s konstantni hustotou g a
% bodovou silou f v bode, ktery se pohybuje od bodu 0 do bodu l konstantni rychlosti v
% tj. vyresi rovnici -k*u''(x) = g + f*delta_vt(x) s okrajovou podminkou u(0)=u(l)=0.
%
% l ... delka mostu
% n ... pocet deleni
% g ... konstantni hustota gravitacni sily pusobici na most
% f ... bodova sila pusobici na most v bode x
% k ... tuhost mostu
% v ... konstantni ryhlost, kterou se pohybuje hmotny bod
%
% U ... vysledna deformace mostu

t=l/v;  %cas, za ktery bod "prejde" most
dt=t/n;

for i=0:n
    pause(dt);
    [U(:,i+1)]=most(l,n,g,f,v*dt*i,k);
    axis([0 l -1 0]); % muzete zkusit i jinou volbu napr. axis([0 l -0.4 0]);
end

%end