      function die(x)

      common/pawpar/p(4)

      if(x.gt.(p(4)+1)) then
         die=p(1)+p(3)*(x-p(4))**p(2)
      else
         die=0.
      endif

      end
