Name: DPOLYN
Purpose: DPOLYN computes the partial derivatives of that polynomial
Files: dpolyn.dc2,polyn.shl
Author: Peter Roelfsema
Call:
CALL DPOLYN(X,PAR,DERIV,NPAR) ( entry point in POLYN )
Input parameters:
X -R -running coordinate.
PAR -R -array defined as follows: PAR(0:NPAR) that
contains the parameters for the polynomial:
PAR(n)-polynomial coefficient for X**n.
NPAR -I -degree of the polynomial.
Output parameters:
DERIV -R -array defined as follows: DERIV(0:NPAR) that
contains the partial derivatives of the
polynomial with respect to the parameters:
DERIV(n)-partial derivative of coefficient for
X**n ( i.e. X**n ).
Use:
DPOLYN gives the partial derivatives for the gaussian.
Note that the polynomial is defined as the following
summation:
PX = 0 FOR N=0 TO NPAR ( PX = PX + PAR(N) * X**N )
DPOLYN does not check for floating overflows in X**N, you
will have to worry about that yourself.