# Materials
glass = Berreman4x4.IsotropicNonDispersiveMaterial(1.55)
front = back = Berreman4x4.IsotropicHalfSpace(glass)
# Liquid crystal oriented along the x direction
(no, ne) = (1.5, 1.6)
Dn = ne - no
LC = Berreman4x4.UniaxialNonDispersiveMaterial(no, ne)
R = Berreman4x4.rotation_v_theta(e_y, pi / 2)
LC = LC.rotated(R)
d = 4.33e-6
TN = Berreman4x4.TwistedMaterial(LC, d)
# Inhomogeneous layer
IL = Berreman4x4.InhomogeneousLayer(TN)
# Structure
s = Berreman4x4.Structure(front, [IL], back)
# Normal incidence:
Kx = 0.0
# Calculation parameters
(lbda_min, lbda_max) = (200e-9, 1) # (m)
k0_list = numpy.linspace(2 * pi / lbda_max, 2 * pi / lbda_min)
# Plot setup
fig = pyplot.figure()
ax = fig.add_subplot("111")
