% Forberedelser
%-------------
h=0.01;      % Skrittlengder i x- og t- retning.
k=0.01;
xend = 1;
tend = 6;
n = xend/h+1;   % Antall diskretiseringspunkter i hhv. x- og t-
m = tend/k+1;   % retning, randpunktene inkludert.
mu=1.e-3
u = zeros(n,m); % Lager en matrise som etterhvert skal inneholde
                % l{\o}sningen
x=[0:h:xend]';
t=[0:k:tend];
u(:,1) = 1.5*x.*(1-x).^2; % Legger startverdiene i 1. kolonne i u.

% Hovedloekke
%-------------
for j=1:m-1
  u(:,j+1) = u(:,j); % Setter startverdien til Newton-iterasjonene
                     % lik verdien i forrige skritt.
  du = ones(n-2,1);  
  
  % Newton-iterasjoner.
  while norm(du,2) > 1.e-4
%====================================================================
%    Denne biten maa du fylle ut selv. Den best{\aa}r av:
%     Regn ut funksjonen f(u):
%     Finn Jacobi-matrisa :
%     Sett: du = -J\f (l{\o}s ligningssystemet) og sett u(2:n-1,j+1) =
%          u(2:n-1,j+1) + du;
%====================================================================
  end;
end;

% Plotter l{\o}sningen
mesh(x,t,u');
view( -15,35);  % Bestemmer vinkelen du ser plottet fra.
xlabel('x');
ylabel('t');
