Path: | ZolotarevModule/sncndnK_kennedy_mod.F90 |
Last Update: | Thu Sep 29 15:16:05 +0900 2011 |
This file provides the following module.
Main Program : |
program et use etf implicit none real(QP) :: v,dv,q,dq,aa,bb integer :: i,ns v = -PI*2 dv = 0.1_QP dq = 0.001_QP do while(v < 2*PI) q = 0.0_QP do while(q < 1.00_QP) aa=etf4(v,q) bb=detf4(v,q) write(*,'(2F24.14,8E26.15E4)')v,q,aa,bb,abs(aa-bb)/abs(bb) q = q + dq enddo write(*,*) write(*,*) v = v + dv enddo stop end program