• Views 888
• Citations 0
• ePub 21
• PDF 645
`Mathematical Problems in EngineeringVolume 2014, Article ID 147497, 6 pageshttp://dx.doi.org/10.1155/2014/147497`
Research Article

## Bezier Curves for Solving Fredholm Integral Equations of the Second Kind

2The Center of Excellence on Modelling and Control Systems (CEMCS), Mashhad, Iran
3Department of Mathematics and Institute for Mathematical Research, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia

Received 10 July 2013; Accepted 11 October 2013; Published 5 January 2014

Copyright © 2014 F. Ghomanjani et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The Bezier curves are presented to estimate the solution of the linear Fredholm integral equation of the second kind. A direct algorithm for solving this problem is given. We have chosen the Bezier curves as piecewise polynomials of degree and determine Bezier curves on [0, 1] by control points. Numerical examples illustrate that the algorithm is applicable and very easy to use.

#### 1. Introduction

Integral equations are often involved in the mathematical formulation of physical phenomena, and they can be encountered in various fields of science such as physics [1], biology [2], and engineering (see [3, 4]). But we can also use it in numerous applications, such as control, biomechanics, elasticity, economics, electrical engineering, electrodynamics, electrostatics, fluid dynamics, game theory, heat and mass transfer, medicine, oscillation theory, plasticity, and queuing theory [5]. Fredholm integral equations of the second kind are shown in studies which include airfoil theory [6], elastic contact problems (see [7, 8]), fracture mechanics [9], combined infrared radiation, and molecular conduction [10]. Many different basic functions have been used to estimate the solution of integral equations, such as orthogonal functions and wavelets (see [11, 12]). Depending on the structure, the orthogonal functions may be widely classified into three families [13]. The first includes sets of piecewise constant orthogonal functions (e.g., Walsh, block-pulse, Haar, etc.). The second consists of sets of orthogonal polynomials (e.g., Laguerre, Legendre, Chebyshev, etc.). The third are the widely used sets of sine-cosine functions in the Fourier series. Fredholm integral equations of the second kind are much more difficult to solve than ordinary differential equations. Therefore, many authors have tried various transform methods to overcome these difficulties (see [11, 12]). Recently, hybrid functions have been applied extensively for solving differential equations or systems, and they proved to be a useful mathematical tool. The pioneering work in system analysis via hybrid functions was led in [14, 15], who first derived an operational matrix for the integrals of the hybrid function vector and paved the way for the hybrid function analysis of the dynamic systems. But they derived the matrix of small order, and the calculations are not enough to achieve high accuracy. Hsiao [16] presented the properties of hybrid functions which consist of block-pulse functions plus the Legendre polynomials. Based upon some useful properties of hybrid functions, integration of the cross product, a special product matrix and a related coefficient matrix with optimal order are applied to solve these integral equations. The main characteristic of this technique is to convert an integral equation into an algebraic one. Maleknejad and Mahmoudi [17] used a simple base, a combination of block-pulse functions on and the Taylor polynomials, that is called the hybrid Taylor block-pulse functions, to solve the linear Fredholm integral equation of the second kind. One of the advantages of this method is that the coefficients of expansion of each function in this base could be computed directly without estimation.

Consider the following integral equation: where , , is an unknown function.

In this paper, we discuss a technique similar to that used in [18] for solving integral equations by using the Bezier control points. There are many papers and books that deal with the Bezier curves or surface techniques. Harada and Nakamae [19] and Nürnberger and Zeilfelder [20] used the Bezier control points in approximated data and functions. Zheng et al. [21] proposed the use of the control points of the Bernstein-Bezier form for solving differential equations numerically, and also Evrenosoglu and Somali [18] used this approach for solving singular-perturbed two-point boundary value problems. The Bezier curves are used in solving partial differential equations; besides, Wave and Heat equations are solved in Bezier form (see [2225] ). Wu [26] presented the least squares method for solving partial differential equations on arbitrary polygon domain by the Bezier control points. Wu [26] used triangular Bezier patches of degree with continuity to approximate the exact solution of partial differential equations. Bezier curves are used for solving dynamical systems (see [27]), also the Bezier control points method is used for solving delay differential equation (see [28]). Some other applications of the Bezier functions and control points are found in ([2931]), that are used in computer-aided geometric design and image compression.

