how can i write attached equation in MATLAB code
1 view (last 30 days)
Show older comments
Accepted Answer
Walter Roberson
on 17 Nov 2021
syms Delta_P_lo Delta_x phi_lo Delta_P_f x LB UB
Delta_P_f = Delta_P_lo * (1/Delta_x .* int(phi_lo.^2, x, LB, UB, 'hold', true))
The 'hold', true is there to prevent int() from immediately doing the integral. is independent of x so int() would treat it as being constant in x and would immediate do the integration, leaving you with something that did not include the integration sign.
If is dependent on x then the code would need to be changed a little to reflect that.
I had to introduce LB and UB as the bounds of integration because it happens that the 'hold' feature would be misinterpreted as being bounds if you were to try to do indefinite integration while 'hold' is on.
3 Comments
Walter Roberson
on 21 Nov 2021
Edited: Walter Roberson
on 21 Nov 2021
syms d g mu_1 mu_v rho_1 rho_v sigma chi
GAMMA = (rho_1/rho_v)^sym(1/2) / (mu_v / mu_1)^sym(1/8)
N_conf = 1/d * (sigma/ (g*(rho_1 - rho_v)))^sym(1/2)
phi_lo__2 = 1 + (sym(4.3) * GAMMA^2 - 1) * (N_conf * chi^sym(7/8) * (1-chi)^sym(7/8) + chi^sym(1/8))
syms Delta_chi Delta_P_f e_lo G L R
f_lo = (sym(0.79) * R * e_lo - sym(1.64))^(-2)
Delta_P_lo = 2*f_lo * G^2 * L / (rho_1 * d)
syms LB UB
Delta_P_f = Delta_P_lo * (1/Delta_chi .* int(phi_lo__2, chi, LB, UB, 'hold', true))
dPf = simplify(release(Delta_P_f))
Note that the Γ that shows up in the result is not the same as the Γ that you define at the bottom left of the equations. The one in the results is the Gamma function -- the extension of factorials.
More Answers (0)
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!