function dphi=dtkdv(t,phihat)

global k k3

%k = 1i*kx; k3=1i*kx^3
e=exp(k3*t);inve=exp(-k3*t);

uhat=e.*phihat;
u=real(ifft(uhat));
%u2=u;
u2=u.*u;
u2f=fft(u2);

dphi=-(k/2).*(u2f.*inve);