function [v] = dctus(u)
%% Gives a better 8x8 1-D DCT

for k=1:8
	for n=1:8
		if (k==1)
			c(k)= 1/sqrt(2)/2;
		else
			c(k)=1/2;
		end
		temp(n)=u(n)* cos( (pi* (2*(n-1)+1) * (k-1) /(2*8) ));
	end
	v(k) = c(k)*sum(temp);
end
