Problems using integral function in Matlab

  • MATLAB
  • Thread starter cacosomoza
  • Start date
  • #1
11
0
Hi there!

i have this function handle:

integrando= @(x,landa,t) 1-exp(-((1/(landa*t))*((exp(-landa*x)/x)+((x-x0).*(1+landa*x0)-x0)*(exp(-landa*x0)/(x0^2)))));

which as you can see using

fplot(@(x) integrando(x,3,0.00063),[0 1.1225])

its value is constant (actually 1, and then starts decreasing exponentially until drops to zero at 1.1225). Landa and t are parameters that change the curvature of this drop. Since this is a positive and finite function, I'm waiting a positive integral. However, integrating this same function on this same interval, I get negative figures!

>> q1=integral(@(x) integrando(x,3,0.00063),0,1.1225)

q1 =

-20.3981



This is getting me nuts! I would very much appreciate any help. Thanks in advance!
 

Answers and Replies

  • #2
kreil
Insights Author
Gold Member
668
68
What value are you using for x0?

I used 2 just to see if the code works, and I got:

integrando= @(x,landa,t) 1-exp(-((1/(landa*t))*((exp(-landa*x)/x)+((x-2).*(1+landa*2)-2)*(exp(-landa*2)/(2^2)))));

fplot(@(x) integrando(x,3,0.00063),[0 1.1225])

q1=integral(@(x) integrando(x,3,0.00063),0,1.1225)

q1 =

1.1225
 
  • #3
11
0
hi kreil,

I'm using x0=1.1225 ,which is the same value as the limit of integration. If you try again using that value on 'integrando', you'll see the result of the integral is a negative number which doesn't make any sense to me
 

Related Threads on Problems using integral function in Matlab

Replies
1
Views
4K
  • Last Post
Replies
2
Views
46K
  • Last Post
Replies
1
Views
5K
Replies
1
Views
9K
  • Last Post
Replies
2
Views
9K
  • Last Post
Replies
2
Views
2K
Replies
3
Views
2K
Replies
0
Views
2K
Replies
1
Views
6K
Replies
3
Views
9K
Top