function [a] = gen_a(theta, par) a = exp(-1j*2*pi*par.delta_d*sind(theta)/par.labmda); end