compute sinuses signal

1 files changed, 25 insertions, 3 deletions

diff --git a/test.f90 b/test.f90
index abc2e58..b5beb0b 100644
--- a/test.f90
+++ b/test.f90
@@ -1,10 +1,32 @@
 PROGRAM test
 
   use, intrinsic :: iso_c_binding 
-  implicit none
+  IMPLICIT none
   
-  include 'fftw3.f03'
+  INCLUDE 'fftw3.f03'
 
-  write(*,*) "ok"
+  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
+
+  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