implicit double precision (a-h,o-z)
      dimension x(100),y(100),sigmay(100),a(10)
1     type 100
      accept 101,npts
      if (npts) 999,999,2
2     type 102
      accept 103,nterms
      type 104
      accept 105,mode
      do 10 i=1,npts
10    read (21,107) x(i),y(i),sigmay(i)
15    call polfit(x,y,sigmay,npts,nterms,mode,a,chisqr)
      type 200,chisqr
      do 30 i=1,nterms
30    type 201,i,a(i)
      type 203
999   stop
100   format(' number of points='$)
101   format(i)
102   format(' nterms='$)
103   format(i)
104   format(' mode='$)
105   format(i)
107   format(3d)
108   format('0input initial parameters:'/)
109   format(d)
110   format(' error='$)
200   format(' chisqr=',1pd20.8)
201   format(1x,i,1pd20.8)
203   format('1')
      end