dg <- function(N) {
  step <- 2.5/N
  beta <- (function(x)(x*exp(-x^2)))(step*1:N)
  beta <- c(-beta[N+1-1:N],0,beta)
  beta <- - beta/(1.4*2*pi)*step^2*10
  beta
}