C23456
      FUNCTION fperiodique(x)

      REAL r,pi
      INTEGER i,N
      PARAMETER(N=25)
      PARAMETER(pi=3.14159)

      r=0.
      do i=1,N
        r=r+(sin(2.*i*pi/N) - sin(2.*x*i*pi/N))**2
      enddo
      fperiodique=r
      END