procedure fourier_transform(var x:complex_array_cpl;k:integer;delta:real);
var i,n:integer;
    sincos:complex_array_cpl;
begin
  n:=2 shl(k-1);
  calc_trigonometrics(sincos,k);
  if ft_direction_cpl=to_time then delta:=1.0/(n*delta);
  fft(x,sincos,k);
  for i:=1 to n do
     begin
       x[i].real:=x[i].real*delta;
       x[i].imag:=x[i].imag*delta;
     end;
end;
