program test
  use inout
  use gamma
  implicit none
  integer                        ::N
  real,dimension(:),allocatable  ::u
  real                           ::t
  integer                        ::nit,i

  real                           ::beta,Td,precision,gam

  N=50
  
!-------------------------------------------------------!
!                                                       !
!         Les lignes suivantes testent l'ecriture       !
!                                                       !
!-------------------------------------------------------!

!!$  allocate(u(0:(N+1)))
!!$
!!$  call Init_ecriture(N)
!!$
!!$  u(:)=0.
!!$
!!$  !initialisation
!!$  do i=0,24
!!$     u(i)=1.
!!$  enddo
!!$
!!$  t=0.
!!$  nit=0
!!$
!!$  call ecriture(u,n,t,nit)
!!$
!!$  do nit=1,25
!!$     do i=(n+1),1,-1
!!$        print*,'i',i
!!$        u(i)=u(i-1)
!!$     enddo
!!$     t=t+0.1
!!$     call ecriture(u,n,t,nit)
!!$  enddo


!-------------------------------------------------------!
!                                                       !
!     Les lignes suivantes testent le calcul de gamma   !
!                                                       !
!-------------------------------------------------------!

  beta=100.
  Td=0.5
  precision=1.d-8


  call calc_gamma(beta,Td,precision,gam)

  print*,'gam=',gam

end program test
