Path: news.mathworks.com!not-for-mail From: <HIDDEN> Newsgroups: comp.soft-sys.matlab Subject: Re: Can I integrate the following function in MATLAB?? Date: Tue, 31 May 2011 04:08:02 +0000 (UTC) Organization: The MathWorks, Inc. Lines: 12 Message-ID: <is1pj2$15a$1@newscl01ah.mathworks.com> References: <is17bu$i8i$1@newscl01ah.mathworks.com> <is1910$lsl$1@newscl01ah.mathworks.com> Reply-To: <HIDDEN> NNTP-Posting-Host: www-06-blr.mathworks.com Content-Type: text/plain; charset=UTF-8; format=flowed Content-Transfer-Encoding: 8bit X-Trace: newscl01ah.mathworks.com 1306814882 1194 172.30.248.38 (31 May 2011 04:08:02 GMT) X-Complaints-To: news@mathworks.com NNTP-Posting-Date: Tue, 31 May 2011 04:08:02 +0000 (UTC) X-Newsreader: MATLAB Central Newsreader 1648643 Xref: news.mathworks.com comp.soft-sys.matlab:729413 "Roger Stafford" wrote in message <is1910$lsl$1@newscl01ah.mathworks.com>... > "kumar vishwajeet" wrote in message <is17bu$i8i$1@newscl01ah.mathworks.com>... > > I want to numerically integrate the following function. Is it possible to do it in MATLAB?? I tried quad2d, dblquad. But there is no scope of symbolic constants in these commands. > > f(x,y) = log(1-x)*exp(-a*log(x) - b*log(1-x) - c*log(y) - d*log(1 - y) - e*x*y) > > > > where a,b,c,d,e are symbolic constants and limit of integration is (0,1) for x and (0,1) for y. > - - - - - - - - - > I can't reconcile the two aspects of your request. On the one hand you seem to be asking for a general formula for the above double integral as a function of a, b, c, d, and e, and yet on the other hand you want to do the integration numerically. It is conceivable that the symbolic toolbox could come up with the first part, (though I have serious doubts,) but I don't really know what you mean by doing the integration numerically in such a way that it incorporates the five undefined symbols. Perhaps you can clarify your question. > > Roger Stafford I want to integrate the function in such a way that the result is still a function of a,b,c,d,e. I tried to do it using symbolic toolbox, but did not get any explicit solution. That's why I want to do it numerically.