MATLAB Answers

Mupad interface- solving equations

6 views (last 30 days)
Mirsad on 14 Jan 2017
Commented: Mirsad on 16 Jan 2017
Hello everyone,
I have an equation that written in Mupad;
F := 38000 * ln((K*n)/(K*n - 250)) - 9500000.0/(K*n - 250) + 7000 = 0;
Here, I want to calculate the K*n value. My purpose is to get this parameter from the equation. When I use solve command such that F:= solve(F,K*n), it gives me only one solution of the equation. How can I see all possible solutions?. I had only one negative solution but there must be some positive solution also but I haven't seen yet.
Note: I also try to limit the parameter range, (K*n = 1..1000) but still doesn't work. Best regards, Mirsad


Sign in to comment.

Answers (1)

Walter Roberson
Walter Roberson on 14 Jan 2017
Delete the ".0" from the constant in order to have the two real-valued solutions returned.
There are an infinite number of solutions corresponding to different branches of the LambertW function. I do not know how to get MuPAD to return the generalized form. The generalized form is
250*LambertW(Z, -exp(-45/38))/(LambertW(Z, -exp(-45/38))+1)
for integer Z, the different roots of LambertW. Z = -1 and Z = 0 are the real-valued solutions.


Mirsad on 14 Jan 2017
Thank you for your reply. The problem is solved in this manner:
S := solve(F1=0, K*n)
When I enter the command, I obtain two solution ({-259.2690041, 589.8981938}). I try to achieve 589.89. My purpose is to prove What I get from the mupad is same with textbook solution at the end. According to your answer, I still do not understand how to get these all possible solutions. Because you say that there are an infinite number of solutions corresponding to different branchses of the lambertw funct. I have to look deeper lamberW function.
Walter Roberson
Walter Roberson on 16 Jan 2017
F := 38000 * ln((K*n)/(K*n - 250)) - 9500000/(K*n - 250) + 7000 = 0;
select( solve(F, K*n), x->float(x)>0);
The description of the LambertW function is not bad; see
Mirsad on 16 Jan 2017
Thank you, the code you write is the better way to get what I need. Yes, lammbert function is more relevant now.

Sign in to comment.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!