function b=factorial(a)
	b=1;	
	for i=1:a	
	    b=b*i;
	end