system of non linear equations
2 views (last 30 days)
Show older comments
Hello
I am dealing with a system of two non-linear equations and I tried to solve it by 'fsolve', but the solver stopped and could not fine the answers.
Would you please help me?
I attached the files
0 Comments
Accepted Answer
Torsten
on 16 Oct 2023
%solving stsyem of nonlinear equation 3.31 and 3.32 Woods paper with fsolve
% eta_l= x1 , eta_u= x2
r=2;
v=20;
F=13/9;
fun= @(x)rooted(x,r,v,F);
x0 = [0 0];
x = fsolve(fun,x0)
function f = rooted (x,r,v,F)
%UNTITLED3 Summary of this function goes here
% Detailed explanation goes here
% F<Fd Fd= 2.6 I put F=1.5
f(1)= (r*v-1)/(v*(r-1)) * ((r*v*x(1)^2-x(2)^2)/(r*v-1))^(3/2) - ...
(v-1)/(v*(r-1)) * ((v*x(1)^2-x(2)^2)/(v-1))^(3/2) - 9 ;
f(2)= x(2)^3/(r*v) + (v-1)/(v*(r-1)) * ((v*x(1)^2-x(2)^2)/(v-1))^(3/2) - ...
(r*v-1)/(r*v*(r-1)) * ((r*v*x(1)^2-x(2)^2)/(r*v-1))^(3/2) - 9*F;
end
7 Comments
Torsten
on 24 Oct 2023
I don't know what the reason for non-convergence is. The equations seem to be correct - and concerning parameters r, v and F and initial values for x - I don't know if they make sense.
More Answers (0)
See Also
Categories
Find more on Systems of Nonlinear Equations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!