DOCSLIB.ORG

Green's Function Solution of Non-Homogenous Regular Sturm

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Green's Function Solution of Non-Homogenous Regular Sturm Computa & tio d n e a li Hassana, J Appl Computat Math 2017, 6:3 l p M p a Journal of A t h DOI: 10.4172/2168-9679.1000362 f e o m l a a n t r ISSN: 2168-9679i c u s o J Applied & Computational Mathematics Research Article Open Access Green’s Function Solution of Non-Homogenous Regular Sturm-Liouville Problem Abdelgabar Adam Hassana* Department of Mathematics, College of Science and Arts at Tabrjal, AlJouf University, Kingdom of Saudi Arabia Abstract In this paper, we propose a new method called exp(−ϕ(ξ)) fractional expansion method to seek traveling wave solutions of the nonlinear fractional Sharma-Tasso-Olver equation. The result reveals that the method together with the new fractional ordinary differential equation is a very influential and effective tool for solving nonlinear fractional partial differential equations in mathematical physics and engineering. The obtained solutions have been articulated by the hyperbolic functions, trigonometric functions and rational functions with arbitrary constants. Introduction d dy d22 ydp() x dy d y =+=+ So px() px() 22p()x p(x) bxz() The series solution of differential equation yields an infinite series dx dx dx dx dx dx which often converges slowly. Thus it is difficult to obtain an insight Thus equation (3) can be written as into over-all behavior of the solution [1,2]. The Green’s function d dy approach would allow us to have an integral representation of the px() ++ qxy()()λ rxy =0 (4) solution instead of an infinite series. dx dx To obtain the filed y, caused by distributed source we calculate the Where q(x)=p(x)c(x) and r(x)=p(x)d(x) effect of each elementary portion of source and add (integral) them all. Equation in form (4) is known as Sturm-Liouville equation. Satisfy If G(r,r0) is the field at the observers point r caused by a unit source the boundary conditions at the source point r0, then the field at r caused by distribution f(r0) is Regular Sturm-Liouville Problem the integral of f(r0) G(r,r0) over the whole range of r0 occupied by the course. The function G is called Green’s function [3-5]. In case p(a)≠0 and p(b)≠0, p(x), q(x), r(x) are continuous, the The Green’s function is powerful tool of mathematical method Sturm-Liouville equation (4) can be expressed as which used I solving linear non-homogenous differential equation L[y]=λr(x)y (5) (ordinary and partial) [6-9]. d dy Where L= px() + qx() (6) Preliminaries dx dx If the above equation is associated with the following boundary Sturm-Liouville problem condition Consider a linear second order differential equation α1y(a)+α2y (a)=0 d2 y dy + ++λ = (1) Ax() 2 Bx() C(x) y D(x) y 0 β1y(b)+β2y′(b)=0 (7) dx∂ x Where α +α ≠0 and β +β ≠0 Where λ is a parameter to be determined by the boundary 1 ′ 2 1 2 conditions? A(x) is positive continuous function, then by dividing The equation (4) and the boundary condition (7) are called regular every term by A(x), equation (1) can be written as Sturm -Liouville problem (RSLP). d2 y dy (a) Singular Sturm-Liouville problem +Bx() ++ C(x) yλ D(x) y = 0 (2) dx2 ∂ x Consider the equation Bx() Cx() Dx() Where b(x)= , c(x)= and d(x)= L[y]+λr(x)y=0 a<x<b (8) Ax() Ax() Ax() Let us define integrating factorp (x) by x *Corresponding author: Abdelgabar Adam Hassana, Department of Mathematics, Px()()= exp∫ bξξ d College of Science and Arts at Tabrjal, AlJouf University, Kingdom of Saudi Arabia, a Tel: +966 14 624 7493; E-mail: [email protected] Multiplying equation (2) by p(x), we have Received April 05, 2017; Accepted April 25, 2017; Published May 11, 2017 d2 y dy Citation: Hassana AA (2017) Green’s Function Solution of Non-Homogenous + ++λ = Regular Sturm-Liouville Problem. J Appl Computat Math 6: 362. doi: 10.4172/2168- px() 2 pxbx()() p(x)c()() x y pxd (x) y 0 (3) dx ∂x 9679.1000362 xx bd()ζζ bd() ζζ x Copyright: © 2017 Hassana AA. This is an open-access article distributed under dp() x dd∫∫ Since = eaa= e b(ζζ )d = pxbx()() the terms of the Creative Commons Attribution License, which permits unrestricted ∫ use, distribution, and reproduction in any medium, provided the original author and dx dx dx a source are credited. J Appl Computat Math, an open access journal Volume 6 • Issue 3 • 1000362 ISSN: 2168-9679 Citation: Hassana AA (2017) Green’s Function Solution of Non-Homogenous Regular Sturm-Liouville Problem. J Appl Computat Math 6: 362. doi: 10.4172/2168-9679.1000362 Page 2 of 6 Where L is defined by (6), p(x) is smooth and r(x) is positive, then n δ = the Sturm- Liouville problem is called singular if one of the following n ( x) 22 π 1+ nx situations is occurred. ( ) n −nx22 (i) If p(a)=0 or p(b)= 0 or both δn ( xe) = (17) π (ii)The interval (a,b) is infinite. sin2 (nx) δ = n ( x) 2 etc (b)Eigenvalue and eigenfunction nxπ The Eigenvalue from equation (4) defining by a Sturm-Liouville (c) Some important properties of Dirac delta function operator can be expressed as Property (1): Symmetry −1 d dy λ = px( ) + qx( ) (9) δ(−x)=δ(x) (18) r( x) dx dx The non-trival solutions that satisfy the equation and boundary Proof conditions are called eigenfunctions. Therefore the eigenfunction of Let ζ=−x, then dx=−dζ the Sturm-Liouville problem from complete sets of orthogonal bases for the function space is which the weight function is r(x). We can write: +∞ −∞ +∞ The Dirac delta function ∫∫∫fxxdxfdfdf( )δ(−=) (ζδζ) ( ) ζ ≈( ζδζ) ( ) ζ =(0) (19) −∞ +∞ −∞ The delta function is defined as +∞ 0 x ≠ξ But, f( x)δ ( x) dx= f (0) (20) (10) ∫ δξ( x −=) −∞ ∞=x ξ Therefore, from (19) and (20), we conclude that δ(−x)=δ(x) But such that the integral of δ (x−ζ) is normalized to unit +∞ Property (2): Scaling ∫ δξ( x−=) dx 1 (11) 1 −∞ δδ(ax)= ( x) (21) a In fact the first operator where Dirac used the delta function is the Proof integration 1 Let ζ=ax, then dx= dζ +∞ a ∫ f( x)δξ( x− ) dx (12) If a>0, then −∞ +∞ −∞ ζζ1 1+∞ 10 Where f(x) is a continuous function, we have to find the value of f( x)δ( ax) dx= f δζ( ) d ζ= f δζ( ) d ζ= f f (0) ∫∫a a a ∫ a aa the integration (12). Since δ (xζ) is zero for x≠ζ, the limit of integration −∞ +∞ −∞ may be change to ζ−ε and ζ+ε, where ε is a small positive number, f(x) +∞ +∞ 11 1 is continuous at x−ζ, it’s values within the interval (ζ−ε,ζ−ε) will not Since ∫∫f(ζδ) ( x) dx= f( x)δ( x) dx= f (0) different much from f(ζ), approximately that: −∞ aa−∞ a +∞ ξξ+∞ +∞ Therefore: f xδξ x−= dx f x δξ x −≈ dx f ξδξ x − dx (13) ∫∫∫( ) ( ) ( ) ( ) ( ) ( ) 1 −∞ ξξ−∞ −∞ δδ(ax)= ( x) a With the approximation improving as ε approaches zero. Similarly for a<0, 1 From (11), we have δδ(ax)= ( x) −a +∞ ξ +∞ then I=∫∫δξ( x −=) dx δξ( x −=) dx 1 (14) −∞ ξ −∞ We can write: From all values of ε, then by letting ε→0, we can exactly have 1 δδ(ax)= ( x) +∞ a ∫ f( x)δξ( x−=) dx f ( ζ) (15) −∞ Property (3) 22 1 Despite the delta function considered as fundamental role in δ(x− a) = δδ( xa +) +( xa −) (22) electrical engineering and quantum mechanics, but no conventional 2a could be found that satisfies (10) and (11), then the delta function Proof sought to be view as the limit of the sequence of strongly peaked The argument of this function goes to zero when x = a and x=−a, function δn(x)such that wherefore +∞ +∞ 22 δδ( xx)= lim n ( ) (16) fζδ x− a dx = f x δ x + a + δ x − a dx n→∞ ∫∫( ) ( ) ( )( ) ( ) −∞ −∞ As Only at the zero of the argument of the delta function that is: J Appl Computat Math, an open access journal Volume 6 • Issue 3 • 1000362 ISSN: 2168-9679 Citation: Hassana AA (2017) Green’s Function Solution of Non-Homogenous Regular Sturm-Liouville Problem. J Appl Computat Math 6: 362. doi: 10.4172/2168-9679.1000362 Page 3 of 6 +∞ −ae + ae+ 22 22 22 αGa( ) += αζ G′( a,0) ∫∫f(ζδ) ( x−= a) dx f( x)δ( x −+ a) dx ∫ f( x)δ( x − a) dx (23) 12 (31) −∞ −ae − ae− β12Gb( )+= βζ G′( b,0) Near the two zeros x2−a2 can be approximated as: Now consider the region a≤x<ζ. 22 (−2,axa)( +) x →− a ( x− a) =−( xaxa)( +=) Let y1(x) be a nontrivial solution at x=a, i.e (+2,axa)( −) x →+ a α y (a)+α y (a)=0 (32) In the limit as ε→0 the integral (23) becomes: 1 1 2 +∞ −ae + ae+ Then α y (a,ζ)+α G (a,ζ)=0 (33) 22 1 1 ′ 2 ∫∫f(ζδ) ( x− a) dx = f(x)δ −( 2 a)( x + a) dx + ∫ f(x)δ ( 2 a)( x − a) dx −∞ −ae − ae− The wronskian of y1′ and G must vanish at x=a or 11ae++ae = f(x)δδ( x) ( x + a) + f(x)δ( x− a) dx ∫∫ y1(a)G (a,ζ)+y 1(a)G(a,ζ)=0 22aaae−−ae +∞ 1 So G(x,′ζ)= u1y1′(x) for a≤x<ζ (34) =∫ fx( ) δδ[ xa−+] ( xadx −) −∞ 2a Where u1 is an arbitrary constant.
Load more

