function qnew=q_after_abcd(q_old,abcd)
% calculates q parameter after propagating through element with given abcd matrix
	qnew=(q_old*abcd(1,1)+abcd(1,2))/(q_old*abcd(2,1)+abcd(2,2));
endfunction