gfortran -o f90_dgint3 f90_dgint3.f90
gfortran -o f90_dsimp f90_dsimp.f90

./f90_dgint3 1 1 > _out_g.txt
./f90_dgint3 2 2 >> _out_g.txt
./f90_dgint3 3 3 >> _out_g.txt
./f90_dgint3 4 4 >> _out_g.txt
./f90_dgint3 5 5 >> _out_g.txt
./f90_dgint3 7 7 >> _out_g.txt
./f90_dgint3 10 10 >> _out_g.txt
./f90_dgint3 20 20 >> _out_g.txt
./f90_dgint3 30 30 >> _out_g.txt
./f90_dgint3 40 40 >> _out_g.txt
./f90_dgint3 50 50 >> _out_g.txt
./f90_dgint3 70 70 >> _out_g.txt
./f90_dgint3 100 100 >> _out_g.txt
./f90_dgint3 200 200 >> _out_g.txt
./f90_dgint3 300 300 >> _out_g.txt
./f90_dgint3 400 400 >> _out_g.txt
./f90_dgint3 500 500 >> _out_g.txt
./f90_dgint3 700 700 >> _out_g.txt
./f90_dgint3 1000 1000 >> _out_g.txt

./f90_dsimp 1 1 > _out_s.txt
./f90_dsimp 2 2 >> _out_s.txt
./f90_dsimp 3 3 >> _out_s.txt
./f90_dsimp 4 4 >> _out_s.txt
./f90_dsimp 5 5 >> _out_s.txt
./f90_dsimp 7 7 >> _out_s.txt
./f90_dsimp 10 10 >> _out_s.txt
./f90_dsimp 20 20 >> _out_s.txt
./f90_dsimp 30 30 >> _out_s.txt
./f90_dsimp 40 40 >> _out_s.txt
./f90_dsimp 50 50 >> _out_s.txt
./f90_dsimp 70 70 >> _out_s.txt
./f90_dsimp 100 100 >> _out_s.txt
./f90_dsimp 200 200 >> _out_s.txt
./f90_dsimp 300 300 >> _out_s.txt
./f90_dsimp 400 400 >> _out_s.txt
./f90_dsimp 500 500 >> _out_s.txt
./f90_dsimp 700 700 >> _out_s.txt
./f90_dsimp 1000 1000 >> _out_s.txt

source a_gmt.txt