
% Mall för Studio II : värmeledningsekvationen
% Löser värmeledningsekvationen på [0,1] med alpha=1,
% Neumann-randvillkor för x=0 och homogent Dirichlet-villkor för x=1

% sätt antal delintervall, steglängd och nodpunkter
n = ...
h = ...
xv = 0:h:1;

% sätt tidsintervall, steglängd i tiden och antal tidssteg
T = ...  % tidsintervall [0,T]
m = ...  % antal tidsintervall
k = T/m; % tidssteg

% skapa matris med lösning för alla rums- och tidssteg (inklusive randvillkor)
U=zeros(n+1,m+1);

% sätt begynnelsevillkor (och randvillkor) vid t=0
U(:,1) = ...
U(n+1,1) = ...

% Skapa matriser
M = ...
S = ...

% Tidsstegning
% OBS att indexeringen i Matlab är från 1, därför l+1 nedan
for l=1:m

   % Skapa högerledet 
   F = ...

   % Lös för nästa tidssteg enl ekv (11)
   xi = ...
   U(1:n,l+1) = xi;

   % Sätt eventuellt randvillkor
   U(n+1,l+1) = ...
end


% Plotta resultat, t.ex. med plot för varje tidssteg, eller surf
