Blame view

fvn_test/test_bskn.f90 439 Bytes
f26a262db   daniau   git-svn-id: https...
1
2
3
4
  program test_bskn
  use fvn_fnlib
  implicit none
  integer :: n,npoints,i
f6bacaf83   cwaterkeyn   ChW 11/09: ANSI c...
5
6
  real(kind=sp_kind) :: xmin,xmax,xstep,x,y
  real(kind=sp_kind),dimension(200,0:5) :: bes
fdbb9bb0a   daniau   git-svn-id: https...
7

f26a262db   daniau   git-svn-id: https...
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
  open(2,file='bskn.dat')
  
  xmin=0.
  xmax=10.
  npoints=200
  xstep=(xmax-xmin)/dble(npoints)
  do i=1,npoints
        x=xmin+i*xstep
        do n=0,5
          bes(i,n)=bskn(n,x)
        end do
        write(2,'(7e22.14)') x,bes(i,0:5)
  end do
  close(2)
  
  ! with gnuplot -> set yrange [-1:1]
  
  end program