matlab特殊函数积分

matlab特殊函数积分
如图所示,int((log(eta*x)+expint(lambda*x))/x+log(lambda.^2*eta.^3*x+expint(-lambda*x))/(x*exp(lambda.*x)),x,0,tau)这是我写的公式,但是expint函数再积分的话,就会报错,请问有没有好的方法.有错的话请改正一下.实在不行的话,数值积分也可以,但最好是直接积分
psii=((1+0.75*gamma.*epsilon)*0.5)*(log(tau)+0.80908+expint(lambda.*tau)+2*S/(1+0.75*gamma.*epsilon))-0.125*epsilon*(pi.^2/3-9*gamma.^2/4+log(tau/100)*log(100*eta.^2*tau))
      -0.25*epsilon.*(pi^2/6-9*gamma.^2/8+int((log(eta*x)+expint(lambda*x))/x+(log(lambda.^2*eta.^3*x/256)-expint(-lambda*x))/(x*exp(lambda.*x)),x,0,tau))