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

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