%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                Losningsforslag til Plotting opg.4 Øv 4 TFY4155             %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear; % Nullstille alle variabler og deres dimensjoner
x = [0 : 1/20 : 1.5];  % Verdier for x (=r/R)
sizex = size(x); %dimensjon til x, (#row, #col)
for ii = 1:sizex(2) 
  if x(ii)<1.0
    rho(ii) = 1 - x(ii);
    Q(ii) = 4*x(ii)^3 - 3*x(ii)^4;
    E(ii) = 4*x(ii) - 3*x(ii)^2;
    V(ii) = 2 - 2*x(ii)^2 + x(ii)^3;
  else
    rho(ii) = 0;  
    Q(ii) = 1;
    E(ii) = x(ii)^(-2);
    V(ii) = x(ii)^(-1);
  end;
end

plot(x,rho,'-', x,Q, '--', x,E,':', x,V,'-.', 'LineWidth',2.5);
title('rho, Q, E og V for inhomogent ladd kule')
xlabel('r/R'); % Tekst på aksene
ylabel('Rel. rho, Q, E og V');
handle=legend('rho', 'Q', 'E', 'V','location','northeast');
set(handle, 'Box', 'off'); % Ikke boks rundt legend
