C --- Inner product with a single loop -----------------
C -------------------------------- version:  May 1996 --

	Double Precision FUNCTION DScapr0(N,X,Y)
	Double Precision X(*),Y(*),S
        Integer   N,I

	S = 0D0
	DO  10  I=1,N
	S = S + X(I) * Y(I)
10	CONTINUE
	DScapr0 = S
	RETURN
	END