http://www.mathworks.com/matlabcentral/newsreader/view_thread/242583
MATLAB Central Newsreader  integration explicit solution could not be found
Feed for thread: integration explicit solution could not be found
enus
©19942014 by MathWorks, Inc.
webmaster@mathworks.com
MATLAB Central Newsreader
http://blogs.law.harvard.edu/tech/rss
60
MathWorks
http://www.mathworks.com/images/membrane_icon.gif

Sat, 17 Jan 2009 00:55:03 +0000
integration explicit solution could not be found
http://www.mathworks.com/matlabcentral/newsreader/view_thread/242583#622208
Dhritiman
i had an earlier post about this with a typo. My problem is that matlab cannot solve the following equation where some function of z multiplied by the pdf of gamma dist is integrated over the range [.3 .9]. I need the answer as a function of w(wage). Infact I did get a scalar answer after i used double(int(f,z,.3,.9)) but this is only working when i substitute w for a number. Is there a way out of this...i.e. I need the soultion as a function of w (to be used later inside a fzero)<br>
<br>
<br>
syms z w<br>
f=138.2528*84.626*z*((10.99*z)^1.32)*((1/w)^3.68)*((1/((2)*.06^3))*(z^(31))*(2.718^(z/.06)))<br>
intfs=int(f,z,.3,.9)<br>
<br>
thanks for your help!!

Sat, 17 Jan 2009 17:45:04 +0000
Re: integration explicit solution could not be found
http://www.mathworks.com/matlabcentral/newsreader/view_thread/242583#622276
Roger Stafford
"Dhritiman " <dhritimanbhattacharya@uiowa.edu> wrote in message <gkra97$24k$1@fred.mathworks.com>...<br>
> i had an earlier post about this with a typo. My problem is that matlab cannot solve the following equation where some function of z multiplied by the pdf of gamma dist is integrated over the range [.3 .9]. I need the answer as a function of w(wage). Infact I did get a scalar answer after i used double(int(f,z,.3,.9)) but this is only working when i substitute w for a number. Is there a way out of this...i.e. I need the soultion as a function of w (to be used later inside a fzero)<br>
> <br>
> syms z w<br>
> f=138.2528*84.626*z*((10.99*z)^1.32)*((1/w)^3.68)*((1/((2)*.06^3))*(z^(31))*(2.718^(z/.06)))<br>
> intfs=int(f,z,.3,.9)<br>
<br>
Dhritiman, I call to your attention again the fact that the factor involving w in your expression need not be included in the integrand here but can be factored outside. What you have can be expressed as<br>
<br>
a/w^3.68 * int(z^3*(1b*z)^c*exp(z/d,z,e,f)<br>
<br>
where a = 138.2528*84.626/(2*.06^3), b = 0.99, c = 1.32, d = .06, and where e and f are the lower and upper limits of integration. The value of the integral part depends only on constants b, c, d and the limits of integration, not on a or w. You can multiply these latter in after the integration is performed.<br>
<br>
As for the integral, the introduction of the power c = 1.32 apparently renders this integral incapable of an explicit solution, but this is a mathematical difficulty and not one that can be blamed on matlab. Consequently it will be necessary for you to do numerical integration to evaluate it for each different set of constants b, c, d and limits of integration.<br>
<br>
About the only saving I can see in this will come from cases where the integration range is taken over various (e,f) intervals that overlap. The overlap portions need be integrated only once (provided b, c, and d remain fixed.)<br>
<br>
Roger Stafford

Sun, 18 Jan 2009 04:27:01 +0000
Re: integration explicit solution could not be found
http://www.mathworks.com/matlabcentral/newsreader/view_thread/242583#622318
Dhritiman
thanks a tonne Roger.....you nailed it!!!!