function y = cos_im2(lambday, n)
	tt=linspace(0, n, n+1);
	row1 = 1 + cos(2* pi * tt/ lambday);
	a1 = ones (n + 1, 1);
	y = a1 * row1;

