function z = cosdiag(freq)
	COS = cos_im2(freq, 1024);
	COS = rot45(COS, 1025);
	z = COS(600:855, 1:256);
