Re: Solving a weakly singular integral equation

*To*: mathgroup at smc.vnet.net*Subject*: [mg53827] Re: Solving a weakly singular integral equation*From*: Paul Abbott <paul at physics.uwa.edu.au>*Date*: Sat, 29 Jan 2005 06:02:35 -0500 (EST)*Organization*: The University of Western Australia*References*: <ctahsj$b14$1@smc.vnet.net>*Sender*: owner-wri-mathgroup at wolfram.com

In article <ctahsj$b14$1 at smc.vnet.net>, Zaeem Burq <Z.Burq at ms.unimelb.edu.au> wrote: > I am trying to solve a non-linear Volterra type-2 integral equation by > using successive approximations method. As far as I can see, your integral equation is not nonlinear (it is certainly linear in f). > I am mainly interested in the behaviour of the solution near zero. I believe that the solution is singular there so successive (numerical) approximation is unlikely to help. > The unknown function is f[t]. Define > > p[x_]:= Exp[- 0.5 * x^2]/(2 Pi) (* Gaussian probability density function*) This is not a pdf (it does not integrate to 1). It would be better to define this as p[x_] = Exp[-x^2/2]/Sqrt[2 Pi] where exact quantities replace approximate ones. > c[t]:= Sqrt[t] + 2t/9 As you write later, this is c[t_] = Sqrt[t] + 2t/9; Note that it is usually more efficient to use = instead of := in function definitions when the right-hand side is explicit. > G[s_,t_]:= (c[t] - c[s]) p[(c[t]-c[s])/Sqrt[t-s]] / (t-s)^1.5 (* The > kernel *) Write this as G[s_,t_] = (c[t] - c[s]) p[(c[t]-c[s])/Sqrt[t-s]] / (t-s)^(3/2); Note, however, that this expression is real only for t >= s. > The integral equation is: > > f[t] = c[t]p[c[t]/Sqrt[t]] / t^1.5 - \int_{0}^{t} G[s,t] f[s] ds In Mathematica notation, you are trying to solve f[t]==c[t] p[c[t]/Sqrt[t]]/t^(3/2) - Integrate[G[s, t] f[s], {s, 0, t}] > As you can see, the perturbation function has an apparant singularity at > 0, and the kernel has singularities along the diagonal (s=t). All the > singularities are appropriately killed by the function p, but Mathematica > seems to have trouble with them. Series expansion of the inhomogenous term about t = 0, i.e., (c[t] p[c[t]/Sqrt[t]])/t^(3/2) + O[t] shows that it goes like 1/(Sqrt[2 E Pi] t) as t -> 0. Hence the singularity is not eliminated by the function p. Morevover, I think it is possible to show that there is _no_ solution to this integral equation: first note that, using the change of variables s -> x t, the integral part reads t Integrate[G[x t, t] f[x t], {x, 0, 1}] Now, if f[t] ~ t^q as t -> 0, then series expansion of the integrand with respect to t leads to computing the following integral: Integrate[(E^((Sqrt[x] - 1)/(2(Sqrt[x] + 1)))* (1 - Sqrt[x]) x^q)/(1 - x)^(3/2), {x, 0, 1}] which diverges if q <= -1. But, to agree with the expansion of the inhomogenous term, we require a term with q = -1. On the other hand, a slight modification of the inhomogenous term, say changing t^(3/2) -> t, leads to f[t] == c[t] p[c[t]/Sqrt[t]]/t - Integrate[G[s, t] f[s], {s, 0, t}] an equation that can be easily solved by appropriate (using fractional powers) series expansion of f. For some relevant references on integral equations see http://mathworld.wolfram.com/IntegralEquation.html Note that Numerical Recipes is online at http://www.library.cornell.edu/nr/bookcpdf.html The following sections are relevant: 18.2 Volterra Equations 794 18.3 Integral Equations with Singular Kernels 797 The method that you were trying to apply to this problem looks like that outlined in the paper by Stan Richardson on successive approximation solutions to integral equations, which appeared in The Mathematica Journal 9(2): 460-482. Cheers, Paul -- Paul Abbott Phone: +61 8 6488 2734 School of Physics, M013 Fax: +61 8 6488 1014 The University of Western Australia (CRICOS Provider No 00126G) 35 Stirling Highway Crawley WA 6009 mailto:paul at physics.uwa.edu.au AUSTRALIA http://physics.uwa.edu.au/~paul