the function 'solve' doesn't give all the solutions.
1 次查看(过去 30 天)
显示 更早的评论
% this does not give all the solutions
syms u;
solve(1-0.99-(u+1)*exp(-u));
double(ans)
see the short 3-line code there, try it out.
The solutions should be u1=-0.996308, and u2=6.63835
But this code only gives me u1, and mentions nothing about u2.
Can any one tells me why or which function should I use to get all the answers.
Thanks!
2 个评论
Roger Stafford
2013-10-9
My very ancient version of 'solve' gives both your solutions. I don't know why your version should give only one of them.
Walter Roberson
2013-10-9
I have noticed the MuPAD's solve() seems to have difficulty in finding solutions involving LambertW.
采纳的回答
Friedrich
2013-10-9
编辑:Friedrich
2013-10-9
Hi,
seems like you are interested in real solutions only. So use:
>> solve(1-0.99-(u+1)*exp(-u),'Real',true)
ans =
- lambertw(0, -exp(-1)/100) - 1
- lambertw(-1, -exp(-1)/100) - 1
If you dont set Real to true MATLAB solves it Complex which would give an infinite amount of solutions. When you use MuPAD Notepad you will obtain:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/172195/image.png)
So only k=0 and k =-1 give Real solutions.
UPDATE: Since you want real solutions only its might be better to use a numeric solver:
>> feval(symengine,'numeric::solve','1-0.99-(u+1)*exp(-u)','AllRealRoots')
ans =
[ [u == -0.99630759688740599414645321772374], [u == 6.6383520679938122693712394050244]]
更多回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!