code 34437;
procedure JACOBNNF(N, X, F, JAC, I, DI, FUNCT); value N;
integer N, I; real DI; array X, F, JAC; procedure FUNCT;
begin integer J; real STEP, AID; array F1[1:N];
    for I:= 1 step 1 until N do 
    begin STEP:= DI; AID:= X[I]; X[I]:= AID + STEP;
        STEP:= 1 / STEP; FUNCT(N, X, F1);
        for J:= 1 step 1 until N do 
        JAC[J,I]:= (F1[J] - F[J]) * STEP; X[I]:= AID
    end 
end JACOBNNF;
        eop