code 34370;
    integer procedure QRIVALHRM(A, N, VAL, EM); value N;
    integer N; array A, VAL, EM;
    begin array B,BB[1:N];
        integer I;
        HSHHRMTRIVAL(A, N, VAL, BB, EM); B[N]:=BB[N]:= 0;
        for I:=1 step 1 until N-1 do B[I]:=SQRT(BB[I]);
        QRIVALHRM:=QRIVALSYMTRI(VAL, BB, N, EM)
    end QRIVALHRM;
        eop