1subroutine nuttal_window(win,n)
2  real win(n)
3
4  pi=4.0*atan(1.0)
5  a0=0.3635819
6  a1=-0.4891775;
7  a2=0.1365995;
8  a3=-0.0106411;
9  do i=1,n
10     win(i)=a0+a1*cos(2*pi*(i-1)/(n))+ &
11          a2*cos(4*pi*(i-1)/(n))+ &
12          a3*cos(6*pi*(i-1)/(n))
13  enddo
14  return
15end subroutine nuttal_window
16