code 34438;
procedure JACOBNMF(N, M, X, F, JAC, I, DI, FUNCT); value N, M;
integer N, M, 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 M do 
    begin STEP:= DI; AID:= X[I]; X[I]:= AID + STEP;
        STEP:= 1 / STEP; FUNCT(N, M, X, F1);
        for J:= 1 step 1 until N do 
        JAC[J,I]:= (F1[J] - F[J]) * STEP; X[I]:= AID
    end 
end JACOBNMF;
        eop