function diff(x,y)
c     c.l.lawson and r.j.hanson, jet propulsion laboratory, 1973 june 7
c     to appear in 'solving least squares problems', prentice-hall, 1974
      diff=x-y
      return
      end