Odd Asphere

import numpy as np

from optiland import optic

lens = optic.Optic()

# add surfaces
lens.surfaces.add(index=0, thickness=np.inf)
lens.surfaces.add(
    index=1,
    thickness=7,
    radius=56.128,
    material="N-SF11",
    is_stop=True,
    surface_type="odd_asphere",
    coefficients=[-3.142e-3, 3.867e-3, -1.602e-4],
)
lens.surfaces.add(index=2, thickness=30, radius=-58.697)
lens.surfaces.add(index=3)

# set aperture
lens.set_aperture(aperture_type="EPD", value=15)

# add field
lens.fields.set_type(field_type="angle")
lens.fields.add(y=0)

# add wavelength
lens.wavelengths.add(value=0.55, is_primary=True)

# draw lens
lens.draw(num_rays=10)