#include?
#include?
#include?

using?namespace?std;

const?double?PI?=?3.141592;

vector?>?DFT（vector?>&?theData）
{
	//?Define?the?Size?of?the?read?in?vector
	const?int?S?=?theData.size（）;

	//?Initalise?new?vector?with?size?of?S
	vector?>?out（S?0）;
	for（unsigned?i=0;?（i?	{
		out[i]?=?complex（0.0?0.0）;
		for（unsigned?j=0;?（j?		{
			out[i]?+=?theData[j]?*?polar（1.0?-?2?*?PI?*?i?*?j?/?S）;
		}
	}

	retur