function T = waermeleitung_loesung(x, t, T0, kappa)
% WAERMELEITUNG_LOESUNG - Analytische L�sung des W�rmeleitungsproblems
%   
%
maxiter = 20;

T = 0;
k = maxiter;
while (k > 0)
  delta = T0 * 4 * sin((2*k - 1) * pi * x) / (2*k - 1) / pi ...
          * exp(- kappa * (2*k - 1)^2 * pi^2 * t);
  T = T - delta;
  %disp(delta)
  %disp(T)
  k = k - 1;
end
T = T + T0;

%% EOF
