
function d2ydx2 = dderiv(u)
% DDERIV calculates the second derivative in Fourier space
  n = length(u);
  k = make_k(n);
  d2ydx2 = real(ifft(-(k.^2).*fft(u)));
end
