Equation involving Incomplete Gamma function
2 views (last 30 days)
Show older comments
How to solve this equation using matlab
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/568924/image.png)
where
and
and here m,
,
are constant.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/568929/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/568934/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/568939/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/568944/image.png)
3 Comments
Star Strider
on 1 Apr 2021
O.K. I missed seeing that m was constant. Code the expression and solve it in terms of only θ.
Answers (1)
Nipun
on 6 Jun 2024
Hi Sudhir,
I understand that you want to solve the given equation in MATLAB comprising of trigonometric and incomplete gamma functions.
I recommend using the "symbolic Math Toolbox" in MATLAB to create the equation and using a MATLAB solver to solve for "theta". MATLAB supports incomplete gamma functions using an in-built function, called "gammainc".
Here are the steps you may take to slve the given equation:
- Create symbolic variables
- Code the equation
- Use solver to solve for "theta"
I am attaching a code snippet with this answer for your reference:
syms theta m tmin tmax
a = 2*sin(theta).^2
b = (cos(theta)-sin(theta)).^2/2
% Equation (split into two parts)
part1 = (4/(cot(theta)-1).^2)^((m+2)/(m+3))
arg1 = (m+1)/(2*m+6);
arg2 = b*tmin/2;
arg3 = a*tmin/2;
arg4 = b*tmax/2;
arg5 = a*tmax/2;
part2_numerator = gammainc(arg1, arg2) - gammainc(arg1, arg4);
part2_denominator = gammainc(arg1, arg3) - gammainc(arg1, arg5);
part2 = part2_numerator/part2_denominator;
% creating equation
eqn = part1*part2
% Solving
ans = solve(eqn, theta);
For more information on the in-built functions, refer to the following MathWorks documentation:
- Symbolic Toolbox: https://www.mathworks.com/help/symbolic/syms.html#d126e328782
- Incomplete Gamma Function: https://www.mathworks.com/help/matlab/ref/gammainc.html
- Equation Solver: https://www.mathworks.com/help/symbolic/sym.solve.html
Hope this helps.
Regards,
Nipun
0 Comments
See Also
Categories
Find more on Nonlinear Optimization 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!