1 2 3 4 5 6
function lens_abcd =abcd_lens(focal_distance) % abcd matrix of the lens with given focal length lens_abcd = [1, 0; -1/focal_distance, 1]; end