Part 2: Ordinary Differential Equations (Odes) (This Is New Material, See Kreyszig, Chapters 1-6, and Related Numerics in Chaps

Total Page:16

File Type:pdf, Size:1020Kb

Part 2: Ordinary Differential Equations (Odes) (This Is New Material, See Kreyszig, Chapters 1-6, and Related Numerics in Chaps New Mexico Tech Hyd 510 Hydrology Program Quantitative Methods in Hydrology Part 2: Ordinary Differential Equations (ODEs) (This is new material, see Kreyszig, Chapters 1-6, and related numerics in Chaps. 19, 20, 20.1-20.7, and 21.1-21.3.) Fundamentals of Differential Equations The calculus problems we’ve reviewed have mostly been involved with finding the numerical value of one magnitude or another. For example, we’ve sought the location x0 and value f(x0) of the maxima of a function, the locations x where f(x)=0 (roots of f(x)), the slope (derivative) of f(x) at x, b or the value of a definite integral f (x)dx . ∫a But now we are going to solve engineering and science problems where the unknown is itself a function, that is, we are going to develop and solve models that express “the dependence of certain variables on others” (Petrovskii, 1956). What kind of questions can these models answer? Here are some examples from hydrology. How will drawdown in an aquifer change with distance from a pumping well? How will the discharge from a lake, flowing over a weir, vary in time? How does the velocity profile in the atmospheric boundary layer vary with distance above the ground surface? How does moisture content vary with depth? In these cases the dependent (unknown) variables are, respectfully, drawdown [L]1, discharge [L3], velocity [L/t], and moisture content [L3/ L3/]=[-]. The respective independent (known) variables are distance [L], time [t], elevation [L], and depth [L]. The problem of finding these functions is most often addressed by solving differential equations, that is, equations in which not only the unknown function(s) occurs, but also its derivatives of various orders. Examples: dh 1. + αh = αh + P dt 0 d 2φ 1 1 Q ' 2. − φ = − φ − w dx 2 λ λ 0 T ∂C ∂C ∂ 2C 3. + v − D = 0 ∂t ∂x ∂x 2 In the first example the unknown is denoted by the letter h, which represents the water level in a lake, and the independent variable by the letter t, which represents time. As there is one independent variable, it is an ordinary differential equation (ODE). As there is only one first order derivative it 1The notation [L] for length and its like (e.g., [M] for mass or [t] for time) is used to indicate units or dimensions of a variable, similar to using standard units like m (meter), s (second), or kg (kilogram). So velocity units can be indicated by [L/t] or [ms-1]. Among other things we use units to keep track of dimensional homogeneity, such that units balance in all terms on each side of an equation. See next page. - 58 - New Mexico Tech Hyd 510 Hydrology Program Quantitative Methods in Hydrology 2 is a first order equation (see below). The other symbols (α, h0, P) denote parameters with no dependence on the unknown, h. Mathematically, its possible that all three vary with the independent variable, t, but in hydrologic practice only h0 and P typically vary with time. In the second example the unknown is denoted by the Greek letter φ , representing hydraulic head in a confined aquifer, and the independent variable by the letter x, representing location. As there is still just one independent variable, it is also an ordinary differential equation (ODE). However, as a second derivative of φ is the highest order derivative it is a second order equation. The symbols λ, φ0, and Qw’ denote parameters that are known. In the third example the unknown is denoted by the letter C, representing solute concentration in a flowing river, and there are two independent variables, denoted by the letters x (space) and t (time). As there are now two (or more) independent variables, it is a partial differential equation (PDE). The highest order derivative in x defines it as a second order equation in space, while the highest order derivative in t defines it as a first order equation in time. (Later we’ll see that this particular type of PDE is called mixed hyperbolic-parabolic PDE.) The symbols v and D denote parameters, the mean river velocity and a “dispersion coefficient.” In each equation the left side contains the unknown and its derivatives, while the right side contains known loadings or forcings. For example, in the first equation the forcing is precipitation, P [L/t]. In the second equation it is pumping Qw’ [L/t] and head in a vertically adjacent aquifer, φ0 [L]. There is no forcing in the third equation (later we’ll see that forcing for that equation comes from the boundary and initial conditions instead). We often rewrite different equations so that they fit this template: unknowns and their derivatives on the left, forcings on the right. As we shall see below, mathematicians take this convention a step further. In short, the first equation is a first order ODE, the second equation is a second order ODE, and the third equation is a second order in space – first order in time PDE. In this section we concentration only on problems with one independent variable, that is, on ODEs. Dimensional homogeneity Each of these DEs is (and must be) dimensionally homogeneous. We can take advantage of this to help detect errors and to improve solutions. For example, the second equation is d 2φ 1 1 Q ' − φ = − φ − w dx 2 λ λ 0 T Here head φ and distance x have units of length [L], parameter λ [L2] includes (=TB’/K’) the influence of aquifer transmissivity (T [L2/t]) and aquitard leakance (B’/K’ [Lt/L]), where B’ is aquitard thickness [L] and K’ is aquitard vertical conductivity[L/t], and Qw’ [L/t] is pumping per unit area. All terms in the second equation have units of [1/L]. Order of equation (text, §1.1) An ordinary differential equation of order n, in one unknown function y, is a relation of the form F[x, y(x), y' (x), y' ' (x),...., y n (x)] = 0 Between the (one) independent variable x and the quantitites 2 The physical meaning of parameters like these will become known to you later. - 59 - New Mexico Tech Hyd 510 Hydrology Program Quantitative Methods in Hydrology y(x), y' (x), y' ' (x),...., y n (x) The order of the ODE is the order to the highest derivatve of the unknown function appearing in the differential equation. We’ll focus on 1st and 2nd order ODEs, respectively Chaps 1 and 2 of the textbook. In our applications, time is usually the independent variable in 1st order ODEs, while it is almost always space in 2nd order ODEs. The major exception to these observations is in geophysics, involving seismic, electromagnetic, and other signals, where 2nd order ODEs in time are also encountered. Operator We can write an equation in the form of an operator L( ), where L(unknown)= forcing or L( )= forcing. In the examples above, the operators are: d( ) 1. L( ) = + α( ) dt d 2 ( ) 1 2. L( ) = − ( ) dx 2 λ ∂( ) ∂( ) ∂ 2 ( ) 3. L( ) = + v − D ∂t ∂x ∂x 2 An operator can be an algebraic, differential or integral operator. Our examples only demonstrate differential operators. Linear dependence There are tremendous advantages of a linear model in deriving and applying solutions. For example, if we have a linear model, then we can scale or “convolute” the model result to represent a new forcing, that is, a new solution, without having to resolve the equation. If we have a linear model we can use certain solution methods that take advantage of linearity, such as LaPlace or Fourier transform methods. Even if the model is nonlinear, if the nonlinearity is mild we might be able to linearize (approximate) the model and still take advantage of linearity. On the other hand there are certain problems that have a sufficient non-linearity that it must be addressed directly. Consider chaotic behavior, as in climate models, which is the result of one kind of non-linearity that cannot be ignored. In vadose zone hydrology the multiphase flow equations are highly non-linear. The land-surface energy balance involves highly non-linear radiation terms. Taking advantage of the operator notation, there are three conditions necessary for an operator (and its differential equation) to be linear. Letting x and y represent unknowns, and a and b constants, these are (a) L(ax) = aL(x) - 60 - New Mexico Tech Hyd 510 Hydrology Program Quantitative Methods in Hydrology (b) L(x + y) = L(x) + L( y) (c) L(ax + by) = aL(x) + bL( y) As you can see the last condition covers the other two. Let’s take a look at two algebraic examples. Example a: bx+c =0 has the operator L( ) = b·( ), so that L(x)=-c. Show that L meets the requirements for a linear operator. Example b: ax2+ bx+c =0 has the operator L( ) = a·( )2 + b·( ), so that L(x)=-c. Show that L does not meets the requirements for a linear operator. That is, this is a non-linear algebraic equation. Aside: if b=0, show that the operator remains non-linear in x, but that a transform of variables to z=x2 leads to a new problem az+c =0 which is linear in z. We often look for such simple transformations to (exactly) linearize a problem (such transformations are common in vadose zone hydrology and infiltration models used for hillslope or watershed hydrology, and for simple models of phreatic aquifers through the so-called linearized Dupuit Approximation). Now, let’s examine operators for differential equations.
Recommended publications
  • Modelling for Aquifer Parameters Determination
    Journal of Engineering Science, Vol. 14, 71–99, 2018 Mathematical Implications of Hydrogeological Conditions: Modelling for Aquifer Parameters Determination Sheriffat Taiwo Okusaga1, Mustapha Adewale Usman1* and Niyi Olaonipekun Adebisi2 1Department of Mathematical Sciences, Olabisi Onabanjo University, Ago-Iwoye, Nigeria 2Department of Earth Sciences, Olabisi Onabanjo University, Ago-Iwoye, Nigeria *Corresponding author: [email protected] Published online: 15 April 2018 To cite this article: Sheriffat Taiwo Okusaga, Mustapha Adewale Usman and Niyi Olaonipekun Adebisi. (2018). Mathematical implications of hydrogeological conditions: Modelling for aquifer parameters determination. Journal of Engineering Science, 14: 71–99, https://doi.org/10.21315/jes2018.14.6. To link to this article: https://doi.org/10.21315/jes2018.14.6 Abstract: Formulations of natural hydrogeological phenomena have been derived, from experimentation and observation of Darcy fundamental law. Mathematical methods are also been applied to expand on these formulations, but the management of complexities which relate to subsurface heterogeneity is yet to be developed into better models. This study employs a thorough approach to modelling flow in an aquifer (a porous medium) with phenomenological parameters such as Transmissibility and Storage Coefficient on the basis of mathematical arguments. An overview of the essential components of mathematical background and a basic working knowledge of groundwater flow is presented. Equations obtained through the modelling process
    [Show full text]
  • The Cauchy-Kovalevskaya Theorem
    CHAPTER 2 The Cauchy-Kovalevskaya Theorem This chapter deals with the only “general theorem” which can be extended from the theory of ODEs, the Cauchy-Kovalevskaya Theorem. This also will allow us to introduce the notion of non-characteristic data, principal symbol and the basic clas- sification of PDEs. It will also allow us to explore why analyticity is not the proper regularity for studying PDEs most of the time. Acknowledgements. Some parts of this chapter – in particular the four proofs of the Cauchy-Kovalevskaya theorem for ODEs – are strongly inspired from the nice lecture notes of Bruce Diver available on his webpage, and some other parts – in particular the organisation of the proof in the PDE case and the discussion of the counter-examples – are strongly inspired from the nice lectures notes of Tsogtgerel Gantumur available on his webpage. 1. Recalls on analyticity Let us first do some recalls on the notion of real analyticity. Definition 1.1. A function f defined on some open subset U of the real line is said to be real analytic at a point x if it is an infinitely di↵erentiable function 0 2U such that the Taylor series at the point x0 1 f (n)(x ) (1.1) T (x)= 0 (x x )n n! − 0 n=0 X converges to f(x) for x in a neighborhood of x0. A function f is real analytic on an open set of the real line if it is real analytic U at any point x0 . The set of all real analytic functions on a given open set is often denoted by2C!U( ).
    [Show full text]
  • Bessel's Differential Equation Mathematical Physics
    R. I. Badran Bessel's Differential Equation Mathematical Physics Bessel's Differential Equation: 2 Recalling the DE y n y 0 which has a sinusoidal solution (i.e. sin nx and cos nx) and knowing that these solutions can be treated as power series, we can find a solution to the Bessel's DE which is written as: 2 2 2 x y xy (x p )y 0 . The solution is represented by a series. This series very much look like a damped sine or cosine. It is called a Bessel function. To solve the Bessel' DE, we apply the Frobenius method by mn assuming a series solution of the form y an x . n0 Substitute this solution into the DE equation and after some mathematical steps you may find that the indicial equation is 2 2 m p 0 where m= p. Also you may get a1=0 and hence all odd a's are zero as well. The recursion relation can be found as a a n n2 (n m 2)2 p 2 with (n=0, 1, 2, 3, ……etc.). Case (1): m= p (seeking the first solution of Bessel DE) We get the solution 2 4 p x x y a x 1 2(2p 2) 2 4(2p 2)(2p 4) This can be rewritten as: x p2n J (x) y a (1)n p 2n n0 2 n!( p n)! 1 Put a 2 p p! The Bessel function has the factorial form x p2n J (x) (1)n p 2n p , n0 n!( p n)!2 R.
    [Show full text]
  • Review of the Methods of Transition from Partial to Ordinary Differential
    Archives of Computational Methods in Engineering https://doi.org/10.1007/s11831-021-09550-5 REVIEW ARTICLE Review of the Methods of Transition from Partial to Ordinary Diferential Equations: From Macro‑ to Nano‑structural Dynamics J. Awrejcewicz1 · V. A. Krysko‑Jr.1 · L. A. Kalutsky2 · M. V. Zhigalov2 · V. A. Krysko2 Received: 23 October 2020 / Accepted: 10 January 2021 © The Author(s) 2021 Abstract This review/research paper deals with the reduction of nonlinear partial diferential equations governing the dynamic behavior of structural mechanical members with emphasis put on theoretical aspects of the applied methods and signal processing. Owing to the rapid development of technology, materials science and in particular micro/nano mechanical systems, there is a need not only to revise approaches to mathematical modeling of structural nonlinear vibrations, but also to choose/pro- pose novel (extended) theoretically based methods and hence, motivating development of numerical algorithms, to get the authentic, reliable, validated and accurate solutions to complex mathematical models derived (nonlinear PDEs). The review introduces the reader to traditional approaches with a broad spectrum of the Fourier-type methods, Galerkin-type methods, Kantorovich–Vlasov methods, variational methods, variational iteration methods, as well as the methods of Vaindiner and Agranovskii–Baglai–Smirnov. While some of them are well known and applied by computational and engineering-oriented community, attention is paid to important (from our point of view) but not widely known and used classical approaches. In addition, the considerations are supported by the most popular and frequently employed algorithms and direct numerical schemes based on the fnite element method (FEM) and fnite diference method (FDM) to validate results obtained.
    [Show full text]
  • The Standard Form of a Differential Equation
    Fall 06 The Standard Form of a Differential Equation Format required to solve a differential equation or a system of differential equations using one of the command-line differential equation solvers such as rkfixed, Rkadapt, Radau, Stiffb, Stiffr or Bulstoer. For a numerical routine to solve a differential equation (DE), we must somehow pass the differential equation as an argument to the solver routine. A standard form for all DEs will allow us to do this. Basic idea: get rid of any second, third, fourth, etc. derivatives that appear, leaving only first derivatives. Example 1 2 d d 5 yx()+3 ⋅ yx()−5yx ⋅ () 4x⋅ 2 dx dx This DE contains a second derivative. How do we write a second derivative as a first derivative? A second derivative is a first derivative of a first derivative. d2 d d yx() yx() 2 dx dxxd Step 1: STANDARDIZATION d Let's define two functions y0(x) and y1(x) as y0()x yx() and y1()x y0()x dx Then this differential equation can be written as d 5 y1()x +3y ⋅ 1()x −5y ⋅ 0()x 4x dx This DE contains 2 functions instead of one, but there is a strong relationship between these two functions d y1()x y0()x dx So, the original DE is now a system of two DEs, d d 5 y1()x y0()x and y1()x +3y ⋅ 1()x −5y ⋅ 0()x 4x⋅ dx dx The convention is to write these equations with the derivatives alone on the left-hand side. d y0()x y1()x dx This is the first step in the standardization process.
    [Show full text]
  • Method of Separation of Variables
    2.2. Linearity 33 2.2 Linearity As in the study of ordinary differential equations, the concept of linearity will be very important for liS. A linear operator L lJy definition satisfies Chapter 2 (2.2.1) for any two functions '/1,1 and 112, where Cl and ('2 are arhitrary constants. a/at and a 2 /ax' are examples of linear operators since they satisfy (2.2.1): Method of Separation iJ at (el'/11 + C2'/1,) of Variables a' 2 (C1'/11 +C2U 2) ax It can be shown (see Exercise 2.2.1) that any linear combination of linear operators is a linear operator. Thus, the heat operator 2.1 Introduction a iJl In Chapter I we developed from physical principles an lmderstanding of the heat at - k a2·2 eqnation and its corresponding initial and boundary conditions. We are ready is also a linear operator. to pursue the mathematical solution of some typical problems involving partial A linear equation fora it is of the form differential equations. \Ve \-vilt use a technique called the method of separation of variables. You will have to become an expert in this method, and so we will discuss L(u) = J, (2.2.2) quite a fev.; examples. v~,fe will emphasize problem solving techniques, but \ve must also understand how not to misuse the technique. where L is a linear operator and f is known. Examples of linear partial dijjerentinl A relatively simple but typical, problem for the equation of heat conduction equations are occurs for a one - dimensional rod (0 ::; x .::; L) when all the thermal coefficients are constant.
    [Show full text]
  • CORSO ESTIVO DI MATEMATICA Differential Equations Of
    CORSO ESTIVO DI MATEMATICA Differential Equations of Mathematical Physics G. Sweers http://go.to/sweers or http://fa.its.tudelft.nl/ sweers ∼ Perugia, July 28 - August 29, 2003 It is better to have failed and tried, To kick the groom and kiss the bride, Than not to try and stand aside, Sparing the coal as well as the guide. John O’Mill ii Contents 1Frommodelstodifferential equations 1 1.1Laundryonaline........................... 1 1.1.1 Alinearmodel........................ 1 1.1.2 Anonlinearmodel...................... 3 1.1.3 Comparingbothmodels................... 4 1.2Flowthroughareaandmore2d................... 5 1.3Problemsinvolvingtime....................... 11 1.3.1 Waveequation........................ 11 1.3.2 Heatequation......................... 12 1.4 Differentialequationsfromcalculusofvariations......... 15 1.5 Mathematical solutions are not always physically relevant . 19 2 Spaces, Traces and Imbeddings 23 2.1Functionspaces............................ 23 2.1.1 Hölderspaces......................... 23 2.1.2 Sobolevspaces........................ 24 2.2Restrictingandextending...................... 29 2.3Traces................................. 34 1,p 2.4 Zero trace and W0 (Ω) ....................... 36 2.5Gagliardo,Nirenberg,SobolevandMorrey............. 38 3 Some new and old solution methods I 43 3.1Directmethodsinthecalculusofvariations............ 43 3.2 Solutions in flavours......................... 48 3.3PreliminariesforCauchy-Kowalevski................ 52 3.3.1 Ordinary differentialequations............... 52 3.3.2 Partial differentialequations...............
    [Show full text]
  • Differential Equations and Linear Algebra
    Chapter 1 First Order Equations 1.1 Four Examples : Linear versus Nonlinear A first order differential equation connects a function y.t/ to its derivative dy=dt. That rate of change in y is decided by y itself (and possibly also by the time t). Here are four examples. Example 1 is the most important differential equation of all. dy dy dy dy 1/ y 2/ y 3/ 2ty 4/ y2 dt D dt D dt D dt D Those examples illustrate three linear differential equations (1, 2, and 3) and a nonlinear differential equation. The unknown function y.t/ is squared in Example 4. The derivative y or y or 2ty is proportional to the function y in Examples 1, 2, 3. The graph of dy=dt versus y becomes a parabola in Example 4, because of y2. It is true that t multiplies y in Example 3. That equation is still linear in y and dy=dt. It has a variable coefficient 2t, changing with time. Examples 1 and 2 have constant coefficient (the coefficients of y are 1 and 1). Solutions to the Four Examples We can write down a solution to each example. This will be one solution but it is not the complete solution, because each equation has a family of solutions. Eventually there will be a constant C in the complete solution. This number C is decided by the starting value of y at t 0, exactly as in ordinary integration. The integral of f.t/ solves the simplest differentialD equation of all, with y.0/ C : D dy t 5/ f.t/ The complete solution is y.t/ f.s/ds C : dt D D C Z0 1 2 Chapter 1.
    [Show full text]
  • Hydrogeology and Groundwater Flow
    Hydrogeology and Groundwater Flow Hydrogeology (hydro- meaning water, and -geology meaning the study of rocks) is the part of hydrology that deals with the distribution and movement of groundwater in the soil and rocks of the Earth's crust, (commonly in aquifers). The term geohydrology is often used interchangeably. Some make the minor distinction between a hydrologist or engineer applying themselves to geology (geohydrology), and a geologist applying themselves to hydrology (hydrogeology). Hydrogeology (like most earth sciences) is an interdisciplinary subject; it can be difficult to account fully for the chemical, physical, biological and even legal interactions between soil, water, nature and society. Although the basic principles of hydrogeology are very intuitive (e.g., water flows "downhill"), the study of their interaction can be quite complex. Taking into account the interplay of the different facets of a multi-component system often requires knowledge in several diverse fields at both the experimental and theoretical levels. This being said, the following is a more traditional (reductionist viewpoint) introduction to the methods and nomenclature of saturated subsurface hydrology, or simply hydrogeology. © 2014 All Star Training, Inc. 1 Hydrogeology in Relation to Other Fields Hydrogeology, as stated above, is a branch of the earth sciences dealing with the flow of water through aquifers and other shallow porous media (typically less than 450 m or 1,500 ft below the land surface.) The very shallow flow of water in the subsurface (the upper 3 m or 10 ft) is pertinent to the fields of soil science, agriculture and civil engineering, as well as to hydrogeology.
    [Show full text]
  • The Calculus of Trigonometric Functions the Calculus of Trigonometric Functions – a Guide for Teachers (Years 11–12)
    1 Supporting Australian Mathematics Project 2 3 4 5 6 12 11 7 10 8 9 A guide for teachers – Years 11 and 12 Calculus: Module 15 The calculus of trigonometric functions The calculus of trigonometric functions – A guide for teachers (Years 11–12) Principal author: Peter Brown, University of NSW Dr Michael Evans, AMSI Associate Professor David Hunt, University of NSW Dr Daniel Mathews, Monash University Editor: Dr Jane Pitkethly, La Trobe University Illustrations and web design: Catherine Tan, Michael Shaw Full bibliographic details are available from Education Services Australia. Published by Education Services Australia PO Box 177 Carlton South Vic 3053 Australia Tel: (03) 9207 9600 Fax: (03) 9910 9800 Email: [email protected] Website: www.esa.edu.au © 2013 Education Services Australia Ltd, except where indicated otherwise. You may copy, distribute and adapt this material free of charge for non-commercial educational purposes, provided you retain all copyright notices and acknowledgements. This publication is funded by the Australian Government Department of Education, Employment and Workplace Relations. Supporting Australian Mathematics Project Australian Mathematical Sciences Institute Building 161 The University of Melbourne VIC 3010 Email: [email protected] Website: www.amsi.org.au Assumed knowledge ..................................... 4 Motivation ........................................... 4 Content ............................................. 5 Review of radian measure ................................. 5 An important limit ....................................
    [Show full text]
  • MTHE / MATH 237 Differential Equations for Engineering Science
    i Queen’s University Mathematics and Engineering and Mathematics and Statistics MTHE / MATH 237 Differential Equations for Engineering Science Supplemental Course Notes Serdar Y¨uksel November 26, 2015 ii This document is a collection of supplemental lecture notes used for MTHE / MATH 237: Differential Equations for Engineering Science. Serdar Y¨uksel Contents 1 Introduction to Differential Equations ................................................... .......... 1 1.1 Introduction.................................... ............................................ 1 1.2 Classification of DifferentialEquations ............ ............................................. 1 1.2.1 OrdinaryDifferentialEquations ................. ........................................ 2 1.2.2 PartialDifferentialEquations .................. ......................................... 2 1.2.3 HomogeneousDifferentialEquations .............. ...................................... 2 1.2.4 N-thorderDifferentialEquations................ ........................................ 2 1.2.5 LinearDifferentialEquations ................... ........................................ 2 1.3 SolutionsofDifferentialequations ................ ............................................. 3 1.4 DirectionFields................................. ............................................ 3 1.5 Fundamental Questions on First-Order Differential Equations ...................................... 4 2 First-Order Ordinary Differential Equations ..................................................
    [Show full text]
  • Linear Differential Equations Math 130 Linear Algebra
    where A and B are arbitrary constants. Both of these equations are homogeneous linear differential equations with constant coefficients. An nth-order linear differential equation is of the form Linear Differential Equations n 00 0 Math 130 Linear Algebra anf (t) + ··· + a2f (t) + a1f (t) + a0f(t) = b: D Joyce, Fall 2013 that is, some linear combination of the derivatives It's important to know some of the applications up through the nth derivative is equal to b. In gen- of linear algebra, and one of those applications is eral, the coefficients an; : : : ; a1; a0, and b can be any to homogeneous linear differential equations with functions of t, but when they're all scalars (either constant coefficients. real or complex numbers), then it's a linear differ- ential equation with constant coefficients. When b What is a differential equation? It's an equa- is 0, then it's homogeneous. tion where the unknown is a function and the equa- For the rest of this discussion, we'll only consider tion is a statement about how the derivatives of the at homogeneous linear differential equations with function are related. constant coefficients. The two example equations 0 Perhaps the most important differential equation are such; the first being f −kf = 0, and the second 00 is the exponential differential equation. That's the being f + f = 0. one that says a quantity is proportional to its rate of change. If we let t be the independent variable, f(t) What do they have to do with linear alge- the quantity that depends on t, and k the constant bra? Their solutions form vector spaces, and the of proportionality, then the differential equation is dimension of the vector space is the order n of the equation.
    [Show full text]