move: text.f90 → src/f1.f90

1 files changed, 32 insertions, 0 deletions

diff --git a/src/f1.f90 b/src/f1.f90
new file mode 100644
index 0000000..38f360f
--- /dev/null
+++ b/src/f1.f90
@@ -0,0 +1,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