%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Desenvolvido por Leonardo M. N. de Mattos. %
% www.mattos.eng.br                          %
% l (at) mattos.eng.br                       %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function dominio_novo = Tprox_ret(dominio,Fo,deltax,deltay,L,H)
% Tprox_ret Calcula a proxima temperatura de um dominio retangular
% Com as condições de contorno dadas anteriormente

dominio_novo = dominio;

for n = 2:H/deltay-1
   for m = 2:L/deltax-1       
       dominio_novo(n,m) = Fo*(dominio(n+1,m) ...
       +dominio(n-1,m) ...
       +dominio(n,m+1) ...
       +dominio(n,m-1)) ...
       +(1-4*Fo)*dominio(n,m);
   end 
end

end %fim do da função