The use of the Bezier curves for solving Fredholm integral equations of the second kind is a novel idea. Although the method is very easy to be used and straightforward, the obtained results are satisfactory (see the numerical results).

We suggest a technique similar to that used in [28] for solving Fredholm integral equations of the second kind. The current paper is organized as follows.

Presented algorithm will be stated in Section 2. In Section 3, the convergence analysis will be presented. Some numerical examples are solved in Section 4 which show the efficiency and reliability of the method. Finally, Section 5 will give a conclusion in brief.

#### 2. The Algorithm

Our strategy is to use Bezier curves to approximate the solutions by where is given below. Define the Bezier polynomial of degree that approximates the values of over the interval as follows: where ; is the Bernstein polynomial of degree over the interval , and is the control point (see [21]). By substituting (2) in (1), one may define for as follows:

In Section 3, the convergence of this method is proven by Bezier curves when the degree of the approximate solution, , tends to infinity.

Now, we define the residual function over the interval as follows: where is the Euclidean norm. Our aim is to solve the following problem over the interval : When the minimization problem (6) is posed, the condition is equivalent to fix the first control point . The mathematical programming problem (6) can be solved by many subroutine algorithms, and we used Maple to solve this optimization problem.

#### 3. Convergence Analysis

In this section without the loss of generality, we analyze the convergence of the control-point-based method when applied to the integral equation (1) with the time interval . So, the following problem is considered: where is a given real number and and are known functions for .

Lemma 1. For a polynomial in Bezier form we have where is the Bezier coefficient of after being degree-elevated to degree .

Proof . See [21]

The convergence of the approximate solution could be done in degree raising of the Bezier polynomial approximation.

Theorem 2. If the integral equation (7) has a unique continuous solution , then the approximate solution obtained by the control-point-based method converges to the exact solution as the degree of the approximate solution tends to infinity.

Proof. Given an arbitrary small positive number , by the Weierstrass Theorem (see [32]), one can easily find polynomial of degree such that , where stands for the -norm over . In particular, we have In general, does not satisfy the boundary conditions. After a small perturbation with constant polynomial , for , we can obtain polynomial such that satisfies the boundary condition . Thus, . By using (10), one has We have Now, let for every . Thus, for , one may find an upper bound for the following residual: where is a constant.
Since the residual can be considered as a polynomial because if it is not a polynomial, we can use the Taylor series for it, we can represent the statement by a Bezier form. Thus, we have Then, by Lemma 1, there exists an integer such that when , we have which gives Suppose that is an approximated solution of (7) obtained by the control-point-based method of degree . Let Define the following norm for the difference-approximated solution and the exact solution : It is easy to show that The last inequality in (20) is obtained by Lemma 1 in which is a constant positive number. Now, by Lemma 1 and (15), it can be shown that where the last inequality in (21) is coming from (17).
Thus, from (21) we have Since the infinite norm and the norm defined in (19) are equivalent, there is a where
This completes the proof.

#### 4. Numerical Examples

In this Section, we present some test problems and apply the method presented in this paper for solving them. The well-known symbolic software “Maple ” has been employed for calculations.

Example 1. Consider the integral equation described by For this integral equation, there exists the exact solution (see [17]). With the method described in the paper, and with one can find the following approximate solution: Figure 1 shows the value of error for Example 1 where the maximum error of hybrid Taylor and Block-Pulse functions [17] is for and .

Figure 1: The error for Example 1.

Example 2. Consider the following integral equation: , (see [17]). With the described method and with , one can find the following solution: The error curve of Example 2 is shown in Figure 2 where the maximum error of hybrid Taylor and Block-Pulse functions is for and (see [17]). In Table 1, analytic, numerical results of the presented method, and the absolute error of the presented method are shown, respectively.

