这其实就是一个一阶线性微分方程。y'+y/(xlnx)=(lnx+1)/lnx，用常数变异法或直接公式解。这里简单叙述一下常数变异法：首先解y'+y/(xlnx)=0，得y=c/lnx，c为任意常数再设原方程的解为y=c(x)/lnx，带入原方程整理求积分可解得c(x)=xlnx+C所以得通解y=(xlnx+C)/lnx=x+ C/lnx
用matlab,很快
