Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=8), | intent(in) | :: | x | |||
real(kind=8), | intent(in), | optional | :: | ts |
function dens_2dsquare(x,ts) result(dos) real(8),intent(in) :: x real(8),intent(in),optional :: ts real(8) :: wband,y,kint,eint,dos real(8),parameter :: pi=acos(-1d0) wband=4.d0;if(present(ts))wband=4.d0*ts dos=0d0 if(abs(x)>wband)return y = sqrt(1d0 - (x/wband)**2) dos= 2/pi**2/wband*EllipticK(y) end function dens_2dsquare