function spectrum=pk(spectrum,ph0,w)
	spectrum(:,2)=spectrum(:,2).*exp(i*ph0);
	plot(spectrum(:,1),real(spectrum(:,2)))
end
