hskl.f

来自「gsac程序包」· F 代码 · 共 21 行

	subroutine hskl(cossh,rsinsh,sinshr,TL,iwat,hl)
	implicit none
	complex*16 cossh, rsinsh, sinshr
	real TL
	integer iwat
	complex*16 hl(2,2)
	if(iwat.eq.0)then	
		hl(1,1) = cossh
		hl(2,1) = TL*rsinsh
		hl(1,2) = sinshr/TL
		hl(2,2) = cossh
	else
		hl(1,1) = dcmplx(1.0d+00,0.0d+00)
		hl(1,2) = dcmplx(0.0d+00,0.0d+00)
		hl(2,1) = dcmplx(0.0d+00,0.0d+00)
		hl(2,2) = dcmplx(1.0d+00,0.0d+00)
	endif
	return
	end