<?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">AM</journal-id><journal-title-group><journal-title>Applied Mathematics</journal-title></journal-title-group><issn pub-type="epub">2152-7385</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/am.2012.35061</article-id><article-id pub-id-type="publisher-id">AM-19048</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Variational Iterative Method Applied to Variational Problems with Moving Boundaries
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>ateme</surname><given-names>Ghomanjani</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Sara</surname><given-names>Ghaderi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Applied Mathematics, Ferdowsi University of Mashhad, Mashhad, Iran</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>fatemeghomanjani@gmail.com(AG)</email>;<email>s_gh333@yahoo.com(SG)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>16</day><month>05</month><year>2012</year></pub-date><volume>03</volume><issue>05</issue><fpage>395</fpage><lpage>402</lpage><history><date date-type="received"><day>January</day>	<month>3,</month>	<year>2012</year></date><date date-type="rev-recd"><day>March</day>	<month>19,</month>	<year>2012</year>	</date><date date-type="accepted"><day>March</day>	<month>26,</month>	<year>2012</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 paper, He’s variational iterative method has been applied to give exact solution of the Euler Lagrange equation which arises from the variational problems with moving boundaries and isoperimetric problems. In this method, general Lagrange multipliers are introduced to construct correction functional for the variational problems. The initial approximations can be freely chosen with possible unknown constant, which can be determined by imposing the boundary conditions. Illustrative examples have been presented to demonstrate the efficiency and applicability of the variational iterative method.
 
