% Mall för Studio 1, FEM för stationärt konvektions-diffusions-problem
% Löser -Du'' + 1/2*u' = 1 på [0,pi] med homogena Dirichlet-randvillkor
% Minor changes from Fardin Saedpanah's code

% Sätt antal obekanta nodvärdena och steglängd
m = ... ;
h = ... ;
% skapa vektor med nodpunkter
xv=0:h:pi;  

% Sätt diffusionskonstant
D = ... ;

% Skapa matriser
% Görs här enkelt med en loop. 
% Ni kan också använda Matlabs kommandon diag eller spdiags
S = zeros(m,m);
C = zeros(m,m);
b = zeros(m,1);
for j=1 : ...
  % Fyll i själva!
  S(j,...) = ... ;
  C(j,...) = ... ;
  b(j) = ... ;
end

A = ... ;

% Lös ekvationssystemet A xi = b
xi = ... ;

% lägg till randvärden till xi
xi = ... ;


% rita lösningen
figure(1)
plot(xv, xi, 'r-')
