Solves a first order differential-algebraic system of equations, g(t, y, yʹ) = 0, possibly with additional constraints.

Approximates the Jacobian of m functions in n unknowns using divided differences

Solves a linear programming problem.

Computes the first, second or third derivative of a user-supplied function.

Computes the determinant of a rectangular matrix, A.

Constructs a square diagonal matrix from a rank-1 array or several diagonal matrices from a rank-2 array.

Extracts a rank-1 array whose values are the diagonal terms of a rank-2 array argument.

Computes the 1-norm distance between two points.

Computes the Euclidean (2-norm) distance between two points.

Computes the infinity norm distance between two points.

Solves a linear programming problem via the revised simplex algorithm.

See AMACH.

DQADD (See Extended Precision Arithmetic, Chapter 9)

Adds a double-precision scalar to the accumulator in extended precision.

DQINI (See Extended Precision Arithmetic, Chapter 9)

Initializes an extended-precision accumulator with a double-precision scalar.

DQMUL (See Extended Precision Arithmetic, Chapter 9)

Multiplies double-precision scalars in extended precision.

DQSTO (See Extended Precision Arithmetic Chapter 9)

Stores a double-precision approximation to an extended-precision scalar.

DSDOT (See Chapter 9)

Computes the single-precision dot product xTy using a double precision accumulator.

This routine handles MATH/LIBRARY and STAT/LIBRARY type DOUBLE PRECISION options.