From f894d4ea34045c97aea2fc0853696431797d9f81 Mon Sep 17 00:00:00 2001
From: Wojtek Kosior <kwojtus@protonmail.com>
Date: Sat, 1 Jun 2019 17:27:42 +0200
Subject: plot original function

---
 Makefile                | 6 ++++++
 src/f1.f90              | 2 +-
 src/f1_original.gnuplot | 4 ++++
 3 files changed, 11 insertions(+), 1 deletion(-)
 create mode 100644 src/f1_original.gnuplot

diff --git a/Makefile b/Makefile
index 772ca48..8c77b48 100644
--- a/Makefile
+++ b/Makefile
@@ -7,6 +7,12 @@ all : f1
 f1 : src/f1.f90
 	$(FC) $(FFLAGS) src/f1.f90 -o $@
 
+res/f1_original.txt : f1
+	./f1 > res/f1_original.txt
+
+res/%.png : src/%.gnuplot res/%.txt
+	gnuplot -c $<
+
 clean :
 	rm f1
 
diff --git a/src/f1.f90 b/src/f1.f90
index e1e621e..7f07afa 100644
--- a/src/f1.f90
+++ b/src/f1.f90
@@ -8,7 +8,7 @@ PROGRAM test
   real(kind=8), parameter :: pi = acos(-1.0) ! also used from f1()
 
   ! rest of variables only used from "main" program code
-  integer, parameter :: npoints = 1024, ncomplex = npoints / 2 + 1
+  integer, parameter :: npoints = 65536, ncomplex = npoints / 2 + 1
   real(kind=8), parameter :: start = 0, end = 1
   integer :: i
   real(kind=8) :: t,                                            &
diff --git a/src/f1_original.gnuplot b/src/f1_original.gnuplot
new file mode 100644
index 0000000..c0ecba8
--- /dev/null
+++ b/src/f1_original.gnuplot
@@ -0,0 +1,4 @@
+set terminal png size 1440, 900
+set output 'res/f1_original.png'
+set xrange [0:0.15]
+plot 'res/f1_original.txt' title 'original' with lines
\ No newline at end of file
-- 
cgit v1.2.3