Off-axis Parabola

import numpy as np

from optiland import optic

# Initialize an optic system
system = optic.Optic()

# add surfaces to the system
system.surfaces.add(index=0, thickness=np.inf)  # object at infinity
system.surfaces.add(
    index=1,
    thickness=-100,
    radius=-200,
    conic=-1,
    material="mirror",
    is_stop=True,
    dy=36,
)
system.surfaces.add(index=2)

# add aperture
system.set_aperture(aperture_type="EPD", value=33)

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

# add wavelength
system.wavelengths.add(value=0.587, is_primary=True)

system.draw(num_rays=4)