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