Solve numerical equation with Y at both sides
7 views (last 30 days)
Show older comments
Hi all,
I have a general equation I would like to solve. However, it contains the dependent variable Y on both sides.
The eqaution is: Y/Y_0 = B * (1 + x*Y/D)^0.5 * (1-x)^alpha
Where x is the independent variable,
and Y_0, B, D and alpha are constants. What would be the most appropriate function in Matlab for this particular problem?
Thanks,
Ron
0 Comments
Answers (3)
KSSV
on 23 Aug 2020
You can use symbolic package. Something like this:
syms Y_0 Y B x Y D alpha
eqn = Y/Y_0 - B * (1 + x*Y/D)^0.5 * (1-x)^alpha==0 ;
s = solve(eqn,Y)
Star Strider
on 23 Aug 2020
A numeric approach:
Y_0 = 3; % Define Constants
B = 5;
D = 7;
alpha = 11;
x = linspace(0, 0.9, 10);
Yfcn = @(Y,x) B * sqrt(1 + x*Y/D) .* (1-x).^alpha - Y/Y_0; % Use Element-Wise Operations
for k = 1:numel(x)
Y(k) = fzero(@(Y)Yfcn(Y,x(k)), 0.1);
end
figure
plot(x, Y)
grid
2 Comments
Alan Stevens
on 23 Aug 2020
Edited: Alan Stevens
on 23 Aug 2020
Your equation can also be expressed as a quadratic in Y
which could be solved using roots (for specified values of x).
You would need to check that the solutions were consistent with the original equation.
See Also
Categories
Find more on Calculus 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!