%% Eksamen 2009, oppgave 3 i bolgefysikk, program for aa regne ut
%% og plotte interferensmonstret med multiple spalter, med spaltebredde a og
%% spalteavstand d

%% Bolgelengde:
lambda = 600e-9;
%% Lager og plotter I(theta) for N=6 spalter
%% og med a=0.01 mm og d=0.100 mm:
%% Litt proving og feiling tilsier at maks vinkel heller ikke her bor vare mer enn pi/90
theta = 1e-10:pi/1e5:0.020;
%% Vi maa ogsaa naa bruke hele funksjonen oppgitt i oppgaveteksten, 
%% men beta og phi trenger naa ikke oppdateres. Funksjonen er
%% [sin(beta)/beta]^2 * [sin(N*phi)/sin(phi)]^2, med
%% beta = pi*a*sin(theta)/lambda og
%% phi = pi*d*sin(theta)/lambda
a=0.009e-3;
d=0.100e-3;
beta=pi*a*sin(theta)/lambda;
phi=pi*d*sin(theta)/lambda;
I6 = ((sin(beta)./beta).^2).*((sin(6*phi)./sin(phi)).^2)/36;
%% Diffraksjonsfaktoren:
D =  (sin(beta)./beta).^2;
%% Deretter plotter vi intensitetsfordelingen som funksjon av 
%% vinkelen theta:
%% plot(theta,I6,theta,D);
plot(theta,I6);
grid on;
xlabel('theta (radianer)');
ylabel('Normert intensitet');
%% Skriver denne figuren til en eps-fil med navn oppgave3.eps
print('oppgave3.eps');