Recommended publications
  • Ordinary Differential Equations
    Ordinary Differential Equations for Engineers and Scientists Gregg Waterman Oregon Institute of Technology c 2017 Gregg Waterman This work is licensed under the Creative Commons Attribution 4.0 International license. The essence of the license is that You are free to: Share copy and redistribute the material in any medium or format • Adapt remix, transform, and build upon the material for any purpose, even commercially. • The licensor cannot revoke these freedoms as long as you follow the license terms. Under the following terms: Attribution You must give appropriate credit, provide a link to the license, and indicate if changes • were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use. No additional restrictions You may not apply legal terms or technological measures that legally restrict others from doing anything the license permits. Notices: You do not have to comply with the license for elements of the material in the public domain or where your use is permitted by an applicable exception or limitation. No warranties are given. The license may not give you all of the permissions necessary for your intended use. For example, other rights such as publicity, privacy, or moral rights may limit how you use the material. For any reuse or distribution, you must make clear to others the license terms of this work. The best way to do this is with a link to the web page below. To view a full copy of this license, visit https://creativecommons.org/licenses/by/4.0/legalcode.
    [Show full text]
  • New Dirac Delta Function Based Methods with Applications To
    New Dirac Delta function based methods with applications to perturbative expansions in quantum field theory Achim Kempf1, David M. Jackson2, Alejandro H. Morales3 1Departments of Applied Mathematics and Physics 2Department of Combinatorics and Optimization University of Waterloo, Ontario N2L 3G1, Canada, 3Laboratoire de Combinatoire et d’Informatique Math´ematique (LaCIM) Universit´edu Qu´ebec `aMontr´eal, Canada Abstract. We derive new all-purpose methods that involve the Dirac Delta distribution. Some of the new methods use derivatives in the argument of the Dirac Delta. We highlight potential avenues for applications to quantum field theory and we also exhibit a connection to the problem of blurring/deblurring in signal processing. We find that blurring, which can be thought of as a result of multi-path evolution, is, in Euclidean quantum field theory without spontaneous symmetry breaking, the strong coupling dual of the usual small coupling expansion in terms of the sum over Feynman graphs. arXiv:1404.0747v3 [math-ph] 23 Sep 2014 2 1. A method for generating new representations of the Dirac Delta The Dirac Delta distribution, see e.g., [1, 2, 3], serves as a useful tool from physics to engineering. Our aim here is to develop new all-purpose methods involving the Dirac Delta distribution and to show possible avenues for applications, in particular, to quantum field theory. We begin by fixing the conventions for the Fourier transform: 1 1 g(y) := g(x) eixy dx, g(x)= g(y) e−ixy dy (1) √2π √2π Z Z To simplify the notation we denote integration over the real line by the absence of e e integration delimiters.
    [Show full text]
  • Second Order Linear Differential Equations Y
    Second Order Linear Differential Equations Second order linear equations with constant coefficients; Fundamental solutions; Wronskian; Existence and Uniqueness of solutions; the characteristic equation; solutions of homogeneous linear equations; reduction of order; Euler equations In this chapter we will study ordinary differential equations of the standard form below, known as the second order linear equations : y″ + p(t) y′ + q(t) y = g(t). Homogeneous Equations : If g(t) = 0, then the equation above becomes y″ + p(t) y′ + q(t) y = 0. It is called a homogeneous equation. Otherwise, the equation is nonhomogeneous (or inhomogeneous ). Trivial Solution : For the homogeneous equation above, note that the function y(t) = 0 always satisfies the given equation, regardless what p(t) and q(t) are. This constant zero solution is called the trivial solution of such an equation. © 2008, 2016 Zachary S Tseng B-1 - 1 Second Order Linear Homogeneous Differential Equations with Constant Coefficients For the most part, we will only learn how to solve second order linear equation with constant coefficients (that is, when p(t) and q(t) are constants). Since a homogeneous equation is easier to solve compares to its nonhomogeneous counterpart, we start with second order linear homogeneous equations that contain constant coefficients only: a y″ + b y′ + c y = 0. Where a, b, and c are constants, a ≠ 0. A very simple instance of such type of equations is y″ − y = 0 . The equation’s solution is any function satisfying the equality t y″ = y. Obviously y1 = e is a solution, and so is any constant multiple t −t of it, C1 e .
    [Show full text]
  • 5 Mar 2009 a Survey on the Inverse Integrating Factor
    A survey on the inverse integrating factor.∗ Isaac A. Garc´ıa (1) & Maite Grau (1) Abstract The relation between limit cycles of planar differential systems and the inverse integrating factor was first shown in an article of Giacomini, Llibre and Viano appeared in 1996. From that moment on, many research articles are devoted to the study of the properties of the inverse integrating factor and its relation with limit cycles and their bifurcations. This paper is a summary of all the results about this topic. We include a list of references together with the corresponding related results aiming at being as much exhaustive as possible. The paper is, nonetheless, self-contained in such a way that all the main results on the inverse integrating factor are stated and a complete overview of the subject is given. Each section contains a different issue to which the inverse integrating factor plays a role: the integrability problem, relation with Lie symmetries, the center problem, vanishing set of an inverse integrating factor, bifurcation of limit cycles from either a period annulus or from a monodromic ω-limit set and some generalizations. 2000 AMS Subject Classification: 34C07, 37G15, 34-02. Key words and phrases: inverse integrating factor, bifurcation, Poincar´emap, limit cycle, Lie symmetry, integrability, monodromic graphic. arXiv:0903.0941v1 [math.DS] 5 Mar 2009 1 The Euler integrating factor The method of integrating factors is, in principle, a means for solving ordinary differential equations of first order and it is theoretically important. The use of integrating factors goes back to Leonhard Euler. Let us consider a first order differential equation and write the equation in the Pfaffian form ω = P (x, y) dy Q(x, y) dx =0 .
    [Show full text]
  • Handbook of Mathematics, Physics and Astronomy Data
    Handbook of Mathematics, Physics and Astronomy Data School of Chemical and Physical Sciences c 2017 Contents 1 Reference Data 1 1.1 PhysicalConstants ............................... .... 2 1.2 AstrophysicalQuantities. ....... 3 1.3 PeriodicTable ................................... 4 1.4 ElectronConfigurationsoftheElements . ......... 5 1.5 GreekAlphabetandSIPrefixes. ..... 6 2 Mathematics 7 2.1 MathematicalConstantsandNotation . ........ 8 2.2 Algebra ......................................... 9 2.3 TrigonometricalIdentities . ........ 10 2.4 HyperbolicFunctions. ..... 12 2.5 Differentiation .................................. 13 2.6 StandardDerivatives. ..... 14 2.7 Integration ..................................... 15 2.8 StandardIndefiniteIntegrals . ....... 16 2.9 DefiniteIntegrals ................................ 18 2.10 CurvilinearCoordinateSystems. ......... 19 2.11 VectorsandVectorAlgebra . ...... 22 2.12ComplexNumbers ................................. 25 2.13Series ......................................... 27 2.14 OrdinaryDifferentialEquations . ......... 30 2.15 PartialDifferentiation . ....... 33 2.16 PartialDifferentialEquations . ......... 35 2.17 DeterminantsandMatrices . ...... 36 2.18VectorCalculus................................. 39 2.19FourierSeries .................................. 42 2.20Statistics ..................................... 45 3 Selected Physics Formulae 47 3.1 EquationsofElectromagnetism . ....... 48 3.2 Equations of Relativistic Kinematics and Mechanics . ............. 49 3.3 Thermodynamics and Statistical Physics
    [Show full text]
  • FIRST-ORDER ORDINARY DIFFERENTIAL EQUATIONS III: Numerical and More Analytic Methods
    FIRST-ORDER ORDINARY DIFFERENTIAL EQUATIONS III: Numerical and More Analytic Methods David Levermore Department of Mathematics University of Maryland 30 September 2012 Because the presentation of this material in lecture will differ from that in the book, I felt that notes that closely follow the lecture presentation might be appreciated. Contents 8. First-Order Equations: Numerical Methods 8.1. Numerical Approximations 2 8.2. Explicit and Implicit Euler Methods 3 8.3. Explicit One-Step Methods Based on Taylor Approximation 4 8.3.1. Explicit Euler Method Revisited 4 8.3.2. Local and Global Errors 4 8.3.3. Higher-Order Taylor-Based Methods (not covered) 5 8.4. Explicit One-Step Methods Based on Quadrature 6 8.4.1. Explicit Euler Method Revisited Again 6 8.4.2. Runge-Trapezoidal Method 7 8.4.3. Runge-Midpoint Method 9 8.4.4. Runge-Kutta Method 10 8.4.5. General Runge-Kutta Methods (not covered) 12 9. Exact Differential Forms and Integrating Factors 9.1. Implicit General Solutions 15 9.2. Exact Differential Forms 16 9.3. Integrating Factors 20 10. Special First-Order Equations and Substitution 10.1. Linear Argument Equations (not covered) 25 10.2. Dilation Invariant Equations (not covered) 26 10.3. Bernoulli Equations (not covered) 27 10.4. Substitution (not covered) 29 1 2 8. First-Order Equations: Numerical Methods 8.1. Numerical Approximations. Analytic methods are either difficult or impossible to apply to many first-order differential equations. In such cases direction fields might be the only graphical method that we have covered that can be applied.
    [Show full text]
  • Phase-Space Matrix Representation of Differential Equations for Obtaining the Energy Spectrum of Model Quantum Systems
    Phase-space matrix representation of differential equations for obtaining the energy spectrum of model quantum systems Juan C. Morales, Carlos A. Arango1, a) Department of Chemical Sciences, Universidad Icesi, Cali, Colombia (Dated: 27 August 2021) Employing the phase-space representation of second order ordinary differential equa- tions we developed a method to find the eigenvalues and eigenfunctions of the 1- dimensional time independent Schr¨odinger equation for quantum model systems. The method presented simplifies some approaches shown in textbooks, based on asymptotic analyses of the time-independent Schr¨odinger equation, and power series methods with recurrence relations. In addition, the method presented here facilitates the understanding of the relationship between the ordinary differential equations of the mathematical physics and the time independent Schr¨odinger equation of physical models as the harmonic oscillator, the rigid rotor, the Hydrogen atom, and the Morse oscillator. Keywords: phase-space, model quantum systems, energy spectrum arXiv:2108.11487v1 [quant-ph] 25 Aug 2021 a)Electronic mail: [email protected] 1 I. INTRODUCTION The 1-dimensional time independent Schr¨odinger equation (TISE) can be solved analyt- ically for few physical models. The harmonic oscillator, the rigid rotor, the Hydrogen atom, and the Morse oscillator are examples of physical models with known analytical solution of the TISE (1). The analytical solution of the TISE for a physical model is usually obtained by using the ansatz of a wavefunction as a product of two functions, one of these functions acts as an integrating factor (2), the other function produces a differential equation solv- able either by Frobenius series method or by directly comparing with a template ordinary differential equation (ODE) with known solution (3).
    [Show full text]
  • Mathematics for Chemistry Prof. Madhav Ranganathan Department of Chemistry Indian Institute of Technology, Kanpur
    Mathematics for Chemistry Prof. Madhav Ranganathan Department of Chemistry Indian Institute of Technology, Kanpur Module - 05 Lecture - 23 Integrating Factors Today I will be talking about the Integrating Factor method, which is a method for solving differential equations; first order differential equations which cannot be expressed as exact differentials. (Refer Slide Time: 00:22) (Refer Slide Time: 00:27) (Refer Slide Time: 02:59) So, just to remind ourselves about exact and in exact differentials. We said that, we can write our first order differential equation as M times dx plus N times dy equal to 0; where M and N are functions of x and y. And what we said is at an exact differential implies. So, exact implies dou N by dou y equal to dou N by dou x. So, then you say it like that. If it is not exact; so if not exact then can we find some function alpha of x y such that alpha times M dx plus alpha times N dy equal to 0 is an exact differential. So, what I mean is that the left hand side is an exact differential. So, what I did is I just took this equation multiplied it by alpha. So, the right hand side will be 0 into alpha which will still be 0. Now, what I have is alpha M and alpha N instead of M and N. Now if this is an exact differential. So, can we find such an alpha? So, this is the main goal of this idea of integrating factor.
    [Show full text]
  • Solving 1St Order Odes
    Solving First Order ODEs Table of contents Solving First Order ODEs ............................................... 1 1. Introduction ..................................... ................. 1 Side: Two ways of writing a semilinear 1st order ODE . .............. 2 2. Integrating factors and separation of variables . ......................... 2 Side: One must be careful when multiplying an ODE by a function .............. 3 Side:Theroadahead .................................. ......... 3 3. Linear1storderequations .......................... ................... 4 Summary: Solving 1st order linear equations . ............... 6 4. Separableequations ............................... .................. 7 Summary: Solvingseparable equations . .............. 9 5. Exactequations ................................... ................ 9 Side:Theoperator“d” ................................ ........... 9 6. Non-exactequations ............................... .................. 12 7. Solvable equations in disguise: Transformation method . ......................... 18 7.1. Homogeneousequation ............................ ............... 18 7.2. dy G a x b y ................................................ 19 dx = ( + ) 7.3. Bernoulliequation .............................. ................ 19 7.4. Equations with linear coefficients (NOT REQUIRED FOR 334 FALL2010) ......... 20 8. Direction field, numerics, and well-posedness theory . ......................... 22 8.1. Directionfield .................................. ............... 22 8.2. Numerics ......................................
    [Show full text]
  • Analytic Solutions for Third Order Ordinary Differential Equations
    2019 HAWAII UNIVERSITY INTERNATIONAL CONFERENCES SCIENCE, TECHNOLOGY & ENGINEERING, ARTS, MATHEMATICS & EDUCATION JUNE 5 - 7, 2019 HAWAII PRINCE HOTEL WAIKIKI, HONOLULU, HAWAII ANALYTIC SOLUTIONS FOR THIRD ORDER ORDINARY DIFFERENTIAL EQUATIONS BECCAR-VARELA, MARIA P. ET AL UNIVERSITY OF TEXAS AT EL PASO EL PASO, TEXAS Dr. Maria P. Beccar-Varela Department of Mathematical Sciences University of Texas at El Paso El Paso, Texas Mr. Md Al Masum Bhuiyan Mr. Osei K. Tweneboah Computational Science Program University of Texas at El Paso El Paso, Texas Dr. Maria C. Mariani Department of Mathematical Sciences and Computational Science Program University of Texas at El Paso El Paso, Texas Analytic Solutions for Third Order Ordinary Differential Equations Synopsis: This work studies an analytic approach for solving higher order ordinary differential equations (ODEs). We develop alternate techniques for solving third order ODEs and discuss possible generalizations to higher order ODEs. The techniques are effective for solving complex ODEs and could be used in other application of sciences such as physics, engineering, and applied sciences. Analytic Solutions for Third Order Ordinary Differential Equations Maria P. Beccar-Varela ∗, Md Al Masum Bhuiyany, Maria C. Marianiz and Osei K. Tweneboahx Abstract This paper focuses on an analytic approach for solving higher order ordinary differential equations (ODEs). We develop a self-adjoint for- mulation and integrating-factor techniques to solve third order ODEs. The necessary conditions for ODEs to be self-adjoint are also pro- vided. Under these conditions, we find the analytic solution of the ODEs. The solutions produced in this work are exact unlike numeri- cal solutions which have approximation errors.
    [Show full text]
  • Solving Linear First-Order Differential Equations Leonard Euler's Integrating Factor Method
    Ursinus College Digital Commons @ Ursinus College Transforming Instruction in Undergraduate Differential Equations Mathematics via Primary Historical Sources (TRIUMPHS) Summer 2020 Solving Linear First-Order Differential Equations Leonard Euler's Integrating Factor Method Adam E. Parker Wittenberg University Follow this and additional works at: https://digitalcommons.ursinus.edu/triumphs_differ Part of the Curriculum and Instruction Commons, Educational Methods Commons, Higher Education Commons, and the Science and Mathematics Education Commons Click here to let us know how access to this document benefits ou.y Recommended Citation Parker, Adam E., "Solving Linear First-Order Differential Equations Leonard Euler's Integrating Factor Method" (2020). Differential Equations. 4. https://digitalcommons.ursinus.edu/triumphs_differ/4 This Course Materials is brought to you for free and open access by the Transforming Instruction in Undergraduate Mathematics via Primary Historical Sources (TRIUMPHS) at Digital Commons @ Ursinus College. It has been accepted for inclusion in Differential Equations by an authorized administrator of Digital Commons @ Ursinus College. For more information, please contact [email protected]. Solving Linear First-Order Differential Equations Leonard Euler’s Integrating Factor Method Adam E. Parker∗ May 26, 2021 1 Introduction In 1926, British mathematician E. L. Ince (1891–1941) described the typical evolution of solution techniques from calculus (and differential equations and science in general).1 1111111111111111111111111111111111111111 The early history of the infinitesimal calculus abounds in instances of problems solved through the agency of what were virtually differential equations; it is even true to say that the problem of integration, which may be regarded as the solution of the simplest of all types of differential equations, was a practical problem even in the middle of the sixteenth century.
    [Show full text]
  • Ordinary Differential Equations Math 22B-002, Spring 2017 Final Exam: Solutions 1. [15 Pts.] Solve the Initial Value Problem
    Ordinary Differential Equations Math 22B-002, Spring 2017 Final Exam: Solutions 1. [15 pts.] Solve the initial value problem 6y00 − y0 − y = 0; y(0) = 10; y0(0) = 0: Solution. • The characteristic equation is 6r2 −r −1 = 0 with roots r = 1=2; −1=3, so the general solution is x=2 −x=3 y(x) = c1e + c2e : • The initial conditions are satisfied if 1 1 c + c = 10; c − c = 0 1 2 2 1 3 2 whose solution is c1 = 4, c2 = 6, so y(x) = 4ex=2 + 6e−x=3: 1 2. [20 pts.] (a) Let y0 is an arbitrary constant. Find the solution y(x) of the initial-value problem 0 y − 2xy = x; y(0) = y0: (b) For what initial value y0 does the solution remain bounded as x ! +1? What is the solution in that case? Solution. • The ODE is first-order, linear, and nohomogeneous, so we use the in- tegrating factor method. • The integrating factor with coefficient p(x) = −2x is given by Z µ(x) = exp − 2x dx = e−x2 : • Multiplying the ODE by the integrating factor and writing the left- hand side as an exact derivative, we get that 2 0 2 e−x y = xe−x ; so 2 Z 2 1 2 e−x y(x) = xe−x dx + c = − ex + c; 2 and the general solution is 1 2 y(x) = − + cex : 2 • The initial condition gives y0 = −1=2 + c, so c = y0 + 1=2 and 1 1 2 y(x) = − + y + ex : 2 0 2 • (b) The solution remains bounded as x ! 1 if c = 0 and y0 = −1=2.
    [Show full text]