Bisection step.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=wp), | intent(in) | :: | x1 | |||
real(kind=wp), | intent(in) | :: | x2 |
point half way between x1 and x2
pure function bisect(x1,x2) result(x3) implicit none real(wp),intent(in) :: x1,x2 real(wp) :: x3 !! point half way between x1 and x2 x3 = (x1 + x2) / 2.0_wp end function bisect