double precision function fun (x)
      double precision x, dsin
      if (x .eq. 0.0) fun = 1.0
      if (x .ne. 0.0) fun = dsin(x)/x
      return
      end

.