SCULIB_FIT_MAKEBETA

calculates beta matrix for SCULIB_FIT_FUNCTION

Description:

This routine calculates the elements of the beta matrix for non-linear fitting as described in ‘Data Reduction and Error Analysis for the Physical Sciences’ by Bevington and Robinson, section 8.6 ‘the Marquardt Method’.

The routine does this by calling the routine SCULIB_FIT_DXISQ_DAJ for each element required. The routine will only execute if entered with good status.

Invocation

CALL SCULIB_FIT_MAKEBETA (XISQ_ROUTINE, N, A, BETA, STATUS)

Arguments

XISQ_ROUTINE (XISQ, N, A, STATUS) = EXTERNAL ROUTINE (Given)
routine to calculate chi-squared of current fit
N = INTEGER (Given)
the number of fit parameters
A (N) = DOUBLE PRECISION (Given)
the values of the fit parameters
BETA (N) = DOUBLE PRECISION (Returned)
the beta vector
STATUS = INTEGER (Given and returned)
global status

Copyright

Copyright ©1995,1996,1997,1998,1999 Particle Physics and Astronomy Research Council. All Rights Reserved.