<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">OJAppS</journal-id><journal-title-group><journal-title>Open Journal of Applied Sciences</journal-title></journal-title-group><issn pub-type="epub">2165-3917</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ojapps.2017.74014</article-id><article-id pub-id-type="publisher-id">OJAppS-75814</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biomedical&amp;Life Sciences</subject><subject> Chemistry&amp;Materials Science</subject><subject> Computer Science&amp;Communications</subject><subject> Engineering</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Numerical Solutions for Quadratic Integro-Differential Equations of Fractional Orders
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Fatheah</surname><given-names>Alhendi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Wafa</surname><given-names>Shammakh</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Hind</surname><given-names>Al-Badrani</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Taibah University, Madinah, Saudi Arabia</addr-line></aff><aff id="aff1"><addr-line>Department of Mathematics, King Abdulaziz University, Jeddah, Saudi Arabia</addr-line></aff><pub-date pub-type="epub"><day>25</day><month>04</month><year>2017</year></pub-date><volume>07</volume><issue>04</issue><fpage>157</fpage><lpage>170</lpage><history><date date-type="received"><day>March</day>	<month>5,</month>	<year>2017</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>April</month>	<year>25,</year>	</date><date date-type="accepted"><day>April</day>	<month>30,</month>	<year>2017</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  In this article, variational iteration method (VIM) and homotopy perturbation method (HPM) solve the nonlinear initial value problems of first-order fractional quadratic integro-differential equations (FQIDEs). We use the Caputo sense in this article to describe the fractional derivatives. The solutions of the problems are derived by infinite convergent series, and the results show that both methods are most convenient and effective.
 
