13 - SPARSE LINEARLY CONSTRAINED NONLINEAR PROGRAMMING
------------------------------------------------------


      CALL EILS13(N,NC,MROWS,NZA,X,PTR,CNR,AU,LENR,
     &            IX,XL,XU,IC,CL,CU,MFP,NEXT,IERR)


         N - INITIAL NUMBER OF NUMBER OF VARIABLES (INPUT/OUTPUT -
             THIS DESCRIPTION MEANS THAT WE CAN SPECIFY IT PRECISELY
             EITHER IN THE INPUT STREAM OR HERE - WE CAN REDEFINE ITS
             VALUE IN THIS SUBROUTINE. THIS POSSIBILITY IS ENABLED BY
             THE SEPARATE SIZES OF DECLARED WORKING NUMBER OF
             ROWS/COLUMNS)
         NC - INITIAL NUMBER OF CONSTRAINTS (INPUT/OUTPUT)
         MROWS - DECLARED NUMBER OF VARIABLES (INPUT)
         NZA - SIZE OF ARRAYS AU/CNR (INPUT)
         X(N) - VECTOR OF VARIABLES (OUTPUT)
         PTR(MROWS) - VECTOR SHOWING BEGINNINGS OF THE INPUT MATRIX ROWS
                (OUTPUT)
         CNR(NZA) - INPUT MATRIX COLUMN INDICES (OUTPUT)
         AU(NZA) - INPUT MATRIX VALUES (OUTPUT)
         LENR(MROWS) - LENGTHS OF THE INPUT MATRIX ROWS (OUTPUT)
         IX(NF) - VECTOR CONTAINING TYPES OF BOUNDS (OUTPUT):
                 =1 FOR LOWER BOUNDED VARIABLE.
                 =2 FOR UPPER BOUNDED VARIABLE.
                 =3 FOR LOWER AND UPPER BOUNDED ONES.
         XL(NF) - VECTOR CONTAINING LOWER BOUNDS FOR VARIABLES (OUTPUT)
         XU(NF) - VECTOR CONTAINING UPPER BOUNDS FOR VARIABLES (OUTPUT)
         IC(NC) - VECTOR CONTAINING TYPES OF CONSTRAINTS (OUTPUT):
                 =1 FOR LOWER BOUNDED ROW.
                 =2 FOR UPPER BOUNDED ROW.
                 =3 FOR LOWER AND UPPER BOUNDED ROW.
         CL(NC) - VECTOR CONTAINING LOWER BOUNDS FOR CONSTRAINT FUNCTIONS.
         CU(NC) - VECTOR CONTAINING UPPER BOUNDS FOR CONSTRAINT FUNCTIONS.
         MFP -  VARIABLE SHOWS WHETHER INITIAL SOLUTION IS(>1) OR IS NOT
                (=0) GIVEN.
         NEXT - NO. OF PROBLEM (INPUT): THERE ARE 6 PROBLEMS
                 (NEXT FROM 1 TO 6)
         IERR - CONFLICT IN INPUT DATA IF IERR <> 0 (OUTPUT)



      CALL EINS13(MCOLS,N,M,HF,JH,IH,NEXT,IERR)


         MCOLS - DECLARED NUMBER OF VARIABLES (INPUT)
         N - NUMBER OF VARIABLES (INPUT/OUTPUT)
         M - SIZE OF WORK SPACE FOR SPARSE HESSIAN (INPUT)
         HF(M) - NUMERICAL VALUES OF THE HESSIAN MATRIX OF THE
                 QUADRATIC MODEL (UNNECESSARY IN MOST CASES)
         JH(M) - INDEX VECTOR OF THE HESSIAN MATRIX (OUTPUT)
         IH(MCOLS+1) - POINTER VECTOR FOR THE ROWS OF HESSIAN (OUTPUT)
         NEXT - NUMBER OF THE TEST PROBLEM (INPUT)
         IERR - CONFLICT IN INPUT DATA IF IERR <> 0 (OUTPUT)



      CALL EFFU13(N,X,F,NEXT)

         N - NUMBER OF VARIABLES (INPUT)
         X(N) - ACTUAL VECTOR OF VARIABLES (INPUT)
         F - VALUE OF THE OBJECTIVE FUNCTION (OUTPUT)
         NEXT - NO. OF PROBLEM (INPUT)


      CALL EFGU13(N,X,G,NEXT)

         N - NUMBER OF VARIABLES (INPUT)
         X(N) - ACTUAL VECTOR OF VARIABLES (INPUT)
         G(N) - GRADIENT OF THE OBJECTIVE FUNCTION (OUTPUT)
         NEXT - NO. OF PROBLEM (INPUT)

