subroutine linfit(x, y, ycalc, a, b, n)
c
c -- generate a straight line for x-y 
c
      integer n, i
      real x(1), y(1), ycalc(1), a, b
c
      a = 2
      b = 5
      do 10 i = 1, n
       ycalc(i) = a + b * x(i)
10    continue
      return
      end