def devi(yy, eta): deveta = yy*eta - scipy.exp(eta) devy = yy*scipy.log(yy) - yy devy[yy == 0] = 0 result = 2*(devy - deveta) return(result)