%% Rovnice vedení tepla
% Rovnice:                         \frac{\partial u}{\partial t}=\alpha\frac{\partial^2 u}{\partial x^2}+f
% Počáteční podmínka:    u(x,0)=u_0(x)
% Okrajové podmínky:    u(0,t)=g_0(t),\ u(1,t)=g_1(t)
%
% Explicitní metoda

help solveHeatExpl

%

a=1;
f=@(x,t)0;
T=0.5;
u0=@(x)sin(pi*x);
g0=@(t)0;
g1=@(t)0;
h=0.1;tau=0.05;
[x,t,U] = solveHeatExpl(a,f,T,u0,g0,g1,h,tau);
mesh(x,t,U)

%
% Zjemněte dělení na ose x

