% zadani vstupnich funkci a parametru
f=@(x,t) zeros(size(x));
T=5;
u0=@(x)sin(2*pi*x);     % pocatecni vychyleni - funkce sinus
du0=@(x)zeros(size(x)); % pocatedni rychlost
g0=@(t)zeros(size(t));  % okrajove podminky
g1=@(t)zeros(size(t));
h=0.05;                 % krok ve smeru x
tau=0.05                % casovy krok
sigma=1;                % implicitni schema
[x,t,U] = solveWaveEq(f,T,u0,du0,g0,g1,h,tau,sigma);
plot(x,U(1,:))
axis([0,1,-1,1])
pause
for k=1:length(U), plot(x,U(k,:)),axis([0,1,-1,1]), pause(0.2);end
