implicit real*8(a-h,o-z)
*
      inpts= 6
      inrand= 6
      wmi= 80.39d0
      rsi= 186.d0
      tc= tc/180.d0*3.141592653589793238462643d0
      alszi= 0.1194d0
      s0i= 0.01d0
      tc= 0.5d0
      tc= tc/180.d0*3.141592653589793238462643d0
*
      scalei= 100.d0
      call ww(inpts,inrand,scalei,wmi,rsi,s0i,tc,alszi,sig,esig)
      re= 1.d2*esig/sig
      print 1,sig,esig
      print 2,re
      scalei= 10.d0
      call ww(inpts,inrand,scalei,wmi,rsi,s0i,tc,alszi,sig,esig)
      re= 1.d2*esig/sig
      print 1,sig,esig
      print 2,re
      scalei= 1.d0
      call ww(inpts,inrand,scalei,wmi,rsi,s0i,tc,alszi,sig,esig)
      re= 1.d2*esig/sig
      print 1,sig,esig
      print 2,re
      scalei= 0.1d0
      call ww(inpts,inrand,scalei,wmi,rsi,s0i,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