Código matlab: desplazamiento de un haz en función del ángulo

Para su reporte,  pueden usar este código que reproduce las curvas de desplazamiento del haz en una ventana, en función del ángulo. 

%% Script para graficar tres curvas relacionadas con la desviacion de un haz
%% de luz al atrevesar una placa transparente, en fucion del ángulo de
%% incidencia

clc % limpia la pantalla
clear % limpia la memoria temporarl
close all % cierra todas las ventanas de figuras

theta =0:0.01:pi/2; % vector que contiene los grados en radianes
thetas = theta*90*2/pi; % vector en grados: desde cero hasta 90 grados, para
% graficar se entiende mejor
h = 1; % espesor unitario
N =1.3; % indice de refraccion de la primera muestra imaginaria
d = h*sin(theta).*(1 - (cos(theta)./(sqrt(N^2 - (sin(theta)).^2)) ));
%calculo de desplazamiento
plot(thetas, d*100, 'b'); % graficando, multiplico por 100 para tener
% una referencia porcentual, creo que se
% entiende mejor
hold on

N =2;
d = h*sin(theta).*(1 - (cos(theta)./(sqrt(N^2 - (sin(theta)).^2)) ));
plot(thetas, d*100, 'r');

N =4;
d = h*sin(theta).*(1 - (cos(theta)./(sqrt(N^2 - (sin(theta)).^2)) ));
plot(thetas, d*100, 'k');

hold off

% fin del guion matlab