Professor Diomar Cesar Lobao

Universidade Federal Fluminense-Volta Redonda, RJ, Brasil

Diomar Cesar


Dept. Ciências Exatas - Exact Science Dept.

Search

linint4.f

c<html>
c<body>
c<pre> 
      subroutine linint(x,y,xtab,ytab,ntab)
c
c     Given a value of x return a value of y based on interpolation
c     within a table of y values (ytab) corresponding to the x values
c     contained in the array xtab.  The subroutine assumes that the
c     values in xtab increase monotonically.  Efficiency is increased
c     by remembering the table points used in the last call (ilast)
c     xtab and ytab are provided through the argument list as is their
c     length ntab.
c
c    John Mahaffy 2/12/95
c
      real xtab(ntab),ytab(ntab)
c<a name="save"><font color="FF0000">
      save ilast
c</a></font>
      data ilast/1/
c
c    Start the search from the last point of table use index
c
      if (x.le.xtab(ilast+1)) then
c
c    Search down the table from point of last use
c
          do 20 i1=ilast,1,-1
              if(x.ge.xtab(i1)) go to 60
  20          continue
          write(6,*) 'x = ', x, '  is below the table range'  
          stop
      else
c
c    Search up the table from point of last use
c
          do 40 i1=ilast+1,ntab-1
              if(x.le.xtab(i1+1)) go to 60
  40          continue
          write(6,*) 'x = ', x, '  is above the table range'
          stop
      endif
c
c   Bounding points found, interpolate
c
  60  wx=(x-xtab(i1))/(xtab(i1+1)-xtab(i1))
      y=(1-wx)*ytab(i1)+wx*ytab(i1+1)
      ilast=i1
      return
      end
c</pre>
c</body>
c</html>
Skip to content