Ordinary Differential Equations Math 22B-002, Spring 2017 Final Exam: Solutions 1. [15 Pts.] Solve the Initial Value Problem

Total Page:16

File Type:pdf, Size:1020Kb

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. In that case, the solution is y(x) = −1=2. Remark. Alternatively, we could use variation of parameters: the solution of the homogeneous equation is y(x) = cex2 , so write y(x) = u(x)ex2 and solve for u(x); or we could observe that yp(x) = −1=2 is a particular solution, and x2 write y(x) = yp(x) + ce . 2 3. [20 pts.] (a) Find the general solution of the ODE y00 − 8y0 + 25y = 0: (b) Find the general solution of the ODE y00 + 25y = 0: Solution. • (a) The characteristic equation is r2 − 8r + 25 = 0 with roots r = 4 ± 3i It follows that the general solution is 4x 4x y(x) = c1e cos 3x + c2e sin 3x where c1, c2 are arbitrary constants. • (b) The characteristic equation is r2 + 25 = 0 with roots r = ±5i: It follows that the general solution is y(x) = c1 cos 5x + c2 sin 5x where c1, c2 are arbitrary constants. 3 4. [20 pts.] (a) Find a particular solution of the ODE y00 + y = ex + x2: (b) Find a particular solution of the ODE y00 − y = ex + x2: Solution. • (a) We look for a particular solution of the form x 2 yp(x) = Ae + Bx + C: • Then 00 x 2 yp + yp = 2Ae + Bx + 2B + C; so we take A = 1=2, B = 1, and 2B + C = 0. It follows that C = −2 and 1 y (x) = ex + x2 − 2: p 2 (b) The function ex is a solution of the homogeneous equation, so we look for a particular solution of the form x 2 yp(x) = Axe + Bx + C: • Then 00 x 2 yp − yp = 2Ae − Bx + 2B − C; so we take A = 1=2, B = −1, and 2B − C = 0. It follows that C = −2 and 1 y (x) = xex − x2 − 2: p 2 4 5. [25 pts.] Suppose that the 2 × 2 matrix A has the following eigenvalues λk and eigenvectors ~rk: 1 2 λ = −2; ~r = ; λ = 1; ~r = : 1 1 3 2 2 −1 (a) Write down the general solution for (x(t); y(t))T of the system x 0 x = A : y y (b) Sketch the phase plane of this system. Include in your sketch the trajec- tory that passes through the point (0; 2). (c) Sketch the graphs of x(t) and y(t) versus t for the solution of this system that satisfies the initial condition x(0) = 0, y(0) = 2. Solution. • (a) The general solution is x 1 2 = c e−2t + c et y 1 3 2 −1 where c1, c2 are arbitrary constants. • (b) The origin is a saddle point, The stable direction, in which (x; y) ! (0; 0) as t ! +1, is (1; 3); the unstable direction, in which (x; y) ! (0; 0) as t ! −∞, is (2; −1). • A numerical plot of the phase plane (using the matlab script pplane8.m) is shown in Figure 1. • (c) Graphs of x(t) and y(t) versus t are shown in Figure 2. 5 Figure 1: Phase plane for 5(b). Figure 2: Graphs for 5(c). 6 Figure 3: Phase plane for 6(b). 6. [25 pts.] (a) Find the general solution of the 2 × 2 system x 0 1 2 x = : y −5 −1 y Write your answer in complex form. (b) Sketch the phase plane of this system. Solution. • The characteristic polynomial of the matrix is 1 − λ 2 2 = λ + 9; −5 −1 − λ so the eigenvalues satisfy λ2 + 9 = 0, and λ = ±3i. • If λ = 3i, then the eigenvector ~r = (p; q)T satisfies (A − 3iI)~r = 0, or 1 − 3i 2 p 0 = : −5 −1 − 3i q 0 7 One solution is p 2 = : q −1 + 3i Any non-zero (complex) scalar multiple of this vector would do just as well. • The eigenvector for λ = −3i is the conjugate of the eigenvector for λ = 3i, so the general solution is x 2 2 = c e3it + c e−3it y 1 −1 + 3i 2 −1 − 3i where c1, c2 are arbitrary (complex) constants. To get real-valued solutions, we take c1 and c2 to be complex conjugates. • The origin is a center with closed (periodic) orbits. Since x0 = x + 2y is positive when x; y > 0 and negative when x; y < 0, the direction of trajectories is clockwise. A numerical plot of the phase plane is shown in Figure 3. 8 2 7. [25 pts.] (a) Let m; c; k; F0 > 0 be positive constants with c > 4mk. Find the general solution for y(t) of the ODE 00 0 my + cy + ky = F0: (b) How does y(t) behave as t ! +1? (c) Give a physical interpretation of this ODE and the behavior of its solu- tions. Solution. • A particular solution of the ODE is the constant yp = F0=k. • The characteristic equation of the homogeneous ODE is mr2+cr+k = 0 with roots p −c ± c2 − 4mk r = = −α; −β: 2 2 pBoth of these roots are negative real numbers, since c − 4mk > 0 and c2 − 4mk < c. • It follows that the general solution of the ODE (particular solution + general solution of the homogeneous equation) is given by F y(t) = 0 + c e−αt + c e−βt k 1 2 where c1, c2 are arbitrary constants. • (b) We have y(t) ! F0=k as t ! +1. • (c) This ODE describes a mass m on a spring with spring constant k and damping proportional to velocity with damping constant c. The mass is subject to a constant external force F0. • The displacement yp = F0=k corresponds to the equilibrium location of the mass, when the external force F0 is balanced by the restoring force −kyp of the spring. As t ! 1, the mass relaxes to its equilibrium state due to damping. 9 8. [25 pts.] Suppose that a hot metal bar with absolute temperature T (t) at time t cools by emitting thermal radiation. According to the Stefan- Boltzmann law, the bar radiates heat energy at rate proportional to T 4. If the initial temperature of the bar is T0 > 0, then T (t) satisfies the initial value problem dT = −kT 4;T (0) = T ; dt 0 where k > 0 is a constant. (a) Solve for T (t). (b) Let t1=2 be the time it takes for the bar to cool to half its initial temper- ature. Find an expression for t1=2 in terms of k and T0. (c) Does t1=2 increase, decrease, or stay the same as T0 increases? How does this compare with the behavior of t1=2 for an object that cools according to Newton's law of cooling (dT=dt = −kT )? Solution. • (a) The ODE is nonlinear and separable. Separating variables, we get Z dT Z = − k; dt + c T 4 Integration gives 1 − = −kt + c: 3T 3(t) 3 The initial condition gives c = −1=3T0 , and solving for T , we find that T T (t) = 0 : 3 1=3 (1 + 3T0 kt) • (b) The temperature of the bar is half the initial temperature when 3 1=3 1 + 3T0 kt1=2 = 2: Solving for t1=2, we get that 7 t1=2 = 3 : 3T0 k • (c) The time t1=2 decreases as T0 increases (the cooling rate is larger at −kt higher temperatures). For Newton's law of cooling, T (t) = T0e and the decay time t1=2 = ln 2=k is independent of T0. 10 9. [25 pts.] Consider the following ODE for y(x): y0 + sin2(x + y) = 0: (1) (a) What is the order of this ODE? Is it linear or nonlinear? Separable or non-separable? (b) What can you say from general theorems about the existence and unique- ness of solutions of the initial value problem for (1) with the initial condition y(x0) = y0? (c) Show that the substitution u(x) = x + y(x) gives the following ODE: u0 = cos2 u: (2) (d) Solve (2) for u(x) and write down the general solution of (1) for y(x).
Recommended publications
  • Ordinary Differential Equation 2
    Unit - II Linear Systems- Let us consider a system of first order differential equations of the form dx = F(,)x y dt (1) dy = G(,)x y dt Where t is an independent variable. And x & y are dependent variables. The system (1) is called a linear system if both F(x, y) and G(x, y) are linear in x and y. dx = a ()t x + b ()t y + f() t dt 1 1 1 Also system (1) can be written as (2) dy = a ()t x + b ()t y + f() t dt 2 2 2 ∀ = [ ] Where ai t)( , bi t)( and fi () t i 2,1 are continuous functions on a, b . Homogeneous and Non-Homogeneous Linear Systems- The system (2) is called a homogeneous linear system, if both f1 ( t )and f2 ( t ) are identically zero and if both f1 ( t )and f2 t)( are not equal to zero, then the system (2) is called a non-homogeneous linear system. x = x() t Solution- A pair of functions defined on [a, b]is said to be a solution of (2) if it satisfies y = y() t (2). dx = 4x − y ........ A dt Example- (3) dy = 2x + y ........ B dt dx d 2 x dx From A, y = 4x − putting in B we obtain − 5 + 6x = 0is a 2 nd order differential dt dt 2 dt x = e2t equation. The auxiliary equation is m2 − 5m + 6 = 0 ⇒ m = 3,2 so putting x = e2t in A, we x = e3t obtain y = 2e2t again putting x = e3t in A, we obtain y = e3t .
    [Show full text]
  • 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]
  • Phase Plane Diagrams of Difference Equations
    PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS TANYA DEWLAND, JEROME WESTON, AND RACHEL WEYRENS Abstract. We will be determining qualitative features of a dis- crete dynamical system of homogeneous difference equations with constant coefficients. By creating phase plane diagrams of our system we can visualize these features, such as convergence, equi- librium points, and stability. 1. Introduction Continuous systems are often approximated as discrete processes, mean- ing that we look only at the solutions for positive integer inputs. Dif- ference equations are recurrence relations, and first order difference equations only depend on the previous value. Using difference equa- tions, we can model discrete dynamical systems. The observations we can determine, by analyzing phase plane diagrams of difference equa- tions, are if we are modeling decay or growth, convergence, stability, and equilibrium points. For this paper, we will only consider two dimensional systems. Suppose x(k) and y(k) are two functions of the positive integers that satisfy the following system of difference equations: x(k + 1) = ax(k) + by(k) y(k + 1) = cx(k) + dy(k): This system is more simply written in matrix from z(k + 1) = Az(k); x(k) a b where z(k) = is a column vector and A = . If z(0) = z y(k) c d 0 is the initial condition then the solution to the initial value prob- lem z(k + 1) = Az(k); z(0) = z0; 1 2 TANYA DEWLAND, JEROME WESTON, AND RACHEL WEYRENS is k z(k) = A z0: Each z(k) is represented as a point (x(k); y(k)) in the Euclidean plane.
    [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]
  • Phase Plane Methods
    Chapter 10 Phase Plane Methods Contents 10.1 Planar Autonomous Systems . 680 10.2 Planar Constant Linear Systems . 694 10.3 Planar Almost Linear Systems . 705 10.4 Biological Models . 715 10.5 Mechanical Models . 730 Studied here are planar autonomous systems of differential equations. The topics: Planar Autonomous Systems: Phase Portraits, Stability. Planar Constant Linear Systems: Classification of isolated equilib- ria, Phase portraits. Planar Almost Linear Systems: Phase portraits, Nonlinear classi- fications of equilibria. Biological Models: Predator-prey models, Competition models, Survival of one species, Co-existence, Alligators, doomsday and extinction. Mechanical Models: Nonlinear spring-mass system, Soft and hard springs, Energy conservation, Phase plane and scenes. 680 Phase Plane Methods 10.1 Planar Autonomous Systems A set of two scalar differential equations of the form x0(t) = f(x(t); y(t)); (1) y0(t) = g(x(t); y(t)): is called a planar autonomous system. The term autonomous means self-governing, justified by the absence of the time variable t in the functions f(x; y), g(x; y). ! ! x(t) f(x; y) To obtain the vector form, let ~u(t) = , F~ (x; y) = y(t) g(x; y) and write (1) as the first order vector-matrix system d (2) ~u(t) = F~ (~u(t)): dt It is assumed that f, g are continuously differentiable in some region D in the xy-plane. This assumption makes F~ continuously differentiable in D and guarantees that Picard's existence-uniqueness theorem for initial d ~ value problems applies to the initial value problem dt ~u(t) = F (~u(t)), ~u(0) = ~u0.
    [Show full text]
  • Calculus and Differential Equations II
    Calculus and Differential Equations II MATH 250 B Linear systems of differential equations Linear systems of differential equations Calculus and Differential Equations II Second order autonomous linear systems We are mostly interested with2 × 2 first order autonomous systems of the form x0 = a x + b y y 0 = c x + d y where x and y are functions of t and a, b, c, and d are real constants. Such a system may be re-written in matrix form as d x x a b = M ; M = : dt y y c d The purpose of this section is to classify the dynamics of the solutions of the above system, in terms of the properties of the matrix M. Linear systems of differential equations Calculus and Differential Equations II Existence and uniqueness (general statement) Consider a linear system of the form dY = M(t)Y + F (t); dt where Y and F (t) are n × 1 column vectors, and M(t) is an n × n matrix whose entries may depend on t. Existence and uniqueness theorem: If the entries of the matrix M(t) and of the vector F (t) are continuous on some open interval I containing t0, then the initial value problem dY = M(t)Y + F (t); Y (t ) = Y dt 0 0 has a unique solution on I . In particular, this means that trajectories in the phase space do not cross. Linear systems of differential equations Calculus and Differential Equations II General solution The general solution to Y 0 = M(t)Y + F (t) reads Y (t) = C1 Y1(t) + C2 Y2(t) + ··· + Cn Yn(t) + Yp(t); = U(t) C + Yp(t); where 0 Yp(t) is a particular solution to Y = M(t)Y + F (t).
    [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]
  • TWO DIMENSIONAL FLOWS Lecture 4: Linear and Nonlinear Systems
    TWO DIMENSIONAL FLOWS Lecture 4: Linear and Nonlinear Systems 4. Linear and Nonlinear Systems in 2D In higher dimensions, trajectories have more room to manoeuvre, and hence a wider range of behaviour is possible. 4.1 Linear systems: definitions and examples A 2-dimensional linear system has the form x˙ = ax + by y˙ = cx + dy where a, b, c, d are parameters. Equivalently, in vector notation x˙ = Ax (1) where a b x A = and x = (2) c d ! y ! The Linear property means that if x1 and x2 are solutions, then so is c1x1 + c2x2 for any c1 and c2. The solutions ofx ˙ = Ax can be visualized as trajectories moving on the (x,y) plane, or phase plane. 1 Example 4.1.1 mx¨ + kx = 0 i.e. the simple harmonic oscillator Fig. 4.1.1 The state of the system is characterized by x and v =x ˙ x˙ = v k v˙ = x −m i.e. for each (x,v) we obtain a vector (˙x, v˙) ⇒ vector field on the phase plane. 2 As for a 1-dimensional system, we imagine a fluid flowing steadily on the phase plane with a local velocity given by (˙x, v˙) = (v, ω2x). − Fig. 4.1.2 Trajectory is found by placing an imag- • inary particle or phase point at (x0,v0) and watching how it moves. (x,v) = (0, 0) is a fixed point: • static equilibrium! Trajectories form closed orbits around (0, 0): • oscillations! 3 The phase portrait looks like... Fig. 4.1.3 NB ω2x2 + v2 is constant on each ellipse.
    [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]
  • Phaser: an R Package for Phase Plane Analysis of Autonomous ODE Systems by Michael J
    CONTRIBUTED RESEARCH ARTICLES 43 phaseR: An R Package for Phase Plane Analysis of Autonomous ODE Systems by Michael J. Grayling Abstract When modelling physical systems, analysts will frequently be confronted by differential equations which cannot be solved analytically. In this instance, numerical integration will usually be the only way forward. However, for autonomous systems of ordinary differential equations (ODEs) in one or two dimensions, it is possible to employ an instructive qualitative analysis foregoing this requirement, using so-called phase plane methods. Moreover, this qualitative analysis can even prove to be highly useful for systems that can be solved analytically, or will be solved numerically anyway. The package phaseR allows the user to perform such phase plane analyses: determining the stability of any equilibrium points easily, and producing informative plots. Introduction Repeatedly, when a system of differential equations is written down, it cannot be solved analytically. This is particularly true in the non-linear case, which unfortunately habitually arises when modelling physical systems. As such, it is common that numerical integration is the only way for a modeller to analyse the properties of their system. Consequently, many software packages exist today to assist in this step. In R, for example, the package deSolve (Soetaert et al., 2010) deals with many classes of differential equation. It allows users to solve first-order stiff and non-stiff initial value problem ODEs, as well as stiff and non-stiff delay differential equations (DDEs), and differential algebraic equations (DAEs) up to index 3. Moreover, it can tackle partial differential equations (PDEs) with the assistance ReacTran (Soetaert and Meysman, 2012).
    [Show full text]
  • Phase-Plane Methods
    NPS ARCHIVE 1969 HSIUNG, Y. PHASE-PLANE METHODS Yun-Lan Hsiung DUDLEY KNOX LIBRARY NAVAL POSTGRADUATE SCHOOL MONTEREY, CA 93943-5101 United States Naval Postgraduate School THESIS PHASE-PLANE METHODS by Yun-Lan Hsiung October 1969 Tlnu document, kcud been app\ove.d fan. pubtic *e- leo*e and 6aJU; itb duVubution <u unlimUzd. T13S516 DUDLEY KNOX LIBRARY SCHOOL N M/AL POSTGRADUATE MONTEREY, CA 93943-5101 Phase -Plane Methods by Yun-Lan Hsiung Lieutenant Commander, Chinese Navy B. S., Chinese Naval Academy, 1954 Submitted in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE IN ELECTRICAL ENGINEERING from the NAVAL POSTGRADUATE SCHOOL October 1969 c \ \<o'=\ ABSTRACT The phase-plane method is a graphical method for linear and non-linear second-order systems. There are many techniques for obtaining the phase portrait which consists of a number of phase trajectories on the phase- plane. From the summary and discussion the isocline method is a most general and useful method. The only difficulty is in the labor required. Solution by digi- tal computer reduces this difficulty and makes the iso- cline method a much more useful method in non-linear systems analysis and design applications. LIBRARY IAYAL POSTGRADUATE SCHOOL MOHTJEBJSy, CALIF. 93940 TABLE OF CONTENTS Page I. INTRODUCTION 15 A. CONDITIONS FOR LINEARITY AND DEFINITION 15 OF NONLINEARITY B. CLASSES OF NONLINEARITY 15 C. NONLINEAR ANALYSIS 19 D. THE GRAPHICAL METHOD - THE PHASE-PLANE 20 (PHASE-PLANE ANALYSIS OF NONLINEAR SYSTEMS) II. SUMMARY OF USEFUL GRAPHICAL METHODS FOR 22 CONSTRUCTING PHASE TRAJECTORIES ON THE PHASE-PLANE A.
    [Show full text]