</p></abstract><kwd-group><kwd>Fractional Quadratic Integro-Differential Equations</kwd><kwd> Variational Iteration Method</kwd><kwd> Homotopy Perturbation Method</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The fractional calculus has appeared in many areas during the recent decades. Some scientists use approximation and numerical methods because there are almost no exact solutions of the fractional differential equations. He has pro- posed the VIM and HPM to solve the problems of linear and nonlinear [<xref ref-type="bibr" rid="scirp.75814-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.75814-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.75814-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.75814-ref4">4</xref>] . VIM is based on Lagrange multiplier. The another method is HPM which defines as a coupling of the traditional perturbation method and homotopy in topology. Many authors successfully apply these methods to find the solutions of functional equations which arise in scientific and engineering problems [<xref ref-type="bibr" rid="scirp.75814-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.75814-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.75814-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.75814-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.75814-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.75814-ref9">9</xref>] . The Adomian decomposition method presents solution of functional equations but perhaps we find some difficulties that will arise during the computation of Adomian polynomials, the VIM and HPM overcome it is difficulties [<xref ref-type="bibr" rid="scirp.75814-ref7">7</xref>] . Fractional differential equations have diverse applications of physical phenomena [<xref ref-type="bibr" rid="scirp.75814-ref10">10</xref>] [<xref ref-type="bibr" rid="scirp.75814-ref11">11</xref>] [<xref ref-type="bibr" rid="scirp.75814-ref12">12</xref>] [<xref ref-type="bibr" rid="scirp.75814-ref13">13</xref>] , for instance, acoustics, electromag- netism, control theory, robotics, viscoelastic materials, diffusion, edge detection, turbulence, signal processing, anomalous diffusion and fractured media [<xref ref-type="bibr" rid="scirp.75814-ref14">14</xref>] . In literature, Momani and Noor [<xref ref-type="bibr" rid="scirp.75814-ref15">15</xref>] used the Adomian decomposition method for solving the fourth order fractional integro-differential Equation. Elbeleze et al. [<xref ref-type="bibr" rid="scirp.75814-ref16">16</xref>] [<xref ref-type="bibr" rid="scirp.75814-ref17">17</xref>] , Kadem and Kilicman [<xref ref-type="bibr" rid="scirp.75814-ref18">18</xref>] applied the HPM and VIM methods for integro-differential Equation of fractional order with initial-boundary con- ditions. Recently, Gaafar [<xref ref-type="bibr" rid="scirp.75814-ref19">19</xref>] studied the existence and nondecreasing solution for the initial value problem of a quadratic integro-differential equations. How- ever, there is little work on nonlinear fractional quadratic integro-differential equations.</p><p>Our goal for this article is extending the analysis of VIM and HPM to construct the approximate solutions of the following nonlinear initial value problems for first-order fractional quadratic integro-differential equations.</p><disp-formula id="scirp.75814-formula364"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x2.png"  xlink:type="simple"/></disp-formula><p>subject to the following initial condition:</p><disp-formula id="scirp.75814-formula365"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x3.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x4.png" xlink:type="simple"/></inline-formula>is the fractional derivative in the caputo sense, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x5.png" xlink:type="simple"/></inline-formula>is any nonlinear function, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x6.png" xlink:type="simple"/></inline-formula>is real constant and g is given and can be approximated by taylor polynomials.</p></sec><sec id="s2"><title>2. Basic Definitions</title><p>In this section, we intend to present some basic definitions and properties of fractional calculus theory which are further used in this article.</p><p>Definition 1.</p><p>Areal function <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x7.png" xlink:type="simple"/></inline-formula> is said to be in space<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x8.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x9.png" xlink:type="simple"/></inline-formula>if there exists a real number<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x10.png" xlink:type="simple"/></inline-formula>, such that<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x11.png" xlink:type="simple"/></inline-formula>, where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x12.png" xlink:type="simple"/></inline-formula>, and it is said to be in the space <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x13.png" xlink:type="simple"/></inline-formula> if and only if<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x14.png" xlink:type="simple"/></inline-formula>,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x15.png" xlink:type="simple"/></inline-formula>.</p><p>Definition 2.</p><p>The Riemann-Liouville fractional integral operator of order <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x16.png" xlink:type="simple"/></inline-formula> of a function<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x17.png" xlink:type="simple"/></inline-formula>, is defined as</p><disp-formula id="scirp.75814-formula366"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x18.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula367"><graphic  xlink:href="http://html.scirp.org/file/5-2310712x19.png"  xlink:type="simple"/></disp-formula><p>Some properties of the operator <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x20.png" xlink:type="simple"/></inline-formula> can be found in [<xref ref-type="bibr" rid="scirp.75814-ref11">11</xref>] , which are needed here, as follows:</p><p>for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x21.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x22.png" xlink:type="simple"/></inline-formula>:</p><p>1. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x23.png" xlink:type="simple"/></inline-formula></p><p>2. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x24.png" xlink:type="simple"/></inline-formula></p><p>3. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x25.png" xlink:type="simple"/></inline-formula></p><p>Definition 3. The fractional derivative of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x26.png" xlink:type="simple"/></inline-formula> in the caputo sense is defined as</p><disp-formula id="scirp.75814-formula368"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x27.png"  xlink:type="simple"/></disp-formula><p>for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x28.png" xlink:type="simple"/></inline-formula>.</p><p>Lemma 1.</p><p>If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x29.png" xlink:type="simple"/></inline-formula> the the following two properties hold:</p><p>1. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x30.png" xlink:type="simple"/></inline-formula></p><p>2. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x31.png" xlink:type="simple"/></inline-formula></p></sec><sec id="s3"><title>3. Analysis of VIM</title><p>The basic concept of the VIM is constructing the correction functional for the fractional quadratic integro-differential equation sees Equation (1) with initial conditions,</p><disp-formula id="scirp.75814-formula369"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x32.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula370"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x33.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x34.png" xlink:type="simple"/></inline-formula>is the Rieman-Liouville fractional integral operator of order<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x35.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x36.png" xlink:type="simple"/></inline-formula>is a general Lagrange multiplier and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x37.png" xlink:type="simple"/></inline-formula> refers to the restricted variation (i.e.) <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x38.png" xlink:type="simple"/></inline-formula>to identify the approximate Lagrange multiplier, construct the correctional function (6) which can be approximately expressed as:</p><disp-formula id="scirp.75814-formula371"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x39.png"  xlink:type="simple"/></disp-formula><p>taking the variation of Equation (7) to the independent variable <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x40.png" xlink:type="simple"/></inline-formula> we find</p><disp-formula id="scirp.75814-formula372"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x41.png"  xlink:type="simple"/></disp-formula><p>to make the previous equation stationary, we gain the following stationary conditions:</p><disp-formula id="scirp.75814-formula373"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x42.png"  xlink:type="simple"/></disp-formula><p>finally, the Lagrange multiplier is:</p><disp-formula id="scirp.75814-formula374"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x43.png"  xlink:type="simple"/></disp-formula><p>We achieve the following iteration formula by substitution of (10) into the functional (6)</p><disp-formula id="scirp.75814-formula375"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x44.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula376"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x45.png"  xlink:type="simple"/></disp-formula><p>the initial approximation <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x46.png" xlink:type="simple"/></inline-formula> can be selected by the following way which satisfies initial conditions</p><disp-formula id="scirp.75814-formula377"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x47.png"  xlink:type="simple"/></disp-formula></sec><sec id="s4"><title>4. Analysis of HPM</title><p>The main concept of the HPM is constructing the homotopy for fractional quadratic integro-differential equation sees Equation (1),</p><disp-formula id="scirp.75814-formula378"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x48.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula379"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x49.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x50.png" xlink:type="simple"/></inline-formula>is an embedding parameter. If<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x51.png" xlink:type="simple"/></inline-formula>, then Equation (15) turns into a linear Equation</p><disp-formula id="scirp.75814-formula380"><graphic  xlink:href="http://html.scirp.org/file/5-2310712x52.png"  xlink:type="simple"/></disp-formula><p>and when<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x53.png" xlink:type="simple"/></inline-formula>, then Equation (15) becomes to be the original problem.</p><p>The solution of Equation (1) can be considered as a power series in p which is the basic assumption of HPM :</p><disp-formula id="scirp.75814-formula381"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x54.png"  xlink:type="simple"/></disp-formula><p>when <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x55.png" xlink:type="simple"/></inline-formula> in (16) the approximate solution of Equation (1) can be as following</p><disp-formula id="scirp.75814-formula382"><graphic  xlink:href="http://html.scirp.org/file/5-2310712x56.png"  xlink:type="simple"/></disp-formula><p>First, substitute the relation (16) in the Equation (15). Second, equate the terms which have the same power’s of p which yield to the following series of equations:</p><disp-formula id="scirp.75814-formula383"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x57.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula384"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x58.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula385"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x59.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula386"><label>(20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x60.png"  xlink:type="simple"/></disp-formula><p>and so on, the functions <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x61.png" xlink:type="simple"/></inline-formula> satisfy the following condition:</p><disp-formula id="scirp.75814-formula387"><graphic  xlink:href="http://html.scirp.org/file/5-2310712x62.png"  xlink:type="simple"/></disp-formula></sec><sec id="s5"><title>5. Applications</title><p>n this section, we apply VIM and HPM to first-order nonlinear (FQIDEs).</p><p>Example 1.</p><p>Consider the following nonlinear first-order (FQIDEs):</p><disp-formula id="scirp.75814-formula388"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x63.png"  xlink:type="simple"/></disp-formula><p>subject to the following initial condition</p><disp-formula id="scirp.75814-formula389"><label>(22)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x64.png"  xlink:type="simple"/></disp-formula><p>According to VIM, the expression of the iteration formula (12) for Equation (21) can be observed in the following form:</p><disp-formula id="scirp.75814-formula390"><label>(23)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x65.png"  xlink:type="simple"/></disp-formula><p>To avoid the difficulty of fractional integration, for the exponential term we</p><p>take the truncated taylor expansion in (23), e.g., <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x66.png" xlink:type="simple"/></inline-formula>to sa-</p><p>tisfy the initial condition (22), we assume that the initial approximation has the following form <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x67.png" xlink:type="simple"/></inline-formula> first-order approximation takes the following form by using iteration Formula (23):</p><disp-formula id="scirp.75814-formula391"><label>(24)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x68.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula392"><label>(25)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x69.png"  xlink:type="simple"/></disp-formula><p><xref ref-type="table" rid="table1">Table 1</xref> and <xref ref-type="fig" rid="fig1">Figure 1</xref> presents the approximate solution for the different values of a, we have noticed that the accuracy is improving. First, by computing more terms of the approximate solutions. The second way is taking more terms in the taylor expansion of the exponential term.</p><p>According to HPM, we build the following homotopy:</p><disp-formula id="scirp.75814-formula393"><label>(26)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x70.png"  xlink:type="simple"/></disp-formula><p>First, substitute the relation (16) in the Equation (26).</p><p>Second, equate the terms which have the same power’s of p which yield to the following series of Equations:</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Approximate solution for Equation (21) is ob- tained by VIM with different values of a</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2310712x71.png"/></fig><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Approximate solution for Equation (21) at different values of a</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >x</th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x72.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x73.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x74.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x75.png" xlink:type="simple"/></inline-formula></th></tr></thead><tr><td align="center" valign="middle" >0.10</td><td align="center" valign="middle" >3.46688592</td><td align="center" valign="middle" >2.27052430</td><td align="center" valign="middle" >1.65826130</td><td align="center" valign="middle" >1.33251002</td></tr><tr><td align="center" valign="middle" >0.20</td><td align="center" valign="middle" >4.93376139</td><td align="center" valign="middle" >3.18886968</td><td align="center" valign="middle" >2.27376989</td><td align="center" valign="middle" >1.74152708</td></tr><tr><td align="center" valign="middle" >0.30</td><td align="center" valign="middle" >6.70750098</td><td align="center" valign="middle" >4.28829701</td><td align="center" valign="middle" >3.00556169</td><td align="center" valign="middle" >2.24764961</td></tr><tr><td align="center" valign="middle" >0.40</td><td align="center" valign="middle" >8.84373442</td><td align="center" valign="middle" >5.65002513</td><td align="center" valign="middle" >3.90802714</td><td align="center" valign="middle" >2.87616667</td></tr><tr><td align="center" valign="middle" >0.50</td><td align="center" valign="middle" >11.37418953</td><td align="center" valign="middle" >7.33008584</td><td align="center" valign="middle" >5.02934761</td><td align="center" valign="middle" >3.65730794</td></tr><tr><td align="center" valign="middle" >0.60</td><td align="center" valign="middle" >14.32676071</td><td align="center" valign="middle" >9.37815747</td><td align="center" valign="middle" >6.41827073</td><td align="center" valign="middle" >4.62649375</td></tr><tr><td align="center" valign="middle" >0.70</td><td align="center" valign="middle" >17.72917203</td><td align="center" valign="middle" >11.84228561</td><td align="center" valign="middle" >8.12594765</td><td align="center" valign="middle" >5.82458502</td></tr><tr><td align="center" valign="middle" >0.80</td><td align="center" valign="middle" >21.60990786</td><td align="center" valign="middle" >14.77061532</td><td align="center" valign="middle" >10.20665676</td><td align="center" valign="middle" >7.29813333</td></tr><tr><td align="center" valign="middle" >0.90</td><td align="center" valign="middle" >25.99846722</td><td align="center" valign="middle" >18.21217939</td><td align="center" valign="middle" >12.71816657</td><td align="center" valign="middle" >9.09963086</td></tr><tr><td align="center" valign="middle" >1.00</td><td align="center" valign="middle" >30.92541601</td><td align="center" valign="middle" >22.21731293</td><td align="center" valign="middle" >15.72195568</td><td align="center" valign="middle" >11.28776041</td></tr></tbody></table></table-wrap><disp-formula id="scirp.75814-formula394"><label>(27)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x76.png"  xlink:type="simple"/></disp-formula><p>and so on, apply the operator <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x77.png" xlink:type="simple"/></inline-formula> to the previous equations, and use the initial condition (22), to gain the following equations:</p><disp-formula id="scirp.75814-formula395"><label>(28)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x78.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula396"><label>(29)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x79.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula397"><label>(30)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x80.png"  xlink:type="simple"/></disp-formula><p>and so on, by taking the truncated taylor expansions for the exponential term in</p><p>(29, 30): e.g., <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x81.png" xlink:type="simple"/></inline-formula>to avoid the difficulty of fractional integra-</p><p>tion, thus by solving Equations (28, 29, 30), we obtain <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x82.png" xlink:type="simple"/></inline-formula></p><disp-formula id="scirp.75814-formula398"><label>(31)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x83.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula399"><label>(32)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x84.png"  xlink:type="simple"/></disp-formula><p>The two terms approximation are formed as the following Equation</p><disp-formula id="scirp.75814-formula400"><label>(33)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x85.png"  xlink:type="simple"/></disp-formula><p><xref ref-type="table" rid="table2">Table 2</xref> and <xref ref-type="fig" rid="fig2">Figure 2</xref> shows the approximate solutions of (33) for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x86.png" xlink:type="simple"/></inline-formula> and for some values of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x86.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x87.png" xlink:type="simple"/></inline-formula>.</p><p><xref ref-type="fig" rid="fig3">Figure 3</xref> represent a comparison between two approximate solutions by using VIM and HPM methods.</p><p>Example 2.</p><p>Consider the following (FQIDEs):</p><disp-formula id="scirp.75814-formula401"><label>(34)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x88.png"  xlink:type="simple"/></disp-formula><p>According to VIM, the expression of the iteration Formula (12) for Equation (34) can be observed in the following form:</p><disp-formula id="scirp.75814-formula402"><label>(35)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x89.png"  xlink:type="simple"/></disp-formula><p>To avoid the difficulty of fractional integration, for the exponential term we</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Approximate solution for Equation (21) is obtained by HPM with different values of a</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2310712x90.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Comparison of approximate solution by using HPM and VIM at<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x92.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2310712x91.png"/></fig><p>take the truncated taylor expansion in (35), e.g., <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x93.png" xlink:type="simple"/></inline-formula>to satisfy</p><p>the initial condition, we assume that the initial approximation has the following form <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x94.png" xlink:type="simple"/></inline-formula> first-order approximation takes the following form by using iteration Formula (35):</p><disp-formula id="scirp.75814-formula403"><graphic  xlink:href="http://html.scirp.org/file/5-2310712x95.png"  xlink:type="simple"/></disp-formula><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Approximate solution for Equation (21) at different values of a</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >x</th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x96.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x97.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x98.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x99.png" xlink:type="simple"/></inline-formula></th></tr></thead><tr><td align="center" valign="middle" >0.10</td><td align="center" valign="middle" >3.34976072</td><td align="center" valign="middle" >2.25905581</td><td align="center" valign="middle" >1.65729056</td><td align="center" valign="middle" >1.33243542</td></tr><tr><td align="center" valign="middle" >0.20</td><td align="center" valign="middle" >4.53980075</td><td align="center" valign="middle" >3.12399413</td><td align="center" valign="middle" >2.26453456</td><td align="center" valign="middle" >1.74033333</td></tr><tr><td align="center" valign="middle" >0.30</td><td align="center" valign="middle" >5.90653782</td><td align="center" valign="middle" >4.10952104</td><td align="center" valign="middle" >2.97106726</td><td align="center" valign="middle" >2.24160625</td></tr><tr><td align="center" valign="middle" >0.40</td><td align="center" valign="middle" >7.51861405</td><td align="center" valign="middle" >5.28303363</td><td align="center" valign="middle" >3.82016538</td><td align="center" valign="middle" >2.85706667</td></tr><tr><td align="center" valign="middle" >0.50</td><td align="center" valign="middle" >9.41603098</td><td align="center" valign="middle" >6.68897807</td><td align="center" valign="middle" >4.84789742</td><td align="center" valign="middle" >3.61067708</td></tr><tr><td align="center" valign="middle" >0.60</td><td align="center" valign="middle" >11.63265252</td><td align="center" valign="middle" >8.36684782</td><td align="center" valign="middle" >6.09010254</td><td align="center" valign="middle" >4.52980000</td></tr><tr><td align="center" valign="middle" >0.70</td><td align="center" valign="middle" >14.20081997</td><td align="center" valign="middle" >10.35549050</td><td align="center" valign="middle" >7.58435402</td><td align="center" valign="middle" >5.64544792</td></tr><tr><td align="center" valign="middle" >0.80</td><td align="center" valign="middle" >17.15275198</td><td align="center" valign="middle" >12.69459789</td><td align="center" valign="middle" >9.37077013</td><td align="center" valign="middle" >6.99253333</td></tr><tr><td align="center" valign="middle" >0.90</td><td align="center" valign="middle" >20.52106372</td><td align="center" valign="middle" >15.42533768</td><td align="center" valign="middle" >11.49244150</td><td align="center" valign="middle" >8.61011875</td></tr><tr><td align="center" valign="middle" >1.00</td><td align="center" valign="middle" >24.33898202</td><td align="center" valign="middle" >18.59065973</td><td align="center" valign="middle" >13.99570082</td><td align="center" valign="middle" >10.54166666</td></tr></tbody></table></table-wrap><disp-formula id="scirp.75814-formula404"><label>(36)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x100.png"  xlink:type="simple"/></disp-formula><p><xref ref-type="table" rid="table3">Table 3</xref> and <xref ref-type="fig" rid="fig4">Figure 4</xref> presents the approximate solution for the different values of a, we have noticed that the accuracy is improving. First, by computing more terms of the approximate solutions. The second way is taking more terms in the taylor expansion of the exponential term.</p><p>According to HPM, we build the following homotopy:</p><disp-formula id="scirp.75814-formula405"><label>(37)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x101.png"  xlink:type="simple"/></disp-formula><p>First, substitute the relation (16) in the Equation (37).</p><p>Second, equate the terms which have the same power’s of p which yield to the following series of Equations:</p><disp-formula id="scirp.75814-formula406"><label>(38)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x102.png"  xlink:type="simple"/></disp-formula><p>and so on, applying the operator <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x103.png" xlink:type="simple"/></inline-formula> to the previous Equations, and use the initial condition (34), to gain the following Equations:</p><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Approximate solution for Equation (34) is obtained by VIM with different values of a</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2310712x104.png"/></fig><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Approximate solution for Equation (34) at different values of α</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >x</th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x105.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x106.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x107.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x108.png" xlink:type="simple"/></inline-formula></th></tr></thead><tr><td align="center" valign="middle" >0.10</td><td align="center" valign="middle" >1.74557818</td><td align="center" valign="middle" >1.39713379</td><td align="center" valign="middle" >1.20769182</td><td align="center" valign="middle" >1.10549607</td></tr><tr><td align="center" valign="middle" >0.20</td><td align="center" valign="middle" >2.06685872</td><td align="center" valign="middle" >1.63644354</td><td align="center" valign="middle" >1.37989631</td><td align="center" valign="middle" >1.22394067</td></tr><tr><td align="center" valign="middle" >0.30</td><td align="center" valign="middle" >2.39681774</td><td align="center" valign="middle" >1.88439871</td><td align="center" valign="middle" >1.56340949</td><td align="center" valign="middle" >1.35821602</td></tr><tr><td align="center" valign="middle" >0.40</td><td align="center" valign="middle" >2.74634447</td><td align="center" valign="middle" >2.15443603</td><td align="center" valign="middle" >1.76686736</td><td align="center" valign="middle" >1.51115026</td></tr><tr><td align="center" valign="middle" >0.50</td><td align="center" valign="middle" >3.11743203</td><td align="center" valign="middle" >2.45118250</td><td align="center" valign="middle" >1.99492282</td><td align="center" valign="middle" >1.68553370</td></tr><tr><td align="center" valign="middle" >0.60</td><td align="center" valign="middle" >3.51027500</td><td align="center" valign="middle" >2.77680334</td><td align="center" valign="middle" >2.25082546</td><td align="center" valign="middle" >1.88413224</td></tr><tr><td align="center" valign="middle" >0.70</td><td align="center" valign="middle" >3.92459385</td><td align="center" valign="middle" >3.13252560</td><td align="center" valign="middle" >2.53716249</td><td align="center" valign="middle" >2.10969920</td></tr><tr><td align="center" valign="middle" >0.80</td><td align="center" valign="middle" >4.35998072</td><td align="center" valign="middle" >3.51915865</td><td align="center" valign="middle" >2.85615873</td><td align="center" valign="middle" >2.36498657</td></tr><tr><td align="center" valign="middle" >0.90</td><td align="center" valign="middle" >4.81599699</td><td align="center" valign="middle" >3.93730820</td><td align="center" valign="middle" >3.20982456</td><td align="center" valign="middle" >2.65275652</td></tr><tr><td align="center" valign="middle" >1.00</td><td align="center" valign="middle" >5.29219083</td><td align="center" valign="middle" >4.38747370</td><td align="center" valign="middle" >3.60003936</td><td align="center" valign="middle" >2.97579365</td></tr></tbody></table></table-wrap><disp-formula id="scirp.75814-formula407"><label>(39)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x109.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula408"><label>(40)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x110.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula409"><label>(41)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x111.png"  xlink:type="simple"/></disp-formula><p>and so on, by taking the truncated taylor expansions for the exponential term in</p><p>(40, 41): e.g.,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x112.png" xlink:type="simple"/></inline-formula>.</p><p>To avoid the difficulty of fractional integration, thus by solving Equations (39, 40, 41), we obtain <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x113.png" xlink:type="simple"/></inline-formula></p><disp-formula id="scirp.75814-formula410"><label>(42)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x114.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula411"><label>(43)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x115.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula412"><label>(44)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x116.png"  xlink:type="simple"/></disp-formula><p>the two terms approximation are formed as the following Equation</p><disp-formula id="scirp.75814-formula413"><label>(45)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x117.png"  xlink:type="simple"/></disp-formula><p><xref ref-type="table" rid="table4">Table 4</xref> and <xref ref-type="fig" rid="fig5">Figure 5</xref> shows the approximate solutions of (34) for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x118.png" xlink:type="simple"/></inline-formula> and for some values of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x118.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x119.png" xlink:type="simple"/></inline-formula>.</p><p><xref ref-type="fig" rid="fig6">Figure 6</xref> represent a comparison between two approximate solutions by using VIM and HPM methods.</p><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Approximate solution for Equation (34) is obtained by HPM with different values of a</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2310712x120.png"/></fig><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Comparison of approximate solution by using HPM and VIM at<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x122.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2310712x121.png"/></fig><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Approximate solution for Equation (34) at different values of a</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >x</th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x123.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x124.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x125.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x126.png" xlink:type="simple"/></inline-formula></th></tr></thead><tr><td align="center" valign="middle" >0.10</td><td align="center" valign="middle" >1.74531342</td><td align="center" valign="middle" >1.39709144</td><td align="center" valign="middle" >1.20768552</td><td align="center" valign="middle" >1.10549521</td></tr><tr><td align="center" valign="middle" >0.20</td><td align="center" valign="middle" >2.06407267</td><td align="center" valign="middle" >1.63581216</td><td align="center" valign="middle" >1.37976305</td><td align="center" valign="middle" >1.22391477</td></tr><tr><td align="center" valign="middle" >0.30</td><td align="center" valign="middle" >2.38608229</td><td align="center" valign="middle" >1.88141089</td><td align="center" valign="middle" >1.56263514</td><td align="center" valign="middle" >1.35803058</td></tr><tr><td align="center" valign="middle" >0.40</td><td align="center" valign="middle" >2.71887729</td><td align="center" valign="middle" >2.14558172</td><td align="center" valign="middle" >1.76420974</td><td align="center" valign="middle" >1.51041158</td></tr><tr><td align="center" valign="middle" >0.50</td><td align="center" valign="middle" >3.06116974</td><td align="center" valign="middle" >2.43083965</td><td align="center" valign="middle" >1.98807470</td><td align="center" valign="middle" >1.68339534</td></tr><tr><td align="center" valign="middle" >0.60</td><td align="center" valign="middle" >3.40994267</td><td align="center" valign="middle" >2.73693555</td><td align="center" valign="middle" >2.23607620</td><td align="center" valign="middle" >1.87906377</td></tr><tr><td align="center" valign="middle" >0.70</td><td align="center" valign="middle" >3.76162440</td><td align="center" valign="middle" >3.06236426</td><td align="center" valign="middle" >2.50903647</td><td align="center" valign="middle" >2.09921280</td></tr><tr><td align="center" valign="middle" >0.80</td><td align="center" valign="middle" >4.11218308</td><td align="center" valign="middle" >3.40478991</td><td align="center" valign="middle" >2.80699405</td><td align="center" valign="middle" >2.34530601</td></tr><tr><td align="center" valign="middle" >0.90</td><td align="center" valign="middle" >4.45689697</td><td align="center" valign="middle" >3.76108758</td><td align="center" valign="middle" >3.12924403</td><td align="center" valign="middle" >2.61840230</td></tr><tr><td align="center" valign="middle" >1.00</td><td align="center" valign="middle" >4.78998433</td><td align="center" valign="middle" >4.12719346</td><td align="center" valign="middle" >3.47425744</td><td align="center" valign="middle" >2.91904762</td></tr></tbody></table></table-wrap><p>Example 3.</p><p>Consider the following nonlinear (FQIDEs)</p><disp-formula id="scirp.75814-formula414"><label>(46)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x127.png"  xlink:type="simple"/></disp-formula><p>For <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x128.png" xlink:type="simple"/></inline-formula> with exact solution<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x128.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x129.png" xlink:type="simple"/></inline-formula>. By using VIM, the iteration formula for Equation (46) is,</p><disp-formula id="scirp.75814-formula415"><label>(47)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x130.png"  xlink:type="simple"/></disp-formula><p>We can take an initial approximation <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x131.png" xlink:type="simple"/></inline-formula></p><p>The first two iterations are easily obtained from (47) and are given by:</p><disp-formula id="scirp.75814-formula416"><label>(48)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x132.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula417"><label>(49)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x133.png"  xlink:type="simple"/></disp-formula><p>Therefore, we obtain the exact solution,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x134.png" xlink:type="simple"/></inline-formula>.</p><p>According to HPM, we construct the following homotopy:</p><disp-formula id="scirp.75814-formula418"><label>(50)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x135.png"  xlink:type="simple"/></disp-formula><p>and continuously trace an implicity defined curve from starting point <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x136.png" xlink:type="simple"/></inline-formula> to a solution function <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x136.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x137.png" xlink:type="simple"/></inline-formula> substituting (16) into (50). Also, we have to equate the terms with the same identical power’s of p, then, we gain these components</p><disp-formula id="scirp.75814-formula419"><label>(51)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x138.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula420"><label>(52)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x139.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula421"><label>(53)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x140.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75814-formula422"><label>(54)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2310712x141.png"  xlink:type="simple"/></disp-formula><p>and so on, we obtain <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x142.png" xlink:type="simple"/></inline-formula> therefore, the approximate solution is</p><p>obtained readily by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2310712x143.png" xlink:type="simple"/></inline-formula> which is the exact solution.</p></sec><sec id="s6"><title>6. Conclusion</title><p>In this paper, we have applied the VIM and HPM to find the solution of nonlinear initial value problem of fractional quadratic integro-differential equations for the first order. The methods do not require any linearization, perturbation or restrictive assumptions, we have observed that the VIM and HPM is a very powerful and effective tool for finding the solutions of the fractional quadratic integro-differential Equation. We use the Maple package (2015) in calculations.</p></sec><sec id="s7"><title>Cite this paper</title><p>Alhendi, F., Shammakh, W. and Al-Badrani, H. (2017) Numerical Solutions for Quadratic Integro- Differential Equations of Fractional Or- ders. Open Journal of Applied Sciences, 7, 157-170. https://doi.org/10.4236/ojapps.2017.74014</p></sec></body><back><ref-list><title>References</title><ref id="scirp.75814-ref1"><label>1</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>He</surname><given-names> J.H. </given-names></name>,<etal>et al</etal>. (<year>1999</year>)<article-title>Homotopy Perturbation Technique</article-title><source> Computer Methods in Applied Mechanics and Engineering</source><volume> 178</volume>,<fpage> 257</fpage>-<lpage>262</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.75814-ref2"><label>2</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>He</surname><given-names> J.H. </given-names></name>,<etal>et al</etal>. (<year>1999</year>)<article-title>Variational Iteration Method—A Kind of Non-Linear Analytical Technique: Some Examples</article-title><source> International Journal of Non-Linear Mechanics</source><volume> 34</volume>,<fpage> 699</fpage>-<lpage>708</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.75814-ref3"><label>3</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>He</surname><given-names> J.H. </given-names></name>,<etal>et al</etal>. (<year>2000</year>)<article-title>A Coupling Method of a Homotopy Technique and a Perturbation Technique for Non-Linear Problems</article-title><source> International Journal of Non-Linear Mechanics</source><volume> 35</volume>,<fpage> 37</fpage>-<lpage>43</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.75814-ref4"><label>4</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>He</surname><given-names> J.H. </given-names></name>,<etal>et al</etal>. (<year>2003</year>)<article-title>Homotopy Perturbation Method: A New Nonlinear Analytical Technique</article-title><source> Applied Mathematics and Computation</source><volume> 135</volume>,<fpage> 73</fpage>-<lpage>79</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.75814-ref5"><label>5</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Abbasbandy</surname><given-names> S. </given-names></name>,<etal>et al</etal>. (<year>2007</year>)<article-title>An Approximation Solution of a Nonlinear Equation with Riemann-Liouville’s Fractional Derivatives by He’s Variational Iteration Method</article-title><source> Journal of Computational and Applied Mathematics</source><volume> 207</volume>,<fpage> 53</fpage>-<lpage>58</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.75814-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Daftardar-Gejji, V. and Jafari, H. (2007) Solving a Multi-Order Fractional Differential Equation Using Adomian Decomposition. Applied Mathematics and Computation, 189, 541-548.</mixed-citation></ref><ref id="scirp.75814-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Odibat, Z. and Momani, S. (2009) The Variational Iteration Method: An Efficient Scheme for Handling Fractional Partial Differential Equations in Fluid Mechanics. Computers &amp; Mathematics with Applications, 58, 2199-2208.</mixed-citation></ref><ref id="scirp.75814-ref8"><label>8</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Sweilam</surname><given-names> N.H. </given-names></name>,<etal>et al</etal>. (<year>2007</year>)<article-title>Fourth Order Integro-Differential Equations Using Variational Iteration Method</article-title><source> Computers &amp; Mathematics with Applications</source><volume> 54</volume>,<fpage> 1086</fpage>-<lpage> 1091</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.75814-ref9"><label>9</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Yildirim</surname><given-names> A. </given-names></name>,<etal>et al</etal>. (<year>2008</year>)<article-title>Solution of BVPs for Fourth-Order Integro-Differential Equations by Using Homotopy Perturbation Method</article-title><source> Computers &amp; Mathematics with Applications</source><volume> 56</volume>,<fpage> 3175</fpage>-<lpage>3180</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.75814-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Benson, D.A., Wheatcraft, S.W. and Meerschaert, M.M. (2000) Application of a Fractional Advection-Dispersion Equation. Water Resources Research, 36, 1403- 1412.  
https://doi.org/10.1029/2000WR900031</mixed-citation></ref><ref id="scirp.75814-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Miller, K.S. and Ross, B. (1993) An Introduction to the Fractional Calculus and Fractional Differential Equations. Wiley-Interscience, Hoboken.</mixed-citation></ref><ref id="scirp.75814-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Podlubny, I. (1998) Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications. Academic Press, Cambridge, MA.</mixed-citation></ref><ref id="scirp.75814-ref13"><label>13</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Podlubny</surname><given-names> I. </given-names></name>,<etal>et al</etal>. (<year>2002</year>)<article-title>Correction to Figure 4 in Geometric and Physical Interpretation of Fractional Integration and Fractional Differentiation</article-title><source> Fractional Calculus &amp; Applied Analysis</source><volume> 5</volume>,<fpage> 367</fpage>-<lpage>386</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.75814-ref14"><label>14</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Nawaz</surname><given-names> Y. </given-names></name>,<etal>et al</etal>. (<year>2011</year>)<article-title>Variational Iteration Method and Homotopy Perturbation Method for Fourth-Order Fractional Integro-Differential Equations</article-title><source> Computers &amp; Mathematics with Applications</source><volume> 61</volume>,<fpage> 2330</fpage>-<lpage>2341</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.75814-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Momani, S. and Noor, M.A. (2006) Numerical Methods for Fourth-Order Fractional Integro-Differential Equations. Applied Mathematics and Computation, 182, 754-760.</mixed-citation></ref><ref id="scirp.75814-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Elbeleze, A.A., Kilicman, A. and Taib, B.M. (2012) Application of Homotopy Perturbation and Variational Iteration Methods for Fredholm Integrodifferential Equation of Fractional Order. Abstract and Applied Analysis, 2012, Article ID: 763139.</mixed-citation></ref><ref id="scirp.75814-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Elbeleze, A.A., Kilicman, A. and Taib, B.M. (2016) Approximate Solution of Integro-Differential Equation of Fractional (Arbitrary) Order. Journal of King Saud University-Science, 28, 61-68.</mixed-citation></ref><ref id="scirp.75814-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Kadem, A. and Kilicman, A. (2012) The Approximate Solution of Fractional Fredholm Integrodifferential Equations by Variational Iteration and Homotopy Perturbation Methods. Abstract and Applied Analysis, 2012, Article ID: 486193.  
https://doi.org/10.1155/2012/486193</mixed-citation></ref><ref id="scirp.75814-ref19"><label>19</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Gaafar</surname><given-names> F.M. </given-names></name>,<etal>et al</etal>. (<year>2014</year>)<article-title>Positive Solutions of a Quadratic Integro-Differential Equation</article-title><source> Journal of the Egyptian Mathematical Society</source><volume> 22</volume>,<fpage> 162</fpage>-<lpage>166</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref></ref-list></back></article>