% 2. gyakorlat: diffúzió 2 dimenzióban (tehát egy síkon)
tic;
dt = 0.1;
dx = 1;
dy = 1;

x = 1:dx:100;
y = 1:dx:100;
nx = length(x);
ny = length(y);
nt = 1500;

Dx=1;
Dy=1;

u = zeros(nx,ny);
u(70:80, 70:80) = 1; %kezdeti állapot (négyzet)

figure

for t = 1:nt
    u2 = u;
    
    for i =1:nx      
        for j =1:ny
            
            %peremek kezelése! ...
  
             u(i, j) =... %dU/dt=D*dU^2/dx^2 %diffúzió, szétfolyik
            
        end
    end
    
    %imagesc(u);
    mesh(u);
    zlim([0 1.5]); ylim([0 ny]); xlim([0 nx]);
     
    caxis([0 1]); %lecsengés normált színekkel jobban látszódik
    %colorbar;
    title(['lépés: ', num2str(t)]);
  %  waitforbuttonpress;
    pause(0.01);
end
toc