1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
module oumods
use, intrinsic :: iso_c_binding
implicit integer(c_int) (i-k), integer(c_int) (m,n), &
& real(c_double) (a-h), real(c_double) (l), real(c_double) (o-z)
TYPE :: ex_type
INTEGER :: A = 0
CONTAINS
FINAL :: del_ex_type
PROCEDURE :: sub => ex_sub
END TYPE ex_type
contains
subroutine zI12(t,c,alpha,beta,r)
complex(c_double_complex) c,r, x,y,z
z = c*t
y = exp(z)
x = (2.0_c_double * cosh((z - cmplx(0._c_double,3.14159265358979324_c_double, kind(1._c_double))) &
& /2._c_double )) / (c / exp((z + cmplx(0._c_double,3.14159265358979324_c_double,kind(1._c_double)))/2._c_double))
r = beta*r+alpha*((t*y - x)/c)
end subroutine
end module
|
