Path: | ZolotarevModule/signfunc_mod.F90 |
Last Update: | Wed Dec 15 15:47:40 +0900 2010 |
This file provides the following module.
Main Program : |
program errorbound implicit none integer :: nss,j real(QP) :: kappa,dkappa,eps real(QP) :: comp_elliptic_integ_K do j=6,14 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_exact(kappa,eps) write(*,'(2ES24.15,2I3)')eps,kappa,nss kappa=kappa*dkappa enddo write(*,*) write(*,*) enddo stop end program