code 34431; procedure QUANEWBND1(N, LW, RW, X, F, FUNCT, IN, OUT); value N, LW, RW; integer N, LW, RW; array X, F, IN, OUT; boolean procedure FUNCT; begin integer I, K; real S; array JAC[1:(LW + RW) * (N - 1) + N]; FUNCT(N, 1, N, X, F); S:= IN[5]; K:= (LW + RW)*(N - 1) + N*2 - ((LW - 1)*LW + (RW - 1)*RW) // 2; IN[4]:= IN[4] - K; JACOBNBNDF(N, LW, RW, X, F, JAC, I, S, FUNCT); QUANEWBND(N, LW, RW, X, F, JAC, FUNCT, IN, OUT); IN[4]:= IN[4] + K; OUT[3]:= OUT[3] + K end QUANEWBND1; eop