% 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

% Sätt antal intervall och steglängd
n = ...
h = ...
% skapa vektor med nodpunkter
xv=0:h:pi;  

% Sätt diffusionskonstant
D = ...

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

A = D*S + 0.5 * C;

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

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


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