p;&nbsp;&nbsp;&nbsp;&nbsp;t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;t
<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;p&nbsp;*(jac&nbsp;*jac)*p&nbsp;=&nbsp;r&nbsp;*r,
<P>
where <I>p</I> is a permutation matrix and <I>jac</I> is the final
calculated Jacobian. column <B>j</B> of <I>p</I> is column
<I>ipvt</I>(<B>j</B>) (see below) of the identity matrix. The lower
trapezoidal part of <I>fjac</I> contains information generated during
the computation o