SUBROUTINE FJACS(X,QR,LENQR,PIVOT)
C INTEGER LENQR,N,PIVOT(N+2)
C DOUBLE PRECISION QR(LENQR),X(N)
C
C Evaluate the N x N symmetric Jacobian matrix of F(X) at X, and return
C the result in packed skyline storage format in QR. LENQR is the length
C of QR, and PIVOT contains the indices of the diagonal elements of the
C Jacobian matrix within QR. PIVOT(N+1) and PIVOT(N+2) are set by
C subroutine FODE.
C
RETURN
END