function tblfit(plab,a,b,en,c,d)
      tblfit=a+b*plab**en+c*(log(plab))**2+           
     a                              d*log(plab)                         
      return
      end