1 real function wid_a1_fitKKpi(qq)
15 real*8 qmax,qmin,qk_min(nq),w_qmin(nq),del_qq
19 real*8 aq,bq,qk_max,w_qmax
37 call inita1tabkkpi(qk_min ,w_qmin )
40 del_qq = (qmax-qmin)/float(nq-1)
46 if(qq.gt.qmin.and.qq.le.qmax)
then
53 aq = (w_qmax-w_qmin(kq))/(qk_max-qk_min(kq))
54 bq = (w_qmax*qk_min(kq) -w_qmin(kq)*qk_max)
55 $ /(qk_min(kq)-qk_max)
56 wid_a1_fitkkpi = aq*qq+bq
58 elseif(qq.ge.qmax)
then
59 wid_a1_fitkkpi = w_qmin(nq)
62 wid_a1_fitkkpi = w_qmin(1)