</p></abstract><kwd-group><kwd>Variational Iterative Method; Variational Problems; Moving Boundaries; Isoperimetric Problems</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>In modeling a large class of problems arising in science, engineering and economics, it is necessary to minimize amounts of a certain functional. Because of the important role of this subject, special attention has been given to these problems. Such problems are called variational problems, see [1,2].</p><p>The simplest form of a variational problem can be considered as</p><disp-formula id="scirp.19048-formula9969"><label>(1)</label><graphic position="anchor" xlink:href="1-7400723\f6b66c41-ecda-41a4-bf30-dfca1a986dd4.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="1-7400723\b0b0ac21-2f6d-4073-8e03-fed3211dcb13.jpg" /> is the functional which its extremum must be found. Functional <img src="1-7400723\3ceb7912-8838-4317-bca5-0e60c38d50c3.jpg" /> can be considered by two kinds of boundary conditions. In the fixed boundary problems, the admissible function <img src="1-7400723\e8bfbc4a-3dc4-41df-804b-b2a39370d004.jpg" /> must satisfy following boundary conditions</p><disp-formula id="scirp.19048-formula9970"><label>(2)</label><graphic position="anchor" xlink:href="1-7400723\f66664aa-d5ca-44a6-9172-b40512eb3067.jpg"  xlink:type="simple"/></disp-formula><p>In moving boundary problems, at least one of the boundary points of the admissible function is movable along a boundary curve. Further more many applications of the calculus of variations lead to problems in which not only boundary conditions, but also a quite different type of conditions known as constraints, are imposed on the admissible function. The necessary condition for the admissible solutions of such problems has to satisfy the Euler-Lagrange equation which is generally nonlinear.</p><p>In this work we consider He’s variational iterative method as a well known method for finding both analytic and approximate solutions of differential equations. Here, the problem is initially approximated with possible unknowns. Then a correction functional is constructed by a general Lagrange multiplier, which can be identified optimally via the variational theory [<xref ref-type="bibr" rid="scirp.19048-ref3">3</xref>].</p><p>Variational iterative method is applied on various kinds of problems [4-31].</p><p>Author of [<xref ref-type="bibr" rid="scirp.19048-ref32">32</xref>] solved variational problems with moving boundaries with Adomian decomposition method. Variational iterative method was applied to solve variational problems with fixed boundaries (see [11,27,30]). In this work we obtain exact solution of variational problems with moving boundaries and isoperimetric problems by variational iterative method. In fact, variational iterative method is applied to solve the Euler-Lagrange equation with prescribed boundary conditions. To present a clear overview of the procedure several illustrative examples are included.</p></sec><sec id="s2"><title>2. Variational Iterative Method</title><p>In variational iterative method which is stated by He [<xref ref-type="bibr" rid="scirp.19048-ref3">3</xref>], solutions of the problems are approximated by a set of functions that may include possible constants to be determined from the boundary conditions. In this method the problem is considered as</p><disp-formula id="scirp.19048-formula9971"><label>, (3)</label><graphic position="anchor" xlink:href="1-7400723\4779ac84-215d-4d68-b4ff-7bac34504be4.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="1-7400723\e7fb85a0-0e34-44c3-800c-6760fe534a82.jpg" /> is a linear operator, and <img src="1-7400723\c97aa89c-0ed7-4308-9617-9da0db51ca59.jpg" /> is a nonlinear operator. <img src="1-7400723\6f2eeef4-2a27-4169-a1cc-f97cd3f56c1e.jpg" />is an inhomogeneos term. By using the variational iterative method, the following correct functional is taken into account</p><disp-formula id="scirp.19048-formula9972"><label>(4)</label><graphic position="anchor" xlink:href="1-7400723\68e73155-d212-4e14-9ec4-23066d0326a4.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="1-7400723\53e952ad-8912-4966-8b8a-81648ce571ad.jpg" /> is Lagrange multiplier [<xref ref-type="bibr" rid="scirp.19048-ref5">5</xref>], the subscript n denotes the n-th approximation, <img src="1-7400723\3eb3a6c5-816f-408f-aeb2-e78bf91a7c0c.jpg" />is as a restricted variation i.e. <img src="1-7400723\151e1d64-69f4-4ad9-8bd5-240f05c5b11c.jpg" />[6-8]. Taking the variation from both sides of the correct functional with respect to y<sub>n</sub> and imposing<img src="1-7400723\cb9f0f09-4a65-40c3-886d-5f860bf3fc66.jpg" />, the stationary conditions are obtained. By using the stationary conditions the optimal value of the <img src="1-7400723\04a0cbf2-c7aa-4fe3-85ec-344a1e055d3f.jpg" /> can be identified.</p><p>The successive approximation <img src="1-7400723\a4d43f9a-4586-464e-91a9-b6d134a3162c.jpg" /> can be established by determining a general lagrangian multiplier <img src="1-7400723\13570fda-4049-4fd7-8acb-03a602ab6041.jpg" /> and initial solution<img src="1-7400723\d96c8a72-e59c-43ad-9df1-f0c35589a2c3.jpg" />. Since this procedure avoids the discretization of the problem, it is possible to find the closed form solution without any round off error.</p><p>In the case of m equations, the equations are rewritten in the form of:</p><disp-formula id="scirp.19048-formula9973"><label>(5)</label><graphic position="anchor" xlink:href="1-7400723\f5ec8b48-ee7f-46da-868f-136af4f37fc0.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="1-7400723\536c313b-98a0-4d24-85d7-fc8d5dfeac71.jpg" /> is a linear with respect to<img src="1-7400723\1fcc5de2-0ee3-4bb7-afcf-1c0ddb3e6670.jpg" />, and <img src="1-7400723\5393b6be-3c6b-459f-b279-a2d2b15c713e.jpg" /> is nonlinear part of the ith equation. In this case the correct functionals are produced as</p><disp-formula id="scirp.19048-formula9974"><label>(6)</label><graphic position="anchor" xlink:href="1-7400723\9154b294-fed3-4ad0-b325-8953609f55ed.jpg"  xlink:type="simple"/></disp-formula><p>and the optimal values of the <img src="1-7400723\dfafd53c-3096-4181-8463-862798f6862e.jpg" /> are obtained by taking the variation from both sides of the correct functionals and finding stationary conditions using<img src="1-7400723\52caa12a-241e-462b-8b55-7294416a7eab.jpg" />。</p></sec><sec id="s3"><title>3. Statement of the Problem</title><sec id="s3_1"><title>3.1. Moving Boundary Problems</title><p>The necessary condition for the solution of problem (1) is to satisfy the Euler-Lagrange equation</p><disp-formula id="scirp.19048-formula9975"><label>(7)</label><graphic position="anchor" xlink:href="1-7400723\dde27fcc-9f0c-4757-84b4-3746291aa3b1.jpg"  xlink:type="simple"/></disp-formula><p>The general form of the variational problem (1) is</p><disp-formula id="scirp.19048-formula9976"><label>(8)</label><graphic position="anchor" xlink:href="1-7400723\d4f43d19-9d6f-4240-badb-bb5aa7f6ec7b.jpg"  xlink:type="simple"/></disp-formula><p>Here, the necessary condition for the extermum of the functional (8) is to satisfy the following system of second-order differential equations</p><disp-formula id="scirp.19048-formula9977"><label>(9)</label><graphic position="anchor" xlink:href="1-7400723\474ac485-9e87-49f9-870d-664f284907fd.jpg"  xlink:type="simple"/></disp-formula><p>In fixed boundary problems, Euler-Lagrange equation must be considered by the boundary conditions, but for the problems with variable boundaries, Euler-Lagrange equation must satisfy natural boundary conditions or transversality conditions which will be described in the following theorems.</p><p>For the problems with variable boundaries, we have two cases:</p><p>Type 1: As the first case, those problems are considered for which at least one of the boundary points move freely along a line parallel to the y-axis. Indeed at this point <img src="1-7400723\02a91a76-f10f-4874-a898-90d7a26725cb.jpg" /> is not specified. In this case all admissible functions have the same domain <img src="1-7400723\c8029b34-f88e-4ab8-8e27-0fa872ddb43f.jpg" /> and satisfy the Euler-Lagrange equation in this interval. Furthermore such functions have to satisfy conditions called natural boundary conditions stated in the following theorem.</p><p>Theorem 3.1. Suppose the function <img src="1-7400723\378d0dec-fc05-4454-82cd-544d43f6fbb2.jpg" /> in<img src="1-7400723\8a09acc0-1080-4e5f-8b92-fb5f9a9ff2a2.jpg" />, yields a relative minimum of the functional (1) that for which<img src="1-7400723\9a33be25-03fc-428e-be00-15e3c263f1d3.jpg" />, <img src="1-7400723\729d4917-c8e9-48e0-8b7b-07e3b560ef2c.jpg" />is arbitrary (free right endpoint) and <img src="1-7400723\0513d845-cb7d-4262-ba13-a61d462d9df9.jpg" /> are arbitrary (free endpoints). Then <img src="1-7400723\18d3cb6c-c2f2-413c-a1ad-a1cfb3fad550.jpg" /> satisfies, the following natural boundary conditions, respectively:</p><disp-formula id="scirp.19048-formula9978"><label>(10)</label><graphic position="anchor" xlink:href="1-7400723\e73425bb-da6c-46d9-a087-f383cb0dfdd1.jpg"  xlink:type="simple"/></disp-formula><p>or</p><disp-formula id="scirp.19048-formula9979"><label>(11)</label><graphic position="anchor" xlink:href="1-7400723\0d1b515e-0b4a-4f3d-a82d-52081c1dbd4b.jpg"  xlink:type="simple"/></disp-formula><p>Type 2: For the second case, the beginning and end points (or only one of them) can move freely on given curves<img src="1-7400723\80e524ca-c2c5-4520-8181-d6b568be10a9.jpg" />. In this case, a function <img src="1-7400723\b41064f3-32f9-4a22-bf7b-e52ec3c3a5f1.jpg" /><img src="1-7400723\00966c41-3ad7-4947-afd1-5943c1d074b0.jpg" /> is required, which emanates at some <img src="1-7400723\253d1f5b-33ae-4bf1-a4fe-c44a1e9c7c8f.jpg" /> from the curve <img src="1-7400723\6db03db8-211f-4550-8200-3e684070969e.jpg" /> and terminates for some <img src="1-7400723\05c85c60-5f22-43ed-94b0-2359171f9a27.jpg" /> on the curve <img src="1-7400723\e8fdccbc-a08d-40d7-a262-d0697e3a5ed5.jpg" /> and minimizes the functional (1). In this problem, the points <img src="1-7400723\6efd3f07-b584-4408-826a-df5a0574c4ef.jpg" /> are not known, and they must satisfy the necessary conditions called transversality conditions, described in the following theorem.</p><p>Theorem 3.2. If the function<img src="1-7400723\0eeaf1c1-d327-4dd5-9f5c-f9f192da005a.jpg" /><img src="1-7400723\bc182cb9-130a-4d9d-a696-96a37ec7480d.jpg" />, which emanates at some<img src="1-7400723\0b26c9cb-3061-4baa-a0b8-0305464c3308.jpg" /> from the curve <img src="1-7400723\78f677d5-da76-4e40-b4c3-fd6e67f08e7e.jpg" /> and terminates for some <img src="1-7400723\68df3b6c-cac8-4ff3-8faa-0cdd9e5a73d3.jpg" /> on the curve <img src="1-7400723\f14c9873-f8aa-46b8-8f6a-ee10d9c7c9e6.jpg" />, yields a relative minimum for functional (1), where<img src="1-7400723\8ae3503c-9933-47dc-a189-0d30b05b1350.jpg" /><img src="1-7400723\ad77216f-bb44-4c35-bf88-d5a3c5e66a45.jpg" />, R being a domain in the <img src="1-7400723\5eff48e9-1177-43ba-a70b-14dab11c4401.jpg" /> space that contains all lineal elements of<img src="1-7400723\b0d7aee5-882c-4512-ae55-3f64ca708f61.jpg" />, then it is necessary that <img src="1-7400723\4c1e83bf-57c7-4419-8a68-36a9b649ac21.jpg" /><img src="1-7400723\4f8a078e-3b36-4e18-ba61-47705c283282.jpg" /> to satisfy the Euler-Lagrange equation in the interval <img src="1-7400723\a091a822-396c-4e54-a396-811ded23b996.jpg" /> and at the point of exit and the point of entrance, the following transversality conditions to be satisfied:</p><disp-formula id="scirp.19048-formula9980"><label>(12)</label><graphic position="anchor" xlink:href="1-7400723\77763ef2-8518-4228-9df8-6da0714e60ff.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.19048-formula9981"><label>(13)</label><graphic position="anchor" xlink:href="1-7400723\2faa892d-6a3b-4bff-a834-a6c837585b85.jpg"  xlink:type="simple"/></disp-formula><p>In the case that one of the points is fixed, then the transversality condition has to be held at the other point.</p><p>One can consider transversality conditions for the problems with more than one unknown functions. For example, in to minimize two dimensional case, a vector function <img src="1-7400723\7ceb6402-e329-4ab7-9b22-1b90de9eb078.jpg" /> is looked for such that</p><disp-formula id="scirp.19048-formula9982"><label>(14)</label><graphic position="anchor" xlink:href="1-7400723\4d377d3b-c9da-40a7-b0ee-1949025966d8.jpg"  xlink:type="simple"/></disp-formula><p>in which <img src="1-7400723\d9f7faa9-ad48-43b6-8b15-32ca4a798bd2.jpg" /> and the endpoint lies on a two-dimensional surface that is given by<img src="1-7400723\78511804-ef48-4eeb-916d-19a41e026bc7.jpg" />. Here the transversality conditions at <img src="1-7400723\35edac0f-5d0c-4a91-896b-7ccd15ee66f1.jpg" /> are:</p><disp-formula id="scirp.19048-formula9983"><label>(15)</label><graphic position="anchor" xlink:href="1-7400723\b14c12d1-50f5-4578-8aa9-b7b2ea211a5a.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.19048-formula9984"><label>(16)</label><graphic position="anchor" xlink:href="1-7400723\fe56fd9e-0747-4a2d-b2ad-8cb336d22e48.jpg"  xlink:type="simple"/></disp-formula><p>In which <img src="1-7400723\a928eac1-4ac4-4ec7-9a30-d627a3da6abf.jpg" /> is an admissible vector function.</p><p>For further information on transversality conditions, specially for the proofs of Theorems 3.1 and 3.2 and conditions (15), (16), see [<xref ref-type="bibr" rid="scirp.19048-ref2">2</xref>].</p><p>Example 3.1. Consider the following functional:</p><disp-formula id="scirp.19048-formula9985"><label>(17)</label><graphic position="anchor" xlink:href="1-7400723\8814a84d-b6f4-499e-affc-9cc2eb26eead.jpg"  xlink:type="simple"/></disp-formula><p>In which <img src="1-7400723\c4fafb32-1aed-4836-a713-20091f25a571.jpg" /> and <img src="1-7400723\bb6c2c41-945f-4a30-ac99-d8082c5222c5.jpg" /> is the amount of a capital at time t (see [<xref ref-type="bibr" rid="scirp.19048-ref1">1</xref>]).</p><p>Here, the capital stock <img src="1-7400723\50a69b60-ec21-4079-a059-a54929c2cc02.jpg" /> at the initial time <img src="1-7400723\485741b0-6bd8-4107-9163-8ed4a39f87da.jpg" /> of the planning period is assumed to be known:<img src="1-7400723\7d8a85f7-0608-4791-98fe-dfb625620137.jpg" />; on the other hand, the planner won’t wish to explain how large the capital would be at time<img src="1-7400723\83922c28-cfbc-4f7d-a6d6-70e45ab7b673.jpg" />. Therefore, there is a variational problem with free right endpoint. Here we let<img src="1-7400723\51b7aa33-5512-45e9-bee8-bc6ee4c32ced.jpg" />, and <img src="1-7400723\96c6a211-07a6-4430-b7c4-7da22fc1bd6a.jpg" /> which has the analytical solution <img src="1-7400723\a18024f1-2551-43a8-8fe4-1005ff119698.jpg" /> The corresponding Euler-Lagrange equation is:</p><p><img src="1-7400723\0d28eb26-d424-4b69-b86d-973ad4d923f5.jpg" /></p><p>Now natural boundary condition at <img src="1-7400723\a823cb22-f873-4d3d-b194-3601fc5386f0.jpg" /> is as following:</p><p><img src="1-7400723\00869b44-650f-44fa-b426-90dac997dc4b.jpg" /></p><p>Therefore, the following boundary conditions are:</p><disp-formula id="scirp.19048-formula9986"><label>(18)</label><graphic position="anchor" xlink:href="1-7400723\7583e453-a365-494f-8e71-04c1585008cc.jpg"  xlink:type="simple"/></disp-formula><p>By using variational iterative method we consider the following functional is considered:</p><p><img src="1-7400723\bd5b6d43-0269-4d31-906d-dd13f452e290.jpg" /></p><p>Taking the variation from both sides of the correct functional with respect to <img src="1-7400723\185af4e8-b090-4a02-a2ec-c39c64f4dafd.jpg" /> given:</p><p><img src="1-7400723\5c8f0d26-ce1c-4974-8224-1dc0a143c3f7.jpg" /></p><p>For all variations <img src="1-7400723\d14c686e-2e63-4c73-b5c9-93c25a2386cb.jpg" /> and<img src="1-7400723\ee543810-19ab-4c03-93e5-d58e293d7874.jpg" />. The following stationary conditions are obtained:</p><p><img src="1-7400723\69c33d36-442f-451c-8806-972dd0ca793a.jpg" />,</p><p><img src="1-7400723\601981aa-6cf0-4ef8-b2f9-088949a4863a.jpg" /></p><p><img src="1-7400723\f589ea0d-1fe3-47ec-80d4-8de220bccf9b.jpg" /></p><p>So that<img src="1-7400723\e477d2dc-9d97-4c3c-89ce-4dc09622c236.jpg" />. Therefore iterative formula can be found as:</p><p><img src="1-7400723\5bc97ccf-fa54-41d1-804a-417a27187796.jpg" /></p><p>If <img src="1-7400723\c509c9ca-db48-4026-9884-de1bf315d992.jpg" /> then</p><p><img src="1-7400723\f66536f1-e0f4-46e9-99d3-f01e81ecd688.jpg" /></p><p>By imposing (18) <img src="1-7400723\d96556fb-6b01-4044-82fb-8c3cc1c7581a.jpg" />are resulted. Which yields the exact solutions of the problem (see <xref ref-type="fig" rid="fig1">Figure 1</xref>).</p><p>Example 3.2. We want to find the shortest distance from the point <img src="1-7400723\e3d5af4e-54b1-46e9-bf82-7d5aeffd3c63.jpg" /> to the sphere</p><p><img src="1-7400723\005c10fa-0344-4f9c-be6d-36c97a45695f.jpg" /></p><p>This problem is reduced to optimize the following functional:</p><disp-formula id="scirp.19048-formula9987"><label>(19)</label><graphic position="anchor" xlink:href="1-7400723\9679efce-90c9-4974-b27a-65b2c2b52710.jpg"  xlink:type="simple"/></disp-formula><p>where the point <img src="1-7400723\0e1cec74-b8ec-47b9-93d1-42b4068a1ac1.jpg" /> must lie on the sphere, with the exact solution<img src="1-7400723\229021af-e3f0-4540-bdd4-35f7fea9a76b.jpg" />, see [<xref ref-type="bibr" rid="scirp.19048-ref33">33</xref>]. The corresponding Euler Lagrange equations for this problem</p><p>are:</p><p><img src="1-7400723\eb55c780-61e0-46d6-8de9-0e766193da9b.jpg" /></p><p>So that</p><p><img src="1-7400723\91716aa1-c89b-4733-93fa-49683b4c6bbd.jpg" /></p><p>In above equations “e” and “f” are constant, so they can be rewritten as:</p><p><img src="1-7400723\b73a7e5b-d1e7-463f-93ea-4deab706bd2f.jpg" /></p><p>The transversality conditions are:</p><disp-formula id="scirp.19048-formula9988"><label>(20)</label><graphic position="anchor" xlink:href="1-7400723\b9acbe5e-1535-42b9-85d0-2fe13029f97e.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.19048-formula9989"><label>(21)</label><graphic position="anchor" xlink:href="1-7400723\225b999e-e0a7-43f4-bb83-3b5e4821595d.jpg"  xlink:type="simple"/></disp-formula><p>By using variational iteration method results:</p><p><img src="1-7400723\acb7d9c7-83bb-4dab-8cde-c348ca41dc17.jpg" /></p><p>and</p><p><img src="1-7400723\7ded91d9-1396-4822-8769-17a5bc6bb810.jpg" /></p><p>The variation from both sides of above equations for finding the optimal value of <img src="1-7400723\a207d701-4ea9-4e5a-8084-7de77f2199b9.jpg" /> is:</p><p><img src="1-7400723\9c6f5596-dbee-47d9-a31d-b5761d61d4c4.jpg" /></p><p>and</p><p><img src="1-7400723\f65fc2c8-7f1a-42aa-80a7-6e969961f15c.jpg" /></p><p>Therefore</p><p><img src="1-7400723\df546e00-fd9a-4997-9c25-5eec47993058.jpg" />.</p><p>and</p><p><img src="1-7400723\4191d51c-b742-4cbe-9d4a-cbd0206121c0.jpg" /></p><p>which yields:</p><p><img src="1-7400723\d4df27f5-df7e-48b7-bced-7905c0c67c1d.jpg" /></p><p>So that the following iterative formulas are obtained:</p><p><img src="1-7400723\54597ce9-2681-4f15-ab7b-1f57076ef616.jpg" /></p><p><img src="1-7400723\fc4dc9f1-0291-4f1c-9a4c-6c1f85de1da3.jpg" /></p><p>If <img src="1-7400723\d9078ed0-eb25-427c-b8e9-801300f8f56c.jpg" /> then we have:</p><p><img src="1-7400723\0ac06fc1-e206-43a8-9b13-2d65e08a33ce.jpg" /></p><p>and</p><p><img src="1-7400723\e1e1ba45-ea83-4a13-bf87-f55490d860eb.jpg" /></p><p>By choosing<img src="1-7400723\3ba18105-fd5a-459e-8f60-44ee104fc60d.jpg" />,</p><p><img src="1-7400723\cf680bc4-a060-45d2-9e84-089faa17eba4.jpg" /></p><p>Imposing (20) and (21) lead to, <img src="1-7400723\f5326649-24ae-4848-b36c-214bc9a70454.jpg" />therefore:</p><p><img src="1-7400723\8dae9095-b90b-490f-9b98-74ca00d66a3d.jpg" /></p><p>which is the exact solution.</p></sec><sec id="s3_2"><title>3.2. Isoperimetric Problems</title><p>Assume that two functions <img src="1-7400723\b43aea4c-d378-4c46-a269-1807459a65e9.jpg" /> and <img src="1-7400723\4f18c283-5311-464d-9112-63ad4ce9ad15.jpg" /> are given. Among all curves <img src="1-7400723\4b76da02-1428-43f0-aad6-6fcc94998094.jpg" /> along which the functional</p><p><img src="1-7400723\dea00be4-62bc-4935-b50b-8d0db55f4a71.jpg" /></p><p>assumes a given value l, determine the one for which the functional</p><p><img src="1-7400723\e6644660-a075-4e45-8f5e-d874e52526c1.jpg" /></p><p>Gives an extermal value. Suppose that F and G have continuous first and second partial derivatives for <img src="1-7400723\78a35081-704f-46b0-a8b3-57c4e76b624d.jpg" /> and for arbitrary values of the variables <img src="1-7400723\a392ec93-7b9a-4231-a4b1-db562bc16ce7.jpg" /> and<img src="1-7400723\fd07eba4-2897-4a22-9c57-00c781a3cc5f.jpg" />.</p><p>Euler’s theorem: If a curve <img src="1-7400723\b7aa525d-35cf-4f57-92a7-3a7e1db145f4.jpg" /> extremizes the functional J <img src="1-7400723\9a45019d-2d4b-4d13-b58f-a39664164a4d.jpg" /> under the conditions</p><p><img src="1-7400723\2404f140-c00d-4cc4-a238-759273eede1b.jpg" /></p><p><img src="1-7400723\7194f084-93a9-4bbb-b3a2-a2f02f045ec7.jpg" /></p><p>and if <img src="1-7400723\87a17979-6f25-4e54-abf3-49b93bf15297.jpg" /> is not an extremal of the functional K, there exists a constant <img src="1-7400723\12ca2d2e-8d8f-40a9-8f63-2fc191f70166.jpg" /> such that the curve <img src="1-7400723\c89795a0-09a3-4518-b253-132fccef8960.jpg" /> is an extremal of the functional</p><p><img src="1-7400723\0a46d0b8-a69a-418b-b103-c843aedf9284.jpg" /></p><p>The necessary condition for the solution of this problem is to satisfy the Euler-Lagrange equation</p><p><img src="1-7400723\72f2af9b-94fc-4d18-a968-9d7106e7305e.jpg" /></p><p>with given boundary conditions in which <img src="1-7400723\74b85a73-3dd0-4914-a01f-a02b5b6d40c6.jpg" /> for further information (see [<xref ref-type="bibr" rid="scirp.19048-ref2">2</xref>]).</p><p>Example 3.3. It is aimed to find the minimum of the functional</p><disp-formula id="scirp.19048-formula9990"><label>(22)</label><graphic position="anchor" xlink:href="1-7400723\fa37e50d-c55c-4375-ace2-bec253695cd5.jpg"  xlink:type="simple"/></disp-formula><p>Such that</p><disp-formula id="scirp.19048-formula9991"><label>(23)</label><graphic position="anchor" xlink:href="1-7400723\aaca0c6b-b2cb-49fa-9376-bc269a355a17.jpg"  xlink:type="simple"/></disp-formula><p>and</p><disp-formula id="scirp.19048-formula9992"><label>(24)</label><graphic position="anchor" xlink:href="1-7400723\38fedd1a-135c-4f25-9d31-33e52a9e5dd0.jpg"  xlink:type="simple"/></disp-formula><p>With exact solution <img src="1-7400723\1334aee2-3fc6-4afe-822b-f013bb5029d9.jpg" /> [<xref ref-type="bibr" rid="scirp.19048-ref19">19</xref>]. According to the following auxiliary functional:</p><p><img src="1-7400723\7e99e25a-db12-4b2c-b838-2052d4619aba.jpg" /></p><p>and the corresponding Euler-Lagrange equation:</p><p><img src="1-7400723\48e50511-247e-4c24-91cf-728ea170b364.jpg" /></p><p>so</p><p><img src="1-7400723\4e3b59ce-704c-461e-ae7c-f058ebb1b491.jpg" /></p><p>By applying He’s variational iterative method results</p><p><img src="1-7400723\ea7e9ba9-b162-482c-ba51-3ee315aab1c3.jpg" /></p><p>To find the optimal value of <img src="1-7400723\f3de0d7a-8996-4491-96c0-b4cb49c6b71f.jpg" /> following equation is required:</p><p><img src="1-7400723\d0cc9b00-0307-4d39-ac07-f43772d5acf3.jpg" /></p><p>Therefore, the stationary conditions are obtained in the following form:</p><p><img src="1-7400723\76e941a4-580f-49fb-9044-5c12f3fedee0.jpg" /></p><p>which yields</p><p><img src="1-7400723\81c6047d-3b2d-4b85-867b-8acccc549315.jpg" /></p><p>and the desired sequence is</p><p><img src="1-7400723\d4420363-230d-4826-a73e-4b4f1c18ffb8.jpg" /></p><p>By choosing <img src="1-7400723\968dc480-8fca-4b05-8f84-19071875fe2d.jpg" /></p><p><img src="1-7400723\36120049-fc1d-4b91-b520-ca7215c0c86a.jpg" /></p><p>Imposing (24) on this function given</p><disp-formula id="scirp.19048-formula9993"><label>(25)</label><graphic position="anchor" xlink:href="1-7400723\5ba36576-c240-438b-953a-697ac511a8ed.jpg"  xlink:type="simple"/></disp-formula><p>If <img src="1-7400723\7a80b2a0-d5b9-47ae-9000-e42602ce6294.jpg" /> then from (24) <img src="1-7400723\f822f68d-7a55-435a-97e2-d2805d04d3f9.jpg" />but from (23) <img src="1-7400723\fa210d48-ca43-432e-a3c5-75b2fbead702.jpg" /> which is a contradiction.</p><p>Now imposing (24), we have: <img src="1-7400723\d963ebf8-18c6-46ea-842e-d469f9ba88bc.jpg" /></p><p>so <img src="1-7400723\fa5dc99d-e6f4-4a46-b3f8-90deb7a7cd4d.jpg" /> and it is known that in this case imposing (24) on the Euler Lagrange equation yields</p><p><img src="1-7400723\0f22ab37-87ac-49a4-b386-47125db0ceb0.jpg" /></p><p>Hence:</p><p><img src="1-7400723\dad0587a-6c19-46c2-a3e3-7d6f29f8e3e4.jpg" /></p><p>and from (23)<img src="1-7400723\881645e8-4bce-4b31-ba8e-747a006a8350.jpg" />. But y must be extremal when <img src="1-7400723\9a18c5a6-2078-48d2-a769-576bd2fcaac1.jpg" /> , therefore:</p><p><img src="1-7400723\fe82d573-b49e-42be-9223-646ad99a89fb.jpg" /></p><p>As it is observed that this solution is equal to exact solution (see <xref ref-type="fig" rid="fig2">Figure 2</xref>).</p><p>Example 3.4. The objective is to find an extremum of the functional</p><disp-formula id="scirp.19048-formula9994"><label>(26)</label><graphic position="anchor" xlink:href="1-7400723\c9881c6f-c14e-4bc5-b5ce-e7d4be95b281.jpg"  xlink:type="simple"/></disp-formula><p>Such that</p><disp-formula id="scirp.19048-formula9995"><label>(27)</label><graphic position="anchor" xlink:href="1-7400723\88b9552f-2ab3-49b3-9e98-7a96bed26d47.jpg"  xlink:type="simple"/></disp-formula><p>and</p><disp-formula id="scirp.19048-formula9996"><label>(28)</label><graphic position="anchor" xlink:href="1-7400723\3b079c57-c590-48d0-8bb0-2662503eb813.jpg"  xlink:type="simple"/></disp-formula><p>With exact solution<img src="1-7400723\c8ecf6fe-f343-4ff3-9e7c-523e8d3c9fc9.jpg" />, <img src="1-7400723\3c2a6b49-a4b7-420e-bdd5-5ab017768a13.jpg" />, see [<xref ref-type="bibr" rid="scirp.19048-ref33">33</xref>]. By having the following auxiliary functional:</p><p><img src="1-7400723\53caab94-4cf5-47d2-8c50-73d7a710da03.jpg" /></p><p>The system of Euler-Lagrange equations is in the form:</p><p><img src="1-7400723\20fbab9b-bc8b-4513-9638-1d774e901d2e.jpg" />.</p><p>So</p><p><img src="1-7400723\22a47644-dba9-4286-8abb-e62751684abd.jpg" /></p><p>By using Homotopy variational iterative method gives:</p><p><img src="1-7400723\49f28346-ff71-46a8-ace2-fe2c3c992614.jpg" /></p><p>Now</p><p><img src="1-7400723\7c3fbdc4-012a-4a26-aa30-4a888101fdc3.jpg" /></p><p>therefore</p><p><img src="1-7400723\864de2bf-84ab-497f-9f82-21e45acba353.jpg" /></p><p>Hence</p><p><img src="1-7400723\1a25589f-31f1-42da-8193-9fa1854d13c9.jpg" /></p><p>and</p><p><img src="1-7400723\a16a9eec-de99-4e62-8ab5-c3cda0ccb2e5.jpg" /></p><p>so</p><p><img src="1-7400723\a9520d52-3ea1-4d7a-9d4e-181538b4b3aa.jpg" /></p><p>So <img src="1-7400723\5c4f6388-dba2-42ed-91d4-037f2e152748.jpg" /> is obtained as:</p><p><img src="1-7400723\a31f3f08-385b-4a64-ad44-b3091976f60c.jpg" /></p><p>and the following iterative equations are obtained:</p><p><img src="1-7400723\0e4a9d50-d6ea-47a0-b411-810a643fd25d.jpg" />,</p><p><img src="1-7400723\a9b306b4-458e-47be-ac88-4b8528430470.jpg" /></p><p>By choosing<img src="1-7400723\dd018307-1240-4b13-9750-51f31c27b450.jpg" />:</p><p><img src="1-7400723\a2e8a85b-c435-4ff0-80cb-425b4b62d2c8.jpg" /></p><p>And by imposing (28) on this functions:</p><p><img src="1-7400723\d77ca39b-6001-4cbe-8b76-296480c4aa54.jpg" /></p><p>from (27):</p><p><img src="1-7400723\4048a6a3-6e85-4aa1-96c8-39f1f0f6a1e0.jpg" /></p><p>And consequently:</p><p><img src="1-7400723\0f0207e1-f3fb-4ecc-a9ee-035cd7f6d9eb.jpg" />.</p><p>which is the exact solution (see <xref ref-type="fig" rid="fig3">Figure 3</xref>).</p></sec></sec><sec id="s4"><title>4. Conclusion</title><p>The He’s variational iterative method is an efficient method for solving various kinds of problems. In this paper variational iterative method is employed for finding the minimum of a functional with moving boundaries and isoperimetric problems. Using He’s variational iterative method the solution of the problem is provided in a closed form. Since this method does not need to the discretize of the variables, there is no computational round off error. Moreover, only a few numbers of iterations are needed to obtain a satisfactory result.</p></sec><sec id="s5"><title>REFERENCES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.19048-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">P. G. Engstrom and U. Brechtken-Manderschied, “Introduction to the Calculus of Variations,” Chapman and Hall/CRC, New York, 1991.</mixed-citation></ref><ref id="scirp.19048-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">A. Saadatmandi and M. Dehghan, “He’s Variational Iteration Method for Solving a Partial Differential Equation Arising in Modeling of Water Waves,” Zeitschrift für Naturforschung, Vol. 64, 2009, pp. 783-787.</mixed-citation></ref><ref id="scirp.19048-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">J. H. He, “Variational Iteration Method a Kind of NonLinear Analytical Technique: Some Examples,” International Journal of Nonlinear Mechanics, Vol. 34, No. 4, 1999, pp. 699-708. doi:10.1016/S0020-7462(98)00048-1</mixed-citation></ref><ref id="scirp.19048-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">M. A. Abdou and A. A. Soliman, “Varitional Iteration Method for Solving Burgers’ and Coupled Burgers’ Equations,” Journal of Computational and Applied Mathmatics, Vol. 181, 2005, pp. 45-251.</mixed-citation></ref><ref id="scirp.19048-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">J. Biazar and H. Ghazvini, “He’s Variational Iteration Method for Solving Hyperbolic Differential Equations,” International Journal of Nonlinear Sciences and Numerical Simulation, Vol. 8, 2007, pp. 311-314.</mixed-citation></ref><ref id="scirp.19048-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">M. Dehghan and M. Tatari, “The Use of He’s Variational Iteration Method for Solving a Fokker Planck Equation,” Physica Scripta, Vol. 74, No. 3, 2006, pp. 310-316. 
doi:10.1088/0031-8949/74/3/003</mixed-citation></ref><ref id="scirp.19048-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">M. Dehghan and A. Saadatmandi, “Variational Iteration Method for Solving the Wave Equation Subject to an Integral Conservation Condition,” Chaos, Solitons and Fractals, Vol. 41, No. 3, 2009, pp. 448-1453. 
doi:10.1016/j.chaos.2008.06.009</mixed-citation></ref><ref id="scirp.19048-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">M. Dehghan and F. Shakeri, “Approximate Solution of a Differential Equation Arising in Astrophysics Using the Variational Iteration Method,” New Astronomy, Vol. 13, No. 1, 2008, pp. 53-59. doi:10.1016/j.newast.2007.06.012</mixed-citation></ref><ref id="scirp.19048-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">M. Dehghan and F. Shakeri, “Application of He’s Variational Iteration Method for Solving the Cauchy Reaction-Diffusion Problem,” Journal of Computational and Applied Mathematics, Vol. 214, No. 2, 2008, pp. 435-446. 
doi:10.1016/j.cam.2007.03.006</mixed-citation></ref><ref id="scirp.19048-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">M. Dehghan and R. Salehi, “The Use of Variational Iteration Method and Adomian Decomposition Method to Solve the Eikonal Equation and Its Application in the Reconstruction Problem,” Communications in Numerical Methods in Engineering, in press.</mixed-citation></ref><ref id="scirp.19048-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">M. Dehghan and M. Tatari, “The Use of Adomian Decomposition Method for Solving Problems in Calculus of Variations,” Mathematical Problems in Engineering, Vol. 2006, 2006, pp. 1-12. doi:10.1155/MPE/2006/65379</mixed-citation></ref><ref id="scirp.19048-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">M. Dehghan and M. Tatari, “Identifying an Unknown Function in a Parabolic Equation with over Specified Data via He’s Variational Iteration Method,” Chaos, Solitons and Fractals, Vol. 36, 2008, pp. 57-166.</mixed-citation></ref><ref id="scirp.19048-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">J. H. He, “Variational Iteration Method for Autonomous Ordinary Differential Systems,” Applied Mathematics and Computation, Vol. 114, No. 2-3, 2000, pp. 115-123. 
doi:10.1016/S0096-3003(99)00104-6</mixed-citation></ref><ref id="scirp.19048-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">J. H. He and X. H. Wu, “Variational Iteration Method: New Development and Applications,” Computers and Mathematics with Applications, Vol. 54, No. 7-8, 2007, pp. 881-894. doi:10.1016/j.camwa.2006.12.083</mixed-citation></ref><ref id="scirp.19048-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">J. H. He, “Variational Iteration Method: Some Recent Results and New Interpretations,” Journal of Computational and Applied Mathematics, Vol. 207, No. 1, 2007, pp. 3-17. doi:10.1016/j.cam.2006.07.009</mixed-citation></ref><ref id="scirp.19048-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">J. H. He, Variational Iteration Method for Delay Differential Equations,” Communications in Nonlinear Science and Numerical Simulation, Vol. 2, No. 4, 1997, pp. 235236. doi:10.1016/S1007-5704(97)90008-3</mixed-citation></ref><ref id="scirp.19048-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">J. H. He, “Approximate Solution of Nonlinear Differential Equations with Convolution Product Nonlinearities,” Computer Methods in Applied Mechanics and Engineering, Vol. 167, No. 1-2, 1998, pp. 69-73. 
doi:10.1016/S0045-7825(98)00109-1</mixed-citation></ref><ref id="scirp.19048-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">J. H. He, “Approximate Analytical Solution for Seepage Flow with Fractional Derivatives in Porous Media,” Computer Methods in Applied Mechanics and Engineering, Vol. 167, No. 1-2, 1998, pp. 57-68. 
doi:10.1016/S0045-7825(98)00108-X</mixed-citation></ref><ref id="scirp.19048-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">M. Inc, “Numerical Simulation of KdV and mKdV Equations with Initial Conditions by the variational Iteration Method,” Chaos, Solitons and Fractals, Vol. 34, No. 4, 2007, pp. 1075-1081. doi:10.1016/j.chaos.2006.04.069</mixed-citation></ref><ref id="scirp.19048-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">S. Momani and S. Abuasad, “Application of He’s Variational Iteration Method to Helmholtz Equation,” Chaos, Solitons and Fractals, Vol. 27, No. 5, 2006, pp. 11191123. doi:10.1016/j.chaos.2005.04.113</mixed-citation></ref><ref id="scirp.19048-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">S. Momani and Z. Odibat, “Analytical Approach to Linear Fractional Partial Differential Equations Arising in Fluid Mechanics. Physics Letters A, Vol. 355, No. 4-5, 2006, pp. 271-279. 
doi:10.1016/j.physleta.2006.02.048</mixed-citation></ref><ref id="scirp.19048-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">H. Ozer, “Application of the Variational Iteration Method to the Boundary Value Problems with Jump Discontinuities Arising in Solid Mechanics,” International Journal of Nonlinear Sciences and Numerical Simulation, Vol. 8, 2007, pp. 513-518.</mixed-citation></ref><ref id="scirp.19048-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Z. M. Odibat and S. Momani, “Application of Variational Iteration Method to Nonlinear Differential Equations of Fractional Order,” International Journal of Nonlinear Sciences and Numerical Simulation, Vol. 7, 2007, pp. 2734. doi:10.1515/IJNSNS.2006.7.1.27</mixed-citation></ref><ref id="scirp.19048-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">H. Sagan, “Introduction to the Calculus of Variations,” Courier Dover Publications, 1992.</mixed-citation></ref><ref id="scirp.19048-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">F. Shakeri and M. Dehghan, “Numerical Solution of the Klein-Gordon Equation via He’s Variational Iteration Method,” Nonlinear Dynamics, Vol. 51, No. 1-2, 2008, pp. 89-97. doi:10.1007/s11071-006-9194-x</mixed-citation></ref><ref id="scirp.19048-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">F. Shakeri and M. Dehghan, “Solution of a Model Describing Biological Species Living Together Using the Variational Iteration Method,” Mathematical and Computer Modeling, Vol. 48, No. 5-6, 2008, pp. 685-699. 
doi:10.1016/j.mcm.2007.11.012</mixed-citation></ref><ref id="scirp.19048-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">M. Tatari and M. Dehghan, “Solution of Problems in Calculus of Variations via He’s Variational Iteration Method,” Physics Letters A, Vol. 362, No. 5-6, 2007, pp. 401-406. doi:10.1016/j.physleta.2006.09.101</mixed-citation></ref><ref id="scirp.19048-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">M. Tatari and M. Dehghan, “Improvement of He’s Variational Iteration Method for Solving Systems of Differential Equations,” Computers and Mathematics with Applications, Vol. 58, No. 11-12, 2009, pp. 2160-2166. 
doi:10.1016/j.camwa.2009.03.081</mixed-citation></ref><ref id="scirp.19048-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">M. Tatari and M. Dehghan, “On the Convergence of He's Variational Iteration Method,” Journal of Computational and Applied Mathematics, Vol. 207, No. 1, 2007, pp. 121128. doi:10.1016/j.cam.2006.07.017</mixed-citation></ref><ref id="scirp.19048-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">S. A. Youse and M. Dehghan, “The Use of He’s Variational Iteration Method for Solving Variational Problems,” International Journal of Computer Mathematics, Vol. 87, No. 6, 2010. doi:10.1080/00207160802283047</mixed-citation></ref><ref id="scirp.19048-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">A. M. Wazwaz, 2007. “A Comparison between the Variational Iteration Method and Adomian Decomposition Method,” Journal of Computational and Applied Mathematics, Vol. 207, No. 1, pp. 129-136. 
doi:10.1016/j.cam.2006.07.018c</mixed-citation></ref><ref id="scirp.19048-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">R. Memarbashi, “Variational Problems with Moving Boundaries Using Decomposition Method,” Mathematical Problems in Engineering, Vol. 2007, 2007, Article ID 10120.</mixed-citation></ref><ref id="scirp.19048-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">M. L. Krasnov, G. I. Makarenko and A. I. Kiselev, “Problems and Exercises in the Calculus of Variations,” Mir Publisher, Moscow, 1964. </mixed-citation></ref></ref-list></back></article>