Problem Calculating Nonlinear Indefinite Integral

1 visualización (últimos 30 días)
Hello,
I am trying to calculate a nonlinear integral t^5.2 *exp((-x^6.2-x*(0.2*4.5^6.2))/4.5^6.2) in order to solve an equation, where the unknown is the lower limit of integrations.Although i use "int" in order to calculate the integral, Matlab returns the equation.Is my equation too complicated or there is another way to get a solution ?I am quite new in Matlab, so any help would be appreciated.
I would like to thank you in advance for the time dealing with my query.
  3 comentarios
Dimitrios Kampitsis
Dimitrios Kampitsis el 1 de Dic. de 2019
My goal is to solve the following equation where the unknown is the limit of integration (k), where a and b are constants
1.jpg
I tried only with int in order to solve the integral and afterwards, if i get a solution from my integration, i am planning to try and solve the equation. However i am not sure if this is the right way to do it.
syms x
g = x^5 *exp((-x^5-x*(4^6))/4^6) ;
int(g,x,k,k+2)
I havent tried yet trapz. Moreover i didnt try integral because i am looking for a symbolic solution.
Thanks again.
darova
darova el 1 de Dic. de 2019
You should use fsolve and integral

Iniciar sesión para comentar.

Respuesta aceptada

David Wilson
David Wilson el 2 de Dic. de 2019
Since you neglected to tell us values for a and b, I'll just set them to a=2 and b=1. You have also two definitions for g(x), so I've used the one in your code. Like the previous poster suggested, I'll embed integral inside fsolve:
g = @(x) x.^5.*exp((-x.^5-x*(4^6))/4^6) ;
Q = @(k) integral(g,k,k+2);
a = 2; b = 1; % Set these to whatever
k0 = 1; % Need a start guess here.
k = fsolve(@(k) Q(k) - a/b, k0)
  1 comentario
Dimitrios Kampitsis
Dimitrios Kampitsis el 4 de Dic. de 2019
Thank you for your help, it works and i get a solution. I have tried it also with a more complex equation and it works fine. I really appreciate your help.
However for a specific value of my constants fsolve returns a value (Equation solved) that it is not a root of my equation. I confirmed that using another software "Mathcad".
Could this be a problem of function tolerance ?

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Numerical Integration and Differentiation en Help Center y File Exchange.

Productos

Community Treasure Hunt

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

Start Hunting!

Translated by