Table 1: Exact and estimated values of for Example 2.
Figure 2: The error for Example 2.

Example 3. Consider the following integral equation (see [17]): (see [17]). With the described method and with , one can find the following approximate solution: Figure 3 shows the value of error for Example 3 where the maximum error of hybrid Taylor and Block-Pulse functions [17] is for and .

Figure 3: The error for Example 3.

Example 4. Consider the nonlinear Fredholm integral equation with exact solution (see [33]): In our method, with one can find the following solution Figure 4 shows the value of error for Example 4.

Figure 4: The error for Example 4.

#### 5. Conclusions

A simple and effective algorithm based on Bezier curves is presented for solving Fredholm integral equations of the second kind. The method is computationally attractive and also reduces the CPU time and the computer memory while at the same time keeping the accuracy of the solution.

#### Conflicts of Interests

The authors declare that they have no conflicts of interests regarding publication of this article.

#### Acknowledgments

The authors are very grateful to the referees for their valuable suggestions and comments that improved the paper. The third author acknowledges that this research was partially supported by the Universiti Putra Malaysia under the research Grant ERGS 1-2013(5527179).

#### References

1. F. Bloom, “Asymptotic bounds for solutions to a system of damped integro-differential equations of electromagnetic theory,” Journal of Mathematical Analysis and Applications, vol. 73, no. 2, pp. 524–542, 1980.
2. K. Holmåker, “Global asymptotic stability for a stationary solution of a system of integro-differential equations describing the formation of liver zones,” SIAM Journal on Mathematical Analysis, vol. 24, no. 1, pp. 116–128, 1993.
3. M. A. Abdou, “On a symptotic methods for Fredholm-Volterra integral equation of the second kind in contact problems,” Journal of Computational and Applied Mathematics, vol. 154, no. 2, pp. 431–446, 2003.
4. L. K. Forbes, S. Crozier, and D. M. Doddrell, “Calculating current densities and fields produced by shielded magnetic resonance imaging probes,” SIAM Journal on Applied Mathematics, vol. 57, no. 2, pp. 401–425, 1997.
5. A. D. Polyanin and A. V. Manzhirov, Handbook of Integral Equations, Chapman & Hall/CRC Press, Boca Raton, Fla, USA, 2nd edition, 2008.
6. M. A. Golberg, “The convergence of a collocation method for a class of Cauchy singular integral equations,” Journal of Mathematical Analysis and Applications, vol. 100, no. 2, pp. 500–512, 1984.
7. E. V. Kovalenko, “Some approximate methods for solving integral equations for mixed problems,” Journal of Applied Mathematics and Mechanics, vol. 53, no. 1, pp. 85–92, 1989.
8. B. I. Smetanin, “On an integral equation of axisymmetric problems for an elastic body containing an inclusion,” Journal of Applied Mathematics and Mechanics, vol. 55, no. 3, pp. 371–375, 1991.
9. J. R. Willis and S. Nemat-Nasser, “Singular perturbation solution of a class of singular integral equations,” Quarterly of Applied Mathematics, vol. 48, no. 4, pp. 741–753, 1990.
10. J. I. Frankel, “A Galerkin solution to a regularized Cauchy singular integro-differential equation,” Quarterly of Applied Mathematics, vol. 53, no. 2, pp. 245–258, 1995.
11. W. F. Blyth, R. L. May, and P. Widyaningsih, “Volterra integral equations solved in Fredholm form using Walsh functions,” The ANZIAM Journal, vol. 45, pp. C269–C282, 2004.
12. M. H. Reihani and Z. Abadi, “Rationalized Haar functions method for solving Fredholm and Volterra integral equations,” Journal of Computational and Applied Mathematics, vol. 200, no. 1, pp. 12–20, 2007.
13. K. B. Datta and B. M. Mohan, Orthogonal Functions in Systems and Control, vol. 9, World Scientific, River Edge, NJ, USA, 1995.
14. H. R. Marzban and M. Razzaghi, “Numerical solution of the controlled duffing oscillator by hybrid functions,” Applied Mathematics and Computation, vol. 140, no. 2-3, pp. 179–190, 2003.
15. M. Razzaghi and H.-R. Marzban, “A hybrid analysis direct method in the calculus of variations,” International Journal of Computer Mathematics, vol. 75, no. 3, pp. 259–269, 2000.
16. C.-H. Hsiao, “Hybrid function method for solving Fredholm and Volterra integral equations of the second kind,” Journal of Computational and Applied Mathematics, vol. 230, no. 1, pp. 59–68, 2009.
17. K. Maleknejad and Y. Mahmoudi, “Numerical solution of linear Fredholm integral equation by using hybrid Taylor and block-pulse functions,” Applied Mathematics and Computation, vol. 149, no. 3, pp. 799–806, 2004.
18. M. Evrenosoglu and S. Somali, “Least squares methods for solving singularly perturbed two-point boundary value problems using Bézier control points,” Applied Mathematics Letters, vol. 21, no. 10, pp. 1029–1032, 2008.
19. K. Harada and E. Nakamae, “Application of the Bézier curve to data interpolation,” Computer-Aided Design, vol. 14, no. 1, pp. 55–59, 1982.
20. G. Nürnberger and F. Zeilfelder, “Developments in bivariate spline interpolation,” Journal of Computational and Applied Mathematics, vol. 121, no. 1-2, pp. 125–152, 2000.
21. J. Zheng, T. W. Sederberg, and R. W. Johnson, “Least squares methods for solving differential equations using Bézier control points,” Applied Numerical Mathematics, vol. 48, no. 2, pp. 237–252, 2004.
22. J. V. Beltran and J. Monterde, “Bézier solutions of the wave equation,” in Computational Science and Its Applications—ICCSA, vol. 3044 of Lecture Notes in Computer Science, pp. 631–640, Springer, Berlin, Germany, 2004.
23. R. Cholewa, A. J. Nowak, R. A. Bialecki, and L. C. Wrobel, “Cubic Bezier splines for BEM heat transfer analysis of the 2-D continuous casting problems,” Computational Mechanics, vol. 28, no. 3-4, pp. 282–290, 2002.
24. C.-H. Chu, C. C. L. Wang, and C.-R. Tsai, “Computer aided geometric design of strip using developable Bézier patches,” Computers in Industry, vol. 59, no. 6, pp. 601–611, 2008.
25. A. T. Layton and M. Van de Panne, “A numerically efficient and stable algorithm for animating water waves,” The Visual Computer, vol. 18, no. 1, pp. 41–53, 2002.
26. J. Wu, “Least squares methods for solving partial differential equations by using Bézier control points,” Applied Mathematics and Computation, vol. 219, no. 8, pp. 3655–3663, 2012.
27. M. Gachpazan, “Solving of time varying quadratic optimal control problems by using Bézier control points,” Computational & Applied Mathematics, vol. 30, no. 2, pp. 367–379, 2011.
28. F. Ghomanjani and M. H. Farahi, “The Bezier control points method for solving delay differential equation,” Intelligent Control and Automation, vol. 3, no. 2, pp. 188–196, 2012.
29. G. Farin, Curves and Surfaces for Computer Aided Geometric Design, Academic Press, New York, NY, USA, 1st edition, 1988.
30. F. Ghomanjani, M. H. Farahi, and M. Gachpazan, “Bézier control points method to solve constrained quadratic optimal control of time varying linear systems,” Computational & Applied Mathematics, vol. 31, no. 3, pp. 433–456, 2012.
31. Y. Q. Shi and H. Sun, Image and Video Compression for Multimedia Engineering, CRC Press, 2000.
32. W. Rudin, Principles of Mathematical Analysis, McGraw-Hill, 1986.
33. S. M. Mirzaei, “Homotopy perturbation method for solving the second kind of non-linear integral equations,” International Mathematical Forum, vol. 5, no. 21–24, pp. 1149–1154, 2010.