PROGRAM TESTV1
C
C	  PROGRAM TO TEST FOLLOWING ROUTINES ON MATLIB
C	  VADD, VCLR, VCOM, VMAG, VMOV, VSCL, VSUB 
C	  (ALL ARE SINGLE PRECISION ROUTINES)
C
	REAL A(5),B(5),C(5),P,Q
	DATA P,Q/ 5.,2./
C
	DO 10 I=1,5
	A(I) = I
	B(I) = 2*I
10	CONTINUE
	TYPE *,' SINGLE PRECISON VECTOR TEST'
	TYPE 20,A,B,P,Q
20	FORMAT(' A=',5F5.1,/' B=',5F5.1,/' P='F5.1,' Q='F5.1)
	CALL VADD(C,A,B,4)
	TYPE *,' VADD ',C
	CALL VCLR(C,4)
	TYPE *,' VCLR ',C
	CALL VCOM(C,P,A,Q,B,4)
	TYPE *,' VCOM ',C
	C(1) = VMAG(A,4)
	TYPE *,' VMAG ',C(1)
	CALL VMOV(C,A,4)
	TYPE *,' VMOV ',C
	CALL VSCL(C,P,A,4)
	TYPE *,' VSCL ',C
	CALL VSUB(C,A,B,4)
	TYPE *,' VSUB ',C
	DO 30 I=1,5
 30	C(I) = 0.
	TYPE *,' 3-VECTOR TEST'
	CALL V3ADD(C,A(3),B(3))
	TYPE *,' V3ADD ',C
	CALL V3CLR(C)
	TYPE *,' V3CLR ',C
	CALL V3COM(C,P,A(3),Q,B(3))
	TYPE *,' V3COM ',C
	C(1) = V3MAG(A(3))
	TYPE *,' V3MAG ',C(1)
	CALL V3SCL(C,P,A(3))
	TYPE *,' V3SCL ',C
	CALL V3SUB(C,A(3),B(3))
	TYPE *,' V3SUB ',C
	STOP
	END