dsolve('x*y-Dy*ln(y)=0','x')
ans =
exp(2^(1/2)*(x^2/2 + C7)^(1/2))
1/exp(2^(1/2)*(x^2/2 + C7)^(1/2))
这是使用Matlab计算结果