blob: 38f360f76f01bc5940c6471573e3d4f29f281403 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
PROGRAM test
use, intrinsic :: iso_c_binding
IMPLICIT none
INCLUDE 'fftw3.f03'
integer, parameter :: npoints = 1024
real(C_DOUBLE), parameter :: start = 0, end = 1
integer :: i
real(C_DOUBLE) :: t, &
delta = (end - start) / real(npoints - 1, C_DOUBLE)
DO i = 0, npoints - 1
t = i * delta + start
write(*,*) t, " ", f1(t)
END DO
CONTAINS
PURE real(C_DOUBLE) FUNCTION f1(t)
real(C_DOUBLE), intent(in) :: t
real(C_DOUBLE), parameter :: pi = acos(-1.0)
f1 = sin(2 * pi * t * 200) + 2 * sin(2 * pi * t * 400)
END FUNCTION f1
END PROGRAM test
|