implicit real*8(a-h,o-z)
      parameter(imax=101,iimax=100)
      dimension asig(imax),as0i(imax),dist(iimax)
*
      inpts= 7
      inrand= 7
      wmi= 80.39d0
      rsi= 186.d0
      alszi= 0.1194d0
      thc= 10.d0
      s0i= 0.01d0
      tc= 0.5d0
      tc= tc/180.d0*3.141592653589793238462643d0
      call ww(inpts,inrand,wmi,rsi,s0i,thc,tc,alszi,sig,esig)
      re= 1.d2*esig/sig
      print 1,sig,esig
      print 2,re
      thc= 8.d0
      call ww(inpts,inrand,wmi,rsi,s0i,thc,tc,alszi,sig,esig)
      re= 1.d2*esig/sig
      print 1,sig,esig
      print 2,re
      thc= 6.d0
      call ww(inpts,inrand,wmi,rsi,s0i,thc,tc,alszi,sig,esig)
      re= 1.d2*esig/sig
      print 1,sig,esig
      print 2,re
      thc= 5.d0
      call ww(inpts,inrand,wmi,rsi,s0i,thc,tc,alszi,sig,esig)
      re= 1.d2*esig/sig
      print 1,sig,esig
      print 2,re
*
    1 format(' sigma (nb)    =',e20.5,3x,'+-',3x,e20.5)
    2 format(/' Rel. error of ',f10.3,' %',///)
*
      stop
      end