function lens_abcd =abcd_lens(focal_distance)
% abcd matrix of the lens with given focal length
	lens_abcd = [1, 0; -1/focal_distance, 1];
endfunction