Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=8), | intent(in) | :: | x | |||
real(kind=8), | intent(in), | optional | :: | t1 |
pure function dens_hyperc(x,t1) real(8),optional,intent(in) :: t1 real(8),intent(in) :: x REAL(8):: dens_hyperc,t1_,pi2,sqrt2 pi2=2.d0*acos(-1.d0) sqrt2=sqrt(2.d0) t1_=sqrt2 ; if(present(t1))t1_=t1 dens_hyperc = (1/(t1_*sqrt(pi2)))*exp(-(x**2)/(2.d0*t1_**2)) return end function dens_hyperc