code 31020;
    procedure MULVEC(L, U, SHIFT, A, B, X); value L,U,SHIFT,X;
    integer L,U,SHIFT; real X; array A,B;
    for L:= L step 1 until U do A[L]:= B[L+SHIFT]*X;
        eop