Program s18aefe ! S18AEF Example Program Text ! Mark 24 Release. NAG Copyright 2012. ! .. Use Statements .. Use nag_library, Only: nag_wp, s18aef ! .. Implicit None Statement .. Implicit None ! .. Parameters .. Integer, Parameter :: nin = 5, nout = 6 ! .. Local Scalars .. Real (Kind=nag_wp) :: x, y Integer :: ifail, ioerr ! .. Executable Statements .. Write (nout,*) 'S18AEF Example Program Results' ! Skip heading in data file Read (nin,*) Write (nout,*) Write (nout,*) ' X Y' Write (nout,*) data: Do Read (nin,*,Iostat=ioerr) x If (ioerr<0) Then Exit data End If ifail = -1 y = s18aef(x,ifail) If (ifail<0) Then Exit data End If Write (nout,99999) x, y End Do data 99999 Format (1X,1P,2E12.3) End Program s18aefe