dy1/dt routine.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=rk), | intent(in) | :: | t | |||
real(kind=rk), | intent(in) | :: | y(:) | |||
real(kind=rk), | intent(out) | :: | yprime(:) |
pure subroutine f(t, y, yprime) !! dy1/dt routine. real(rk), intent(in) :: t real(rk), intent(in) :: y(:) real(rk), intent(out) :: yprime(:) yprime(1) = ((2*log(y(1)) + 8.0_rk)/t - 5.0_rk)*y(1) end subroutine f