%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 	Solves 1d difussion equation using FTCS scheme
% 	Input data
	dx = 0.1;    % Horizontal resolution, Dominio [0,1]
	alfa = 1e-4; % Diffusivity
 	s = 0.2;	 % alfa*deltat/(dx^2)
    tmax = 1000;  % Maximum time
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%
xi=(0:dx:1);      % Dominio [0,1] con resolución dx
dt=s*dx*dx/alfa;  % dt condicional a dx, alfa y s
ti=(0:dt:tmax);
%
T=zeros(length(ti),length(xi));
T(:,1)=1;T(:,end)=1;             % CB: T=1 en x=0 y x=1

for iT=2:length(ti)
    
    T(iT,2:end-1)=T(iT-1,2:end-1)+s*(T(iT-1,1:end-2)-2*T(iT-1,2:end-1)+T(iT-1,3:end));
    
end

plot(xi,T),grid