subroutine ruffen1( conj, n, pp, qq ) integer i, conj, n real pp(n), qq( n-1) call conjnull( conj, 0, pp,n, qq, n-1) if( conj == 0 ) do i= 1, n-1 qq(i) = qq(i) + pp(i+1) - pp(i) else do i= 1, n-1 { pp(i+1) = pp(i+1) + qq(i) pp(i ) = pp(i ) - qq(i) } return; end