Errorbound.F90

Path: ZolotarevModule/Errorbound.F90
Last Update: Wed Jan 18 16:14:48 +0900 2012
dot/f_2.png

Required files

Methods

Included Modules

signfunc_mod

Public Instance methods

Main Program :

[Source]

program errorbound
  use signfunc_mod
  implicit none
  integer, parameter :: QP=KIND(1.0q0)
  integer :: nss,j
  real(QP) :: kappa,dkappa,eps
  do j=6,14,2
    eps=10.0_QP**(-j)
    kappa  = 1.0_QP/1000
    dkappa = exp(1.0_QP/1000)
    do while(kappa<1.0_QP)
      nss = get_ns_exact(kappa,eps)
      write(*,'(2ES24.15,2I3)')eps,kappa,nss
      kappa=kappa*dkappa
    enddo
    write(*,*)
    write(*,*)
  enddo
  stop
end program