c-----------------------------------------------------------------------
	subroutine ARBINT(n,x,y,yint)
c finds the integral of a function on
c a non-uniformly spaced set of points, using two-point formulas
c
c author: Kieron Burke, 9.15.99
c     in: n - no of points
c	  x - vector of x-values
c	  y - vector of y-values
c    out: yint - integral
c  calls: --
	implicit real*8(a-h,o-z)
	parameter(m=1000000)
	dimension x(0:m),y(0:m),yint(0:m)