Program g01tffe
! G01TFF Example Program Text
! Mark 29.1 Release. NAG Copyright 2023.
! .. Use Statements ..
Use nag_library, Only: g01tff, nag_wp
! .. Implicit None Statement ..
Implicit None
! .. Parameters ..
Integer, Parameter :: nin = 5, nout = 6
! .. Local Scalars ..
Real (Kind=nag_wp) :: tol
Integer :: i, ifail, la, lb, lout, lp, ltail
! .. Local Arrays ..
Real (Kind=nag_wp), Allocatable :: a(:), b(:), g(:), p(:)
Integer, Allocatable :: ivalid(:)
Character (1), Allocatable :: tail(:)
! .. Intrinsic Procedures ..
Intrinsic :: max, mod, repeat
! .. Executable Statements ..
Write (nout,*) 'G01TFF Example Program Results'
Write (nout,*)
! Skip heading in data file
Read (nin,*)
! Read in the tolerance
Read (nin,*) tol
! Read in the input vectors
Read (nin,*) ltail
Allocate (tail(ltail))
Read (nin,*) tail(1:ltail)
Read (nin,*) lp
Allocate (p(lp))
Read (nin,*) p(1:lp)
Read (nin,*) la
Allocate (a(la))
Read (nin,*) a(1:la)
Read (nin,*) lb
Allocate (b(lb))
Read (nin,*) b(1:lb)
! Allocate memory for output
lout = max(ltail,lp,la,lb)
Allocate (g(lout),ivalid(lout))
! Calculate deviates (inverse CDF)
ifail = -1
Call g01tff(ltail,tail,lp,p,la,a,lb,b,tol,g,ivalid,ifail)
If (ifail==0 .Or. ifail==1) Then
! Display titles
Write (nout,*) ' TAIL P A B G IVALID'
Write (nout,*) repeat('-',55)
! Display results
Do i = 1, lout
Write (nout,99999) tail(mod(i-1,ltail)+1), p(mod(i-1,lp)+1), &
a(mod(i-1,la)+1), b(mod(i-1,lb)+1), g(i), ivalid(i)
End Do
End If
99999 Format (5X,A,4X,F6.3,2(4X,F6.2),3X,F7.3,4X,I3)
End Program g01tffe