%
%   Ex3e.m
%
%  Dette skriptet gj?r en del av plottet i oppgave 3e i Matlab?vingen
%  
%
hold off
n=20;
% Vi begynner med ytre kuleflate (konstant radius rho=2)
theta = (0:n)/n*pi/2;     % Asimuth-vinkel theta [0,pi/2]
phi = (0:n)'/n*pi/2;      % Elevasjonsvinkel phi [0,pi/2]
cosphi = cos(phi);
sinphi = sin(phi);
sintheta = sin(theta);
costheta = cos(theta);
x = 2*cosphi*costheta;     
y = 2*cosphi*sintheta;
z = 2*sinphi*ones(1,n+1);
h1=surf(x,y,z);

colormap(gray)
shading flat


hold on;

% Plotting av sidevegg y=0
rr = 2*(0:1:n)/n;
x = cosphi*rr;
z = sinphi*rr;
y = zeros(size(x));
h2=surf(x,y,z);

% Plotting av sidevegg x=0
y = cosphi*rr;
z = sinphi*rr;
x = zeros(size(y));
h3=surf(x,y,z);
 
% F?lgende kommando gj?r at flatene vi plottet blir gjennomsiktige
% Jo mindre alpha, jo mer gjennomsiktig.
alpha(0.7);    
axis equal
xlabel('x'); ylabel('z'); zlabel('z');


% For ? fullf?re oppgaven m? man plotte en passende del av den indre
% kuleflaten. Greier du dette? Du trenger 4-5 kommandoer


