% 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
% Mindre ändringar från Fardin Saedpanah's koda

% sätt antal obekanta nodvärdena, steglängd och nodpunkter
m = ... ;
h = ... ;
xv = 0:h:1;

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

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

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

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

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

   % Skapa högerledet 
   F = ... ;

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

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


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