Real, Allocatable :: X(:,:), Y(:), X1(:),Beta(:)
     Integer :: N,M,i,j,k
     Real :: S,Q,U,RR,F

     Open (1,file="data_reg_test.dat",status="Old")
     Read (1,*) N,M
     Allocate (X(N,M-1),Y(N),X1(N*M))
     Do i=1,N
        Read (1,*) Y(i), (X(i,j),j=1,M-1)
     Enddo
     ! Convert to 1D-array:
     Do i=1,N; X1(i)=Y(i); Enddo   ! Y
     k=N
     Do j=1,M-1
        Do i=1,N
           k=k+1; X1(k)=X(i,j)
        Enddo
     Enddo

     Allocate (Beta(M))
     Call REGRE(X1,N,M,Beta, S,Q,U,RR,F)

     print"(10F10.3)", Beta
     print"(A,10F10.3)", "R & F: ", RR,F
     CONTAINS

!     END
     Include "/home/tanpv/Mylibs/Sub_HQTT.f90"
     END