Wie rechne ich... ein optisches System mit mehreren Grenzflächen?

clear all;
%Konstanten
c = 3e8
c = 300000000
%Gegeben
n0 = 1
n0 = 1
n1 = 1.5
n1 = 1.5000
n2 = 1.2
n2 = 1.2000
n3 = 1.1
n3 = 1.1000
n4 = 1.8
n4 = 1.8000
alpha0 = 30
alpha0 = 30

Winkel

alpha1 = asind(n0*sind(alpha0)/n1)
alpha1 = 19.4712
beta = 90 - alpha1
beta = 70.5288
Thetak = asind(n2/n1)
Thetak = 53.1301
Thetak < beta
ans = logical
1
gamma1 = alpha1
gamma1 = 19.4712
gamma2 = asind(n1*sind(gamma1)/n3)
gamma2 = 27.0357
delta1 = gamma2
delta1 = 27.0357
delta2 = asind(n3*sind(delta1)/n4)
delta2 = 16.1276
epsi1 = delta2
epsi1 = 16.1276
epsi2 = asind(n4 * sind(epsi1) / n0)
epsi2 = 30.0000