Sub MRQCOF(X(), Y(), SIG(), NDATA, A(), MA, LISTA(), MFIT, ALPHA(), BETA(), NALP, CHISQ)
    Dim DYDA(20)
    For J = 1 To MFIT
        For K = 1 To J
            ALPHA(J, K) = 0#
        Next K
        BETA(J) = 0#
    Next J
    CHISQ = 0#
    For I = 1 To NDATA
        Call FGAUSS(X(I), A(), YMOD, DYDA(), MA)
        SIG2I = 1# / (SIG(I) * SIG(I))
        DY = Y(I) - YMOD
        For J = 1 To MFIT
            WT = DYDA(LISTA(J)) * SIG2I
            For K = 1 To J
                ALPHA(J, K) = ALPHA(J, K) + WT * DYDA(LISTA(K))
            Next K
            BETA(J) = BETA(J) + DY * WT
        Next J
        CHISQ = CHISQ + DY * DY * SIG2I
    Next I
    For J = 2 To MFIT
        For K = 1 To J - 1
            ALPHA(K, J) = ALPHA(J, K)
        Next K
    Next J
End Sub