function gleba
close all
m = 0;
z = [0 0.005 0.01 0.05 0.1 0.2 0.5 0.7 0.9 0.95 0.99 0.995 1];
t = [0:0.01:24];
T0=10;
sol = pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,z,t);

figure
plot(t,sol)
%ylabel('Deph [m]')
%xlabel('Time [h]')

figure
%pcolor(sol)
imagesc(t,z,sol')
colorbar
%axis xy
shading flat
% --------------------------------------------------------------
function [c,f,s] = pdex4pde(x,t,u,DuDx)
k=1e-2;
c = 1; 
f = k*DuDx; 
s = 0; 
% --------------------------------------------------------------
function u0 = pdex4ic(x);
T0=10;
u0 = T0; 
% --------------------------------------------------------------
function [pl,ql,pr,qr] = pdex4bc(xl,ul,xr,ur,t)
T0=10;
pr = ur-T0; 
qr = 0; 
pl = ul-T0*(1+0.5*sin(pi*(t-6)/12));  ;
ql = 0;   

%pl = T0*sin(pi*t/12); 
%ql = -T0*sin(pi*t/12); 