% FY1003/TFY4155 Elmag vår 2015.
% Øving 3, opg. 3, graf av E og rho
clear all;
% x = r/R deles opp i tre ulike intervall:
x1 = 0:0.001:0.5;
  E1 = (8/5)*(4/3)*x1;
  rho1 = 1+x1-x1; % Lite knep for at rho1 skal bli vektor
x2 = 0.50:0.001:1.0;
  E2 = (8/5)*(-2*x2.^2 + (8/3)*x2 - (1/24)*x2.^(-2));
  rho2 = 2*(1-x2);
x3 = 1.0:0.001:2.0;
  E3 = x3.^(-2);
  rho3 = x3-x3;
plot(x1,E1,'b', x2,E2,'g', x3,E3,'r', x1,rho1,'k--', x2,rho2,'k--', x3,rho3,'k--', 'linewidth', 2);
title('E/(kQ/R^2) og \rho/\alpha som funksjon av x=r/R','fontsize',12)
handle=legend('E(r/R<1/2)','E(1/2<r/R<1)','E(r>R)','\rho','location','northeast');
grid on;
xlabel('r/R','fontsize',12);
ylabel('\rho/\alpha og E/(kQ/R^2)','fontsize',12);
axis([0 2 0 1.5]);