real*16 xx(10),xg,wg,xgl,wgl,a,b,x,Kr,Ki,canr,cani
c      data xx/1.q0,5.q0,10.q0,20.q0,50.q0,100.q0/
      data xx/5.q-1,1.q-1,5.q-2,1.q-2,1.q-3,1.q-4,
     *1.q-5,1.q-6,1.q-7,1.q-8/
      common/gauss/xg(55),wg(55),xgl(60),wgl(60)
      write(*,1)
    1 format(/)
      read(*,2) (xg(k),wg(k),k=1,55)
    2 format(1x,2d38.30)
      read(*,2) (xgl(k-55),wgl(k-55),k=56,115)
c      do 110 k=1,55
c        write(*,101) xg(k),wg(k)
c  101   format(1x,2d24.16)
c  110 continue
c      write(*,1)
c      do 120 k=1,60
c        write(*,101) xgl(k),wgl(k)
c  120 continue
      a=0.q0
      write(*,3)
    3 format(3x,'beta alpha',5x,'x',17x,'Kr/Ki',16x,'canr/cani'/) 
      do 20 ib=0,10
        b=qreal(ib)*1.q0
        do 10 ix=7,10
          x=xx(ix)
          call qmacdonald(x,a,b,Kr,Ki,canr,cani)
          write(*,4) snglq(b),snglq(a),snglq(x),Kr,snglq(canr)
    4     format(1x,2f6.2,f8.2,d34.26,d12.4)
          write(*,5) Ki,cani
    5     format(21x,d34.26,d12.4)
   10   continue
        write(*,1)
   20 continue 
      stop
      end