Mine

%% Task 1 = Animatie in de tijd van Gaussian.
L=1;
kappa=1;
x=-L/2:L/100:L/2;
vf=20; %vertragingsfactor
dt = 0.01/vf;
t=0;
%y=(4*pi*kappa*t)^-(1/2)*exp(-x.^2/(4*kappa*t));
y=100*dirac(x);
h = plot(x,y);
axis([-L/2 L/2 0 10]);
title('Task1, animatie van de warmte diffusie, Analytical solution'); xlabel('Plaats'); ylabel('Graden 

Celsius');
set(h,'EraseMode','xor','MarkerSize',18);

while t<=0.01
   drawnow
   t=t+dt;
   x=x;
   y=(4*pi*kappa*t)^-(1/2)*exp(-x.^2/(4*kappa*t));
   set(h,'XData',x,'YData',y)
   pause(dt);
end;

y_ana=y;