Η κλίση της αντικειμενικής συνάρτησης.Το υποπρόγραμμα: SUBROUTINE GRANAL
subroutine granal ( n, p, grad )
implicit double precision (a-h,o-z)
parameter(mx=8000,nx=250)
common/xydata/x(mx),y(mx),mm
dimension p(n),grad(n),ff(mx),g(nx)
call subsum ( mm, n, p, ff )
call gann(p,nodes,x(l),g)
grad(j)= grad(j) + 2*ff(l)*g(j)