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

clear all;
%Konstanten
c = 3e8
c = 300000000
%Gegeben
r = [0.1 -0.08 0.08 0.04]
r = 1×4
    0.1000   -0.0800    0.0800    0.0400
n = [1 1.3 1.38 1.4 1.36]
n = 1×5
    1.0000    1.3000    1.3800    1.4000    1.3600
g = 2
g = 2

Brechkräfte

deltaN = diff(n)
deltaN = 1×4
    0.3000    0.0800    0.0200   -0.0400
D = deltaN ./ r 
D = 1×4
    3.0000   -1.0000    0.2500   -1.0000
Dtot = sum(D)
Dtot = 1.2500

Brennweiten

f_g = n(1) / Dtot
f_g = 0.8000
f_b = n(end) / Dtot
f_b = 1.0880

Bildweite

b = n(end) / (Dtot - n(1)/g)
b = 1.8133