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.