Path: news.mathworks.com!not-for-mail From: <HIDDEN> Newsgroups: comp.soft-sys.matlab Subject: Re: Integration Date: Wed, 11 Aug 2010 02:01:07 +0000 (UTC) Organization: The MathWorks, Inc. Lines: 9 Message-ID: <i3t093$mk5$1@fred.mathworks.com> References: <i3rs7s$a7p$1@fred.mathworks.com> <i3s1pr$da7$1@fred.mathworks.com> <i3spl3$t1v$1@fred.mathworks.com> Reply-To: <HIDDEN> NNTP-Posting-Host: webapp-02-blr.mathworks.com Content-Type: text/plain; charset=UTF-8; format=flowed Content-Transfer-Encoding: 8bit X-Trace: fred.mathworks.com 1281492067 23173 172.30.248.37 (11 Aug 2010 02:01:07 GMT) X-Complaints-To: news@mathworks.com NNTP-Posting-Date: Wed, 11 Aug 2010 02:01:07 +0000 (UTC) X-Newsreader: MATLAB Central Newsreader 1187260 Xref: news.mathworks.com comp.soft-sys.matlab:660871 "Anna Kaladze" <anna.kaladze@gmail.com> wrote in message <i3spl3$t1v$1@fred.mathworks.com>... > Yes, thanks a lot for the answer and sorry, F(u) and f(u) are the SAME functions -- just a typo. > Is there any code I can to write to solve the problem? I mean technically it is possible to solve the problem in Excel (one column for the inner integral where t argument will take the value from 0 to whatever), and then sum-up the values in that column using the trapezoidal rule. A smaller step size would give a reasoanble degree of approximation). But is there a way to do something like that in MATLAB? Thanks a lot. - - - - - - - - - - - - It's the statement "I have a non-integrable function, f(u)" that you made in the first post that is the stumbling block here. In mathematics and in matlab circles too, when you say a function is non-integrable, it is because that function is sufficiently ill-behaved over the desired integration range that it is impossible to obtain an integral for it. Perhaps it ascends to infinity in the wrong way, the range is infinite and it doesn't get small fast enough, or it is seriously discontinuous. An example is the integral of 1/x^2*sin(1/x) from 0 to 2/pi which is not well-behaved as it approaches x = 0. The integrated value keeps oscillating endlessly back and forth more and more rapidly from -1 to +1 as the lower limit approaches zero and consequently is non-integrable over that full range. I am guessing since you are still talking about trying to get your function's integral that this isn't what you meant by "non-integrable". If so, you should take the advice you were given more seriously. The double quadrature routine 'quad2d' allows for varying limits of integration in its inner integral, which is what it sounds like you are faced with when you say, "The inner integral (where integrand is F(u)) has a low limit 0, but the upper limit is t (in principle, t takes the value from 0 to 1)." I suggest you look into it. Roger Stafford