```Path: news.mathworks.com!not-for-mail
From: <HIDDEN>
Newsgroups: comp.soft-sys.matlab
Subject: Re: integral equation with singularities and double integral
Date: Mon, 31 Oct 2011 13:31:28 +0000 (UTC)
Organization: Saab Underwater Systems AB
Lines: 35
Message-ID: <j8m7vg\$1m\$1@newscl01ah.mathworks.com>
References: <j8c3fp\$8kv\$1@newscl01ah.mathworks.com>
NNTP-Posting-Host: www-02-blr.mathworks.com
Content-Type: text/plain; charset=UTF-8; format=flowed
Content-Transfer-Encoding: 8bit
X-Trace: newscl01ah.mathworks.com 1320067888 54 172.30.248.47 (31 Oct 2011 13:31:28 GMT)
X-Complaints-To: news@mathworks.com
NNTP-Posting-Date: Mon, 31 Oct 2011 13:31:28 +0000 (UTC)
Xref: news.mathworks.com comp.soft-sys.matlab:747704

"Sulaiman Rabbaa" wrote in message <j8c3fp\$8kv\$1@newscl01ah.mathworks.com>...
> Hello
> I need help to solve the following integral equation:
> double integral (with respect to x and y) of 137.03.*y.^2./((0.238.*exp(0.067.*y.^2)+1).*(w-5.26.*x.*y-2.63).*(w-5.26.*x.*y+2.63))=1+8478./(10828-w.^2-1.13.*j.*w)
> xmin=-1, xmax.=1, ymin=0, ymax=inf
> I want to find w which is complex number.
> I tried the following code:
> V=fsolve(f,160)
> There are some singularities. How can I solve the problem

Hi,
The Newton-Raphson solution below gives w = (78.1093 - 0.2275i). Another solution is w = (-78.1093 - 0.2275i).

f = @(x,y,w) 137.03.*y.^2./(0.238.*exp(0.067.*y.^2)+1) ...
./(w-5.26.*x.*y-2.63)./(w-5.26.*x.*y+2.63);

df = @(x,y,w) 137.03.*y.^2./(0.238.*exp(0.067.*y.^2)+1) ...
./(w-5.26.*x.*y-2.63).^2./(w-5.26.*x.*y+2.63).^2 ...
.*2.*(5.26.*x.*y-w);