program test_bsin
use fvn_fnlib
use Kind_Definition
implicit none
integer :: n,npoints,i
real(kind=sp_kind) :: xmin,xmax,xstep,x,y
real(kind=sp_kind),dimension(200,0:5) :: bes

open(2,file='bsin.dat')

xmin=-5.
xmax=5.
npoints=200
xstep=(xmax-xmin)/dble(npoints)
do i=1,npoints
      x=xmin+i*xstep
      do n=0,5
        bes(i,n)=bsin(n,x)
      end do
      write(2,'(7e22.14)') x,bes(i,0:5)
end do
close(2)

end program