mod_odeint Module

This module packages odeint from numerical recipes.

Uses

- mod_comm_lib
Help

Graph Key

Nodes of different colours represent the following:

Solid arrows point from a submodule to the (sub)module which it is descended from. Dashed arrows point from a module or program unit to modules which it uses.

Used by

Help

Graph Key

Nodes of different colours represent the following:

Solid arrows point from a submodule to the (sub)module which it is descended from. Dashed arrows point from a module or program unit to modules which it uses.

Variables

Type	Visibility	Attributes		Name		Initial
integer,	public,	parameter	::	KMAXX	=	200
integer,	public,	parameter	::	MAXSTP	=	10000
integer,	public,	parameter	::	NMAX	=	50
double precision,	public,	parameter	::	TINY	=	1.d-30

Subroutines

public subroutine odeint(ystart, nvar, x1, x2, eps, h1, hmin, nok, nbad, derivs, rkqs, ierror)

Arguments

Type	Intent	Optional	Attributes		Name
double precision				::	ystart(nvar)
integer				::	nvar
double precision				::	x1
double precision				::	x2
double precision				::	eps
double precision				::	h1
double precision				::	hmin
integer				::	nok
integer				::	nbad
real				::	derivs
real				::	rkqs
integer,	intent(out)			::	ierror

public subroutine rkck(y, dydx, n, x, h, yout, yerr, derivs)

Arguments

Type	Intent	Optional	Attributes		Name
double precision				::	y(n)
double precision				::	dydx(n)
integer				::	n
double precision				::	x
double precision				::	h
double precision				::	yout(n)
double precision				::	yerr(n)
real				::	derivs

public subroutine rkqs(y, dydx, n, x, htry, eps, yscal, hdid, hnext, derivs)

Arguments

Type	Intent	Optional	Attributes		Name
double precision				::	y(n)
double precision				::	dydx(n)
integer				::	n
double precision				::	x
double precision				::	htry
double precision				::	eps
double precision				::	yscal(n)
double precision				::	hdid
double precision				::	hnext
real				::	derivs

public subroutine stop(text)

Arguments

Type	Intent	Optional	Attributes		Name
character(len=*),	intent(in)			::	text