intsolve - Maple Help

intsolve

solve a linear integral equation

 Calling Sequence intsolve( Inteqn, funcn ) intsolve( Inteqn, funcn, method = ..., order = ... )

Parameters

 Inteqn - linear integral equation in funcn funcn - unknown function to be solved for method = ... - optional, the right-hand side can be any of: collocation, differentialequation, eigenfunction, Neumann, Laplace order = ... - optional, the right-hand side can be a positive integer, specifying the order of the series approximation to use with the Neumann method, or the degree of polynomial to use with the collocation method

Description

 • The intsolve command is able to find exact closed-form solutions to various types of integral equations as explained below. The solution is returned as an explicit equation with the unknown function in the left-hand side and the solution in the right-hand side. In some cases arbitrary constants are returned in the solution; they are represented as $\mathrm{_C1},\mathrm{_C2},\dots ,\mathrm{_Cn}$, as in the output of dsolve. When the solution contains integrals, they are represented with the inert Int.
 • By default, intsolve converts, using the differentialequation method, the linear integral equation into an equivalent ODE problem with initial values (ODE-IVP), then attempts solving this problem; if it fails, it uses an eigenfunction approach. The method Laplace is optional and tackles the problem using a Laplace transform. The collocation and Neumann methods are also optional, and compute approximations to the solution.
 • intsolve can recognize Fredholm and Volterra integral equations of the following kinds:
 Fredholm equation of the first kind (F1):

${\int }_{a}^{b}k\left(x,y\right)f\left(y\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}ⅆy=g\left(x\right)$

 Fredholm equation of the second kind (F2):

$f\left(x\right)+{\int }_{a}^{b}k\left(x,y\right)f\left(y\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}ⅆy=g\left(x\right)$

 Fredholm equation of the third (homogeneous) kind (F3):

$f\left(x\right)+{\int }_{a}^{b}k\left(x,y\right)f\left(y\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}ⅆy=0$

 Volterra equation of the first kind (V1):

${\int }_{a}^{x}k\left(x,y\right)f\left(y\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}ⅆy=g\left(x\right)$

 Volterra equation of the second kind (V2):

$f\left(x\right)+{\int }_{a}^{x}k\left(x,y\right)f\left(y\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}ⅆy=g\left(x\right)$

 Volterra equation of the third (homogeneous) kind (V3):

$f\left(x\right)+{\int }_{a}^{x}k\left(x,y\right)f\left(y\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}ⅆy=0$

 • The Laplace method is applicable to equations of V1, V2 and V3 type. One advantage to using this option is that the integral equation may contain the Dirac or Heaviside functions. This method is not suitable for finding a solution when the equation does not possess a convolution type kernel.
 • The collocation method for approximating solutions is applicable to Fredholm equations. This technique first partitions the interval $\left[a,b\right]$ uniformly into $n+1$ points, where positive integer $n$ is 10 by default or set by passing order=n. Then, the coefficients of an undetermined polynomial of degree $n$ are determined so that the polynomial satisfies the integral equation at the collocation points. The solve command is used to solve the collocation equations, so if a numeric approximation is desired, you can include floats inside the integral equation.
 • The Neumann method is applicable to equations of V2, F2, and some of V1 type. The default order of the Neumann series in the solution is Order - this value can be overridden by passing the optional argument order = n, where n is a positive integer.
 • When method is not specified, intsolve will check if the integral within the integral equation is of the form of a Fourier-Sine, Fourier-Cosine, or Laplace transform. That is, it will check if the integral equation can be written

${{\int }}_{0}^{\mathrm{\infty }}f\left(y\right)a\left(y\right)b\left(xy\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}y=c\left(x\right)$

 where $b\left(z\right)$ is one of $\mathrm{sin}\left(z\right)$, $\mathrm{cos}\left(z\right)$, or ${ⅇ}^{-z}$. In this case, intsolve will apply the appropriate inverse integral transform.
 • infolevel[intsolve] := 1,2,3 or 5 gives descriptive information on the progress of the computation. 5 may result in excessive printing.

Examples

Example 1

 • Consider the following example:
 > $\mathrm{eq1}≔p\left(x\right)-\frac{1}{2}\left(\mathrm{Int}\left(xyp\left(y\right),y=0..1\right)\right)=\frac{5}{6}x$
 ${\mathrm{eq1}}{≔}{p}{}\left({x}\right){-}\frac{\left({{\int }}_{{0}}^{{1}}{x}{}{y}{}{p}{}\left({y}\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{y}\right)}{{2}}{=}\frac{{5}{}{x}}{{6}}$ (1)
 > $\mathrm{intsolve}\left(\mathrm{eq1},p\left(x\right)\right)$
 ${p}{}\left({x}\right){=}{x}$ (2)
 If we instead use the 'Neumann' method, an approximation to this result is returned:
 > $\mathrm{intsolve}\left(\mathrm{eq1},p\left(x\right),\mathrm{method}=\mathrm{Neumann}\right)$
 ${p}{}\left({x}\right){=}\frac{{279935}{}{x}}{{279936}}$ (3)
 > $\mathrm{intsolve}\left(\mathrm{eq1},p\left(x\right),\mathrm{method}=\mathrm{Neumann},\mathrm{order}=20\right)$
 ${p}{}\left({x}\right){=}\frac{{21936950640377855}{}{x}}{{21936950640377856}}$ (4)

Example 2

 Computing the solution and verifying the result:
 > $\mathrm{eq2}≔f\left(x\right)=x+1+\mathrm{Int}\left(\left(1+2\left(x-y\right)\right)f\left(y\right),y=0..x\right)$
 ${\mathrm{eq2}}{≔}{f}{}\left({x}\right){=}{x}{+}{1}{+}{{\int }}_{{0}}^{{x}}\left({1}{+}{2}{}{x}{-}{2}{}{y}\right){}{f}{}\left({y}\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{y}$ (5)
 > $\mathrm{intsolve}\left(\mathrm{eq2},f\left(x\right)\right)$
 ${f}{}\left({x}\right){=}{{ⅇ}}^{{2}{}{x}}$ (6)
 To verify this solution, transform it first into a procedure:
 > $f=\mathrm{unapply}\left(\mathrm{rhs}\left(\right),x\right)$
 ${f}{=}\left({x}{↦}{{ⅇ}}^{{2}{\cdot }{x}}\right)$ (7)
 > $\mathrm{eval}\left(\mathrm{eq2},\right)$
 ${{ⅇ}}^{{2}{}{x}}{=}{x}{+}{1}{+}{{\int }}_{{0}}^{{x}}\left({1}{+}{2}{}{x}{-}{2}{}{y}\right){}{{ⅇ}}^{{2}{}{y}}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{y}$ (8)
 > $\mathrm{value}\left(\right)$
 ${{ⅇ}}^{{2}{}{x}}{=}{{ⅇ}}^{{2}{}{x}}$ (9)

Example 3

 For some type of problems the solution obtained using the Laplace method contains Dirac functions:
 > $\mathrm{eq3}≔\mathrm{Int}\left(\mathrm{exp}\left(a\left(x-y\right)\right)f\left(y\right),y=0..x\right)=1$
 ${\mathrm{eq3}}{≔}{{\int }}_{{0}}^{{x}}{{ⅇ}}^{{a}{}\left({x}{-}{y}\right)}{}{f}{}\left({y}\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{y}{=}{1}$ (10)
 > $\mathrm{intsolve}\left(\mathrm{eq3},f\left(x\right),\mathrm{method}=\mathrm{Laplace}\right)$
 ${f}{}\left({x}\right){=}{\mathrm{Dirac}}{}\left({x}\right){-}{a}$ (11)

Example 4

 Another type of integral equation automatically solved using the eigenfunction approach:
 > $\mathrm{eq4}≔f\left(x\right)+b\left(\mathrm{Int}\left(\left(xy+{x}^{2}{y}^{2}\right)f\left(y\right),y=-1..1\right)\right)=d\left(x\right)$
 ${\mathrm{eq4}}{≔}{f}{}\left({x}\right){+}{b}{}\left({{\int }}_{{-1}}^{{1}}\left({{x}}^{{2}}{}{{y}}^{{2}}{+}{x}{}{y}\right){}{f}{}\left({y}\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{y}\right){=}{d}{}\left({x}\right)$ (12)
 > $\mathrm{intsolve}\left(\mathrm{eq4},f\left(x\right)\right)$
 ${f}{}\left({x}\right){=}{{\int }}_{{-1}}^{{1}}\left({-}\frac{{5}{}{{x}}^{{2}}{}{b}{}{{y}}^{{2}}}{{2}{}{b}{+}{5}}{-}\frac{{3}{}{x}{}{b}{}{y}}{{2}{}{b}{+}{3}}\right){}{d}{}\left({y}\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{y}{+}{d}{}\left({x}\right)$ (13)

Example 5

 An example where an arbitrary constant $\mathrm{_C1}$ is introduced with the solution:
 > $\mathrm{eq5}≔3{x}^{2}+4x=\mathrm{Int}\left(\left(6{x}^{2}y+4x{y}^{2}\right)p\left(y\right),y=-1..1\right)$
 ${\mathrm{eq5}}{≔}{3}{}{{x}}^{{2}}{+}{4}{}{x}{=}{{\int }}_{{-1}}^{{1}}\left({6}{}{{x}}^{{2}}{}{y}{+}{4}{}{x}{}{{y}}^{{2}}\right){}{p}{}\left({y}\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{y}$ (14)
 > $\mathrm{intsolve}\left(\mathrm{eq5},p\left(x\right)\right)$
 ${p}{}\left({x}\right){=}{{x}}^{{2}}{}{\mathrm{_C1}}{+}\frac{{3}}{{4}}{}{x}{+}\frac{{3}}{{2}}{-}\frac{{3}}{{5}}{}{\mathrm{_C1}}$ (15)

Example 6

 An example which is in the form of an integral transform:
 > $\mathrm{eq6}≔\mathrm{Int}\left(\frac{f\left(y\right)}{1+{y}^{2}}\mathrm{sin}\left(xy\right),y=0..\mathrm{\infty }\right)=\mathrm{exp}\left(-3x\right)$
 ${\mathrm{eq6}}{≔}{{\int }}_{{0}}^{{\mathrm{\infty }}}\frac{{f}{}\left({y}\right){}{\mathrm{sin}}{}\left({x}{}{y}\right)}{{{y}}^{{2}}{+}{1}}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{y}{=}{{ⅇ}}^{{-}{3}{}{x}}$ (16)
 > $\mathrm{intsolve}\left(\mathrm{eq6},f\left(x\right)\right)$
 ${f}{}\left({x}\right){=}\frac{{2}{}{x}{}\left({{x}}^{{2}}{+}{1}\right)}{{\mathrm{\pi }}{}\left({{x}}^{{2}}{+}{9}\right)}$ (17)

Example 7

 The solution of a Fredholm integral equation can be approximated using the collocation method:
 > $a≔0:$
 > $b≔5:$
 > $g≔x↦-x:$
 > $k≔\left(x,y\right)↦\mathrm{piecewise}\left(x
 > $\mathrm{eq7}≔f\left(x\right)+\mathrm{int}\left(f\left(y\right)k\left(x,y\right),y=a..b\right)+g\left(x\right)$
 ${\mathrm{eq7}}{≔}{f}{}\left({x}\right){+}{{\int }}_{{0}}^{{5}}\left(\left\{\begin{array}{cc}{x}{}\left({5}{-}{y}\right)& {x}{<}{y}\\ {y}{}\left({5}{-}{x}\right)& {\mathrm{otherwise}}\end{array}\right\\right){}{f}{}\left({y}\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{y}{-}{x}$ (18)
 > $p≔\mathrm{unapply}\left(\mathrm{rhs}\left(\mathrm{intsolve}\left(\mathrm{eq7},f\left(x\right),\mathrm{method}=\mathrm{collocation},\mathrm{order}=10\right)\right),x\right)$
 ${p}{≔}{x}{↦}\frac{{220000}}{{8122882553}}{\cdot }{{x}}^{{10}}{-}\frac{{33081098947000}}{{60853654989179049}}{\cdot }{{x}}^{{9}}{+}\frac{{98331769350000}}{{20284551663059683}}{\cdot }{{x}}^{{8}}{-}\frac{{501031314879600}}{{20284551663059683}}{\cdot }{{x}}^{{7}}{+}\frac{{1594256521707500}}{{20284551663059683}}{\cdot }{{x}}^{{6}}{-}\frac{{6502492915408455}}{{40569103326119366}}{\cdot }{{x}}^{{5}}{+}\frac{{8429866033778125}}{{40569103326119366}}{\cdot }{{x}}^{{4}}{-}\frac{{9902325865977605}}{{60853654989179049}}{\cdot }{{x}}^{{3}}{+}\frac{{1401159656806875}}{{20284551663059683}}{\cdot }{{x}}^{{2}}{-}\frac{{231864099964442}}{{20284551663059683}}{\cdot }{x}$ (19)
 Plot of the approximate solution:
 > $\mathrm{plot}\left(p,a..b,\mathrm{title}="Approximate Solution of Integral Equation",\mathrm{color}=\mathrm{blue}\right)$
 To test how close the approximation matches the actual solution, we can sample the absolute value of the expression for the integral equation over the interval:
 > $n≔100:$
 > $q≔\mathrm{unapply}\left(\mathrm{abs}\left(\mathrm{eval}\left(\mathrm{eq7},f=p\right)\right),x\right):$
 > $X≔\mathrm{Vector}\left(n,i↦\mathrm{evalf}\left(a+\frac{\left(i-1\right)\cdot \left(b-a\right)}{n-1}\right)\right):$
 > $E≔\mathrm{~}\left[q\right]\left(X\right):$
 > $\mathrm{dataplot}\left(X,E,\mathrm{style}=\mathrm{line},\mathrm{title}="Error of Approximate Solution of Integral Equation"\right)$
 > $\mathrm{avg_error}≔\frac{\mathrm{add}\left(E\right)}{n}$
 ${\mathrm{avg_error}}{≔}{0.0001670862960}$ (20)
 > $\mathrm{max_error}≔\mathrm{max}\left(E\right)$
 ${\mathrm{max_error}}{≔}{0.001720013949}$ (21)

Example 8

 The Laplace method can be used to solve singular integral equations:
 > $\mathrm{eq8}≔\mathrm{Int}\left(\frac{f\left(y\right)}{\mathrm{sqrt}\left(x-y\right)},y=0..x\right)-{x}^{2}$
 ${\mathrm{eq8}}{≔}{{\int }}_{{0}}^{{x}}\frac{{f}{}\left({y}\right)}{\sqrt{{x}{-}{y}}}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{y}{-}{{x}}^{{2}}$ (22)
 > $\mathrm{intsolve}\left(\mathrm{eq8},f\left(x\right),\mathrm{method}=\mathrm{Laplace}\right)$
 ${f}{}\left({x}\right){=}\frac{{8}{}{{x}}^{{3}}{{2}}}}{{3}{}{\mathrm{\pi }}}$ (23)

References

 Chambers, LL. G. Integral Equations: A Short Course. London: International Textbook Company Ltd., 1976.
 Tricomi, F. G. Integral Equations. New York: Interscience, 1957.

Compatibility

 • The intsolve command was updated in Maple 2021.