ADAS Subroutine h4fasy
C subroutine h4fasy( istdim , & x , xa , n , ya , y , dy , & c1 , c2 , c3 , c4 , & form , iforms & ) c----------------------------------------------------------------------- c c **************** fortran77 subroutine: h4spl3.for ******************* c c purpose: provide a spline interpolate making use of specified c asymptotic behaviour c c calling program: various c c notes: (1) uses labelled common /espl3/ c c input : (i*4) istdim = dimensionality for splining arrays c c input : (r*8) x = required x-value c input : (i*4) x(i) = knots c input : (r*8) n = number of knots c c input : (r*8) c1(i,j) = 1st spline coefficient precursor c input : (r*8) c2(i,j) = 2nd spline coefficient precursor c input : (r*8) c3(i,j) = 3rd spline coefficient precursor c input : (r*8) c4(i,j) = 4th spline coefficient precursor c input : (r*8) form = external function specifying asymptotic forms c input : (i*4) iforms = index of required form c c output: (r*8) y=returned y-value c output: (r*8) dy=returned derivative c c c routines: c routine source brief description c ------------------------------------------------------------- c h4fspl adas c i4unit adas fetch unit number for output of messages c c author: Hugh P. Summers, University of Strathclyde c JA7.08 c Tel.: +44 (0)141-548-4196 c c date: 24 July 2002 c c c version: 1.1 Hugh Summers 24/07/02 c modified: first release c c----------------------------------------------------------------------- INTEGER IFORMS, ISTDIM, N REAL*8 C1(ISTDIM,ISTDIM-1), C2(ISTDIM,ISTDIM-1) REAL*8 C3(ISTDIM,ISTDIM-1), C4(ISTDIM,ISTDIM-1) REAL*8 DY, X, XA(ISTDIM), Y REAL*8 YA(ISTDIM)