Υπολογισμός της Ιακωβιανής.Το υποπρόγραμμα: SUBROUTINE JANAL
subroutine janal ( m, n, p, fj, ld )
implicit double precision (a-h,o-z)
parameter(mx=8000,nx=250)
common/xydata/x(mx),y(mx),mm
dimension p(n), fj(ld,n),gper(nx)
call gann(p,nodes,x(i),gper)