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

%% Bolgelengde (brukte 600 nm i 2009):
lambda = 500e-9;
%% Lager og plotter I(theta) for N=8 spalter
%% og med a=0.001 mm og d=0.100 mm:
%% Litt proving og feiling tilsier at maks vinkel heller ikke her bor vare mer enn pi/90
theta2 = -0.05:pi/1e5:0.05;
theta = -0.0125:pi/1e5:0.0125;
%% 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.001e-3;
a2=0.020e-3;
d=0.100e-3;
beta=pi*a*sin(theta)/lambda;
beta2=pi*a2*sin(theta2)/lambda;
phi=pi*d*sin(theta)/lambda;
phi2=pi*d*sin(theta2)/lambda;
I8 = ((sin(beta)./beta).^2).*((sin(8*phi)./sin(phi)).^2)/64;
I8_2 = ((sin(beta2)./beta2).^2).*((sin(8*phi2)./sin(phi2)).^2)/64;
%% Diffraksjonsfaktoren:
D =  (sin(beta)./beta).^2;
D2 =  (sin(beta2)./beta2).^2;
%% Deretter plotter vi intensitetsfordelingen som funksjon av 
%% vinkelen theta:
%% plot(theta,I8,theta,D);
plot(theta,I8);
grid on;
xlabel('theta (radianer)');
ylabel('Normert intensitet');
%% Skriver denne figuren til en eps-fil med navn oppgave3.eps
print('oppgave3.eps');
figure;
%% Deretter plotter vi intensitetsfordelingen som funksjon av 
%% vinkelen theta for storre verdi paa spaltebredden:
plot(theta2,I8_2,theta2,D2);
%% plot(theta,I8_2);
grid on;
xlabel('theta (radianer)');
ylabel('Normert intensitet');
%% Skriver denne figuren til en eps-fil med navn oppgave3.eps
print('oppgave3los.eps');