%--------------------------------------------------------------
% Konturlinje- og gradientvektorplott av fjelltoppen Berenberg
%--------------------------------------------------------------

h0 = 2277;                          % Høyden av toppen av fjellet (m)
R = 4;                              % Mål for radius av fjellet (km)

[x,y] = meshgrid(-10:.1:10);        % Grid for x- og y-verdiene (km)
h = h0./(1+(x.^2+y.^2)./(R.^2));    % Beregn høyden h (m)
contour(x,y,h,'b')                  % Plott konturlinjene (blå farge)

hold on                             % Behold konturlinjene og akse-egenskapene

[x,y] = meshgrid(-10:2:10);         % Definer et grovere grid til vektorfeltet
h = h0./(1+(x.^2+y.^2)./(R.^2));    % Beregn høyden med det nye griddet

[dbx,dby] = gradient(h);            % Beregn gradientvektoren (dh/dx, dh/dy)
quiver(x,y,dbx,dby,0.75, 'r')          % Plott vektorfeltet (rød farge) og skaler
                                        % vektorlengden med en faktor 0.75
hold off                            % Trenger ikke flere plott i denne figuren
axis square                         % Sett aksestørrelsene like