! D05AAF Example Program Text ! Mark 24 Release. NAG Copyright 2012. Module d05aafe_mod ! D05AAF Example Program Module: ! Parameters and User-defined Routines ! .. Use Statements .. Use nag_library, Only: nag_wp ! .. Implicit None Statement .. Implicit None ! .. Parameters .. Real (Kind=nag_wp), Parameter :: a = 0.0_nag_wp Real (Kind=nag_wp), Parameter :: b = 1.0_nag_wp Real (Kind=nag_wp), Parameter :: lambda = 1.0_nag_wp Real (Kind=nag_wp), Parameter :: xval = 0.1_nag_wp Integer, Parameter :: ind = 2, n = 5, nout = 6 Integer, Parameter :: ldw1 = n Integer, Parameter :: ldw2 = 2*n + 2 Contains Function k1(x,s) ! .. Function Return Value .. Real (Kind=nag_wp) :: k1 ! .. Scalar Arguments .. Real (Kind=nag_wp), Intent (In) :: s, x ! .. Executable Statements .. k1 = s*(1.0_nag_wp-x) Return End Function k1 Function k2(x,s) ! .. Function Return Value .. Real (Kind=nag_wp) :: k2 ! .. Scalar Arguments .. Real (Kind=nag_wp), Intent (In) :: s, x ! .. Executable Statements .. k2 = x*(1.0_nag_wp-s) Return End Function k2 Function g(x) ! .. Use Statements .. Use nag_library, Only: x01aaf ! .. Function Return Value .. Real (Kind=nag_wp) :: g ! .. Scalar Arguments .. Real (Kind=nag_wp), Intent (In) :: x ! .. Local Scalars .. Real (Kind=nag_wp) :: pi ! .. Intrinsic Procedures .. Intrinsic :: sin ! .. Executable Statements .. pi = x01aaf(pi) g = sin(pi*x)*(1.0_nag_wp-1.0_nag_wp/(pi*pi)) Return End Function g End Module d05aafe_mod Program d05aafe ! D05AAF Example Main Program ! .. Use Statements .. Use nag_library, Only: c06dcf, d05aaf, nag_wp Use d05aafe_mod, Only: a, b, g, ind, k1, k2, lambda, ldw1, ldw2, n, & nout, xval ! .. Implicit None Statement .. Implicit None ! .. Local Scalars .. Integer :: ifail, is ! .. Local Arrays .. Real (Kind=nag_wp) :: ans(1), x(1) Real (Kind=nag_wp), Allocatable :: c(:), f(:), w1(:,:), w2(:,:), & wd(:) ! .. Executable Statements .. Write (nout,*) 'D05AAF Example Program Results' Allocate (c(n),f(n),w1(ldw1,ldw2),w2(ldw2,4),wd(ldw2)) ifail = 0 Call d05aaf(lambda,a,b,k1,k2,g,f,c,n,ind,w1,w2,wd,ldw1,ldw2,ifail) Write (nout,99999) Write (nout,99998) c(1:n) x(1) = xval Select Case (ind) Case (1) is = 3 Case (2) is = 2 Case Default is = 1 End Select ifail = 0 Call c06dcf(x,1,a,b,c,n,is,ans,ifail) Write (nout,99997) 'X=', x, ' ANS=', ans 99999 Format (/1X,'Kernel is centro-symmetric and G is even so the ', & 'solution is even'//1X,'Chebyshev coefficients'/) 99998 Format (1X,5E14.4/) 99997 Format (1X,A,F5.2,A,1F10.4) End Program d05aafe