dy/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) !! dy/dt routine. real(rk), intent(in) :: t real(rk), intent(in) :: y(:) real(rk), intent(out) :: yprime(:) yprime(1) = y(2) yprime(2) = 100*(one - y(1)**2)*y(2) - y(1) end subroutine f