Find a General Solution to the Cauchy-euler Equation Order 3
Uncertain Input Data Problems and the Worst Scenario Method
Ivan Hlaváěek Dr. , ... Ivo Babuška Dr. , in North-Holland Series in Applied Mathematics and Mechanics, 2004
Starshaped domains
This particular family of domains enables us to calculate explicitly the constants emerging in the estimates we are dealing with. In our setting, domains Ω1, Ω2, and Ω3 have the following properties: Ω1 is a domain starshaped with respect to the origin of the coordinate system,
(23.23)
where α > 1 is a given constant, and . The domain Ω2 can be N-unstable.
The domains Ω1 and Ω3 are related by the mapping ϰ(x) = αx, which maps Ω1 onto Ω3. Therefore, ϰ can be used to transform functions defined on Ω1 into functions defined on Ω3; namely u 1α (y) = u 1α (ϰ(x)) = u 1(x) and G 1α (y) = G 1α (ϰ(x)) = G 1(x), where .
By defining an auxiliary equation for u 1α and by comparing u 3 and u 1 α, we can, after some labor, infer the estimate.
Theorem 23.3
Let ∇G ∈ [L ∞(Ω13)]2, G ∈ H 2(Ω13) and let α ∈ (1, α 0]. Then
where C is a bounded positive parameter. Also,
with the same parameter C.
Proof. See (Babuška and Chleboun, 2002, Section 4). ◇
Remark 23.9
The formula for C can be evaluated so that the parameter can be calculated. To form an idea of what C looks like, let us specify that it includes , and α 0. □
Read full chapter
URL:
https://www.sciencedirect.com/science/article/pii/S0167593104800147
Ordinary differential equations
Brent J. Lewis , ... Andrew A. Prudil , in Advanced Mathematics for Engineering Students, 2022
2.2.3 Euler–Cauchy equation
The Euler–Cauchy equations of the form
(2.52)
can also be solved by an algebraic method by substituting into Eq. (2.52), yielding
.
Omitting (where ) yields the auxiliary equation
(2.53)
Similarly, three cases arise for the general solution of Eq. (2.52), as shown in Table 2.3.
Table 2.3. General solution of the Euler–Cauchy equation.
Case | Roots (Eq. (2.53) ) | Basis (Eq. (2.52) ) | General solution |
---|---|---|---|
I | Distinct real m 1,m 2 | ||
II | Real double root | ||
III | Complex conjugate | ||
m 1 =μ +iν | |||
m 2 =μ −iν |
Example 2.2.6
(Case I) Solve .
Solution. The auxiliary equation is . Hence, and from Eq. (2.53). Therefore, the general solution is . [answer]
Example 2.2.7
(Case II) Solve .
Solution. The auxiliary equation has a double root . Therefore, the general solution is . [answer]
Example 2.2.8
(Case III) Solve .
Solution. The auxiliary equation is , which implies . The general solution is . [answer]
Read full chapter
URL:
https://www.sciencedirect.com/science/article/pii/B9780128236819000101
Ordinary Differential Equations
Xin-She Yang , in Engineering Mathematics with Examples and Applications, 2017
12.3 Second-Order Equations
For second-order ordinary differential equations (ODEs), it is generally more tricky to find their general solutions. However, a special case with significantly practical importance and mathematical simplicity is the second-order linear differential equation with constant coefficients in the following form
(12.16)
where the coefficients b and c are constants, and is a known function of x. Obviously, the more general form is
(12.17)
However, if we divide both sides by a, we will reach our standard form. Here we assume . In a special case of , it reduces to a first-order linear differential equation which has been discussed in the previous section. So we will start our discussion from (12.16).
A differential equation is said to be homogeneous if . For a given generic second-order differential equation (12.16), a function that satisfies the homogeneous equation
(12.18)
is called the complementary function, denoted by . Obviously, the complementary function alone cannot satisfy the original equation (12.16) because there is no way to produce the required on the right-hand side. Therefore, we have to find a specific function, called the particular integral, so that it indeed satisfies the original equation (12.16). The combined general solution
(12.19)
will automatically satisfy the original equation (12.16). The general solution of (12.16) consists of two parts: the complementary function and the particular integral . We can obtain these two parts separately, and simply add them together because the original equation is linear, so their solutions are linear combinations.
First things first, how to obtain the complementary function? The general technique is to assume that it takes the form
(12.20)
where A is a constant, and λ is an exponent to be determined. Substituting this assumed form into the homogeneous equation and using both and , we have
(12.21)
Since should not be zero (otherwise, we have a trivial solution everywhere), we can divide all the terms by , and we have
(12.22)
which is the characteristic equation for the homogeneous equation. It is also called the auxiliary equation of the ODE. The solution of λ in this case is simply
(12.23)
For simplicity, we can take as it does not affect the results.
From quadratic equations, we know that there are three possibilities for λ. They are: I) two real distinct roots, II) two identical roots, and III) two complex roots.
In the case of two different roots: . Then, both and satisfy the homogeneous equation, so their linear combination forms the complementary function
(12.24)
where A and B are constants.
In the special case of identical roots , or
(12.25)
we cannot simply write
(12.26)
because it still only one part of the complementary function where is just another constant. In this case, we should try a different combination, say, to see if it satisfies the homogeneous equation or not. Since , and , we have
(12.27)
where we have used (identical roots) and (the auxiliary equation). This indeed implies that also satisfies the homogeneous equation. Therefore, the complementary function for the identical roots is
(12.28)
Now let us use this technique to solve a second-order homogeneous ODE.
Example 12.4
The second-order homogeneous equation
has a corresponding auxiliary equation
It has two real roots
So the complementary function is
But for the differential equation
its auxiliary equation becomes
which has two identical roots . The complementary function in this case can be written as
As complex roots always come in pairs, the case of complex roots would give
(12.29)
where α and β are real numbers. The complementary function becomes
(12.30)
where we have used the Euler formula and also absorbed the constants A and B into and .
A special case is when , so the roots are purely imaginary. We have , and . Equation (12.16) in this case becomes
(12.31)
which is a differential equation for harmonic motions such as the oscillations of a pendulum or a small-amplitude seismic detector. Here β is the angular frequency of the system.
Example 12.5
For a simple pendulum of mass m shown in Fig. 12.2, we now try to derive its equation of oscillations and its period.
Since the motion is circular, the tension or the centripetal force T is thus given by
where is the angular velocity.
Forces must be balanced both vertically and horizontally. The component of T in the vertical direction is which must be equivalent to mg, though in the opposite direction. Here g is the acceleration due to gravity. That is
Since θ is small or , we have . This means that .
In the horizontal direction, Newton's second law implies that the horizontal force must be equal to the mass m times the acceleration . Now we have
Dividing both sides by mL, we have
Since θ is small, we have . Therefore, we finally have
This is the equation of motion for a simple pendulum. From equation (12.31), we know that the angular frequency is or . Thus the period of the pendulum is
(12.32)
We can see that the period is independent of the bob mass. For m and , the period is approximately seconds.Up to now, we have found only the complementary function. Now we will try to find the particular integral for the original non-homogeneous equation (12.16). For particular integrals, we do not intend to find the general form; any specific function or integral that satisfies the original equation (12.16) will do. Before we can determine the particular integral, we have to use some trial functions, and such functions will have strong similarity to function . For example, if is a polynomial such as , we will try a similar form and try to determine the coefficients. Let us demonstrate this by an example.
Example 12.6
In order to solve the differential equation
we first find its complementary function. From the earlier example, we know that the complementary function can be written as
For the particular integral, we know that , so we try the form
Thus and . Substituting them into the original equation, we have
or
As this equality must be true for any x, so the coefficients of the same power of x on both sides of the equation should be equal. That is
which gives , and . So the general solution becomes
Similarly, if , we will try to so as to determine a. In addition, for or , we will attempt the general form .
Again let us demonstrate this through an example.
Example 12.7
The motion of a damped pendulum is governed by
Its auxiliary equation becomes
which has two complex solutions
Therefore, the complementary function becomes
where A and B are two undetermined constants.
Now we try to get the particular integral. Since , so we try the similar form
Since
we have
which leads to
which requires
or
This gives
The particular solution is
Therefore, the general solution becomes
As time is sufficiently long (i.e., ), we can see that the first two terms will approach zero (since ). Therefore, the solution will be dominated by the last two terms. That is, the long-term behavior becomes
In an earlier chapter about trigonometry, we have shown that the following identity holds
where
Therefore, we have
where and .
Read full chapter
URL:
https://www.sciencedirect.com/science/article/pii/B9780128097304000161
Differential equations and difference equations
Mary Attenborough , in Mathematics for Electrical Engineering and Computing, 2003
Example 14.8
Solve the differential equation 5y″ + 6y′ +5y = 6 cos(t), where y (0) = 0 and y′ (0) = 0.
Step 1
To find the complementary function we solve the homogeneous equation 5y″ + 6y′ + 5y = 0. Trying solutions of the form y = A e λt leads to the auxiliary equation 5λ 2 + 6λ + 5 = 0. Notice that a quick way to get the auxiliary equation is to 'replace' y″ by λ2, y′ by A, and y by 1. The auxiliary equation has solutions
Comparing this with λ = k ± jω0 gives k = −0.6, ω0 = 0.8. This means that the complementary function has the form
(given as Case (2) in the general method), with k = −0.6 and ω0 = 0.8. This gives
Step 2
As f(t) = 6 cos(t), from Table 14.1 we decide to try a particular solution of the form y = c cos(t) + d sin(t). Then
and
Substituting in
we find
As we want this to be an identity we equate the coefficients of cos(t) and the coefficients of sin (t) and get the two equations
Hence, a particular solution is y = sin(t).
Step 3
The general solution is given by the sum of the complementary function and a particular solution, so
Step 4
Use the given initial conditions to find values for the constants A and B. Substituting y = 0 when t = 0 into
then
To use the other condition, that y′ = 0 when t = 0, we need to differentiate the general solution to find an expression for y'. Differentiating (using the product rule):
and using y′ = 0 when t = 0 gives
We have already found A = 0 from the first condition, so
Therefore, the solution is
Check
Substitute into
giving
which is true for all values of t.
We can also check that the solution satisfies the given initial conditions.
Read full chapter
URL:
https://www.sciencedirect.com/science/article/pii/B9780750658553500401
Vibrations and the Harmonic Oscillator
J.E. House , in Fundamentals of Quantum Mechanics (Third Edition), 2018
6.2 Linear Differential Equations with Constant Coefficients
Because differential equations are required when representing motion, some familiarity with that branch of mathematics is beneficial when studying quantum mechanics. This section provides a very brief introduction to such mathematics. However, in this section, the results of several important theorems in differential equations will be presented. Because of the nature of this book, the theorems will be presented in an operational manner and used without proof. The interested reader should consult a text on differential equations for more details.
A linear differential equation with constant coefficients is an equation of the form
(6.7)
in which the constants a 0(x), a 1(x), …, and F(x) have values that change only with changes in x. A particularly important equation of this type is the second-order case:
(6.8)
The differential operator D is defined as
(6.9)
When an operator meets the conditions represented as
(6.10)
and that which can be shown by the relationship
(6.11)
the operator is called a linear operator.
From Eq. (6.8), it can be seen that a second-order linear differential equation can be written in the operator notation as
(6.12)
The solution of an equation of this form is obtained by considering an auxiliary equation , which is obtained by writing an equation in the form
(6.13)
when the general differential equation is written as
(6.14)
The auxiliary equation is called the complementary equation, and its solution is known as the complementary solution. The complete solution of the differential equation is the sum of the particular and general solution of the complementary equation. These principles will be illustrated by working through the following example:
Suppose it is necessary to find the general solution of
(6.15)
In operator form, this equation can be written as
(6.16)
A solution of this type of equation is frequently of the form
(6.17)
in which a and b are to be determined by the solutions of the complementary equation,
(6.18)
Therefore, factoring the polynomial gives
(6.19)
from which it can be seen that m = 4 and m = 1. In this case, the general solution of Eq. (6.15) is
(6.20)
It can easily be verified that this is the solution by using it in the complementary equation. If this is the solution, then
and the second derivative can be written as
Therefore, the auxiliary equation becomes
(6.21)
By substituting for the first and second derivatives, it is found that
(6.22)
By simplification, it can be seen that this equation reduces to 0 = 0. However, it can also be shown that a particular solution is
(6.23)
It will now be shown that this particular solution also satisfies the general equation. From Eq. (6.23) it is found that
(6.24)
When these values are substituted in Eq. (6.15) the result is
Therefore, the complete solution of Eq. (6.15) is the sum of the two expressions,
(6.25)
In most problems, it is sufficient to obtain a general solution, and "singular" solutions that do not describe the physical behavior of the system are ignored.
It must be mentioned that there are two arbitrary constants that characterize the solution that was obtained. Of course, the solution of an nth order equation would result in n constants. In quantum mechanics, these constants are determined by the physical constraints of the system (known as boundary conditions), as was shown in Chapter 3.
The equation
(6.26)
has the auxiliary equation that can be written as
(6.27)
From this equation, it can be seen that m 2 = − 1 and m = ± i. Therefore, the general solution is
(6.28)
At this point, it is useful to remember that
(6.29)
and
(6.30)
From these relationships, it is evident that
(6.31)
and the solution y = sin x satisfies the equation. In fact, if it is assumed that a solution to Eq. (6.26) is of the form
(6.32)
then
(6.33)
and
(6.34)
Therefore,
(6.35)
Therefore, Eq. (6.35) shows the solution represented by Eq. (6.32) satisfies Eq. (6.26). A differential equation can have only one general solution, so the solution shown in Eqs. (6.28) and (6.32) must be equal:
(6.36)
It can be seen that when x = 0, C 1 + C 2 = B. Differentiating Eq. (6.36) gives
(6.37)
However, when x = 0, sin x = 0. It is apparent, then, that
(6.38)
By substituting for A and B and simplifying, the result can be shown as
(6.39)
If C 2 = 0 and C 1 = 1, then
(6.40)
Therefore, if C 2 = 1 and C 1 = 0, the result can be shown as
(6.41)
The relationships shown in Eqs. (6.40) and (6.41) are known as Euler's formulas.
Suppose it is necessary to solve the differential equation
(6.42)
When written in operator form the result is
(6.43)
Therefore, the auxiliary equation is written as
(6.44)
and its roots are found by using the quadratic formula,
(6.45)
Therefore, the solution of Eq. (6.42) is
(6.46)
By expanding the exponentials, this equation can also be written as
(6.47)
Using Euler's formulas to express the exponential functions in terms of sin and cos, the result is
(6.48)
In general, if the auxiliary equation has roots a ± bi, the solution of the differential equation has the form
(6.49)
An equation having the form
(6.50)
was obtained in solving the particle in the one-dimensional box problem in quantum mechanics (see Chapter 3). The auxiliary equation is
(6.51)
and it has the solutions m = ia. Therefore, the solution of Eq. (6.50) can be written as
(6.52)
This is exactly the form of the solution found when solving the particle in a one-dimensional box. The boundary conditions of a particular system make it possible to evaluate the constants in the solution, as was shown in Chapter 3.
Read full chapter
URL:
https://www.sciencedirect.com/science/article/pii/B9780128092422000061
Differential equations
Huw Fox , Bill Bolton , in Mathematics for Engineers and Technologists, 2002
Exceptional cases of particular integrals
There are situations when the obvious form of function to be tried to obtain the particular integral yields no result because when it is substituted in the differential equation we obtain 0 = 0. This occurs when the right-hand side of the non-homogeneous differential equation consists of a function that is also a term in the complementary function. To illustrate this, consider the differential equation:
The complementary function is y = A e−2x + B e x . For the particular integral, if we try the solution y = A e kx we obtain:
and so no solution for A. In such cases we have to try something different.
Thus we try y = Ax e kx . This gives, for the above differential equation:
Thus k = −2 and −2A + 4Ax − 2A + A − 2Ax − 2Ax = 1. Equating constants gives 3A = 1, equating the x coefficients gives 0 = 0, and so the particular integral is:
Example
Determine the general solution of the differential equation:
The corresponding homogeneous differential equation is:
Trying y = A e sx as a solution gives the auxiliary equation:
This can be factored as (s − 5)(s + 2) = 0 and so the complementary function is y c = A e5x + B e−2x . The right-hand side of the non-homogeneous differential equation is the sum of two terms for which the trial functions would be C and Dx e kx . We thus try the sum of these. Thus:
Equating exponential terms gives k = −2, 4Dx − 2D − 2D + 6Dx − 3D − 10Dx = −1 and so D = 1/7. Equating constants gives −10C = 4 and so C = −4/10. Thus the particular integral is −(4/10) + (1/7) e2x . The general solution is therefore:
Read full chapter
URL:
https://www.sciencedirect.com/science/article/pii/B9780750655446500060
Advances in Atomic, Molecular, and Optical Physics
Erik Torrontegui , ... Juan Gonzalo Muga , in Advances In Atomic, Molecular, and Optical Physics, 2013
2.1 Invariant-Based Inverse Engineering
Lewis-Riesenfeld invariants.The Lewis and Riesenfeld (1969) theory is applicable to a quantum system that evolves with a time-dependent Hermitian Hamiltonian , which supports a Hermitian dynamical invariant satisfying
(1)
Therefore its expectation values for an arbitrary solution of the time-dependent Schrödinger equation do not depend on time. can be used to expand as a superposition of "dynamical modes" ,
(2)
where are time-independent amplitudes, and are orthonormal eigenvectors of the invariant ,
(3)
The are real constants, and the Lewis-Riesenfeld phases are defined as (Lewis and Riesenfeld, 1969)
(4)
We use for simplicity a notation for a discrete spectrum of but the generalization to a continuum or mixed spectrum is straightforward. We also assume a non-degenerate spectrum. Non-Hermitian invariants and Hamiltonians have been considered for example in Gao et al. (1991, 1992), Lohe (2009), Ibáñez et al. (2011).
Inverse engineering.Suppose that we want to drive the system from an initial Hamiltonian to a final one , in such a way that the populations in the initial and final instantaneous bases are the same, but admitting transitions at intermediate times. To inverse engineer a time-dependent Hamiltonian and achieve this goal, we may first define the invariant through its eigenvalues and eigenvectors. The Lewis-Riesenfeld phases may also be chosen as arbitrary functions to write down the time-dependent unitary evolution operator
(5)
obeys ,where the dot means time derivative. Solving formally this equation for , we get
(6)
According to Eq. (6), for a given invariant there are many possible Hamiltonians corresponding to different choices of phase functions . In general does not commute with , so the eigenstates of , do not coincide with the eigenstates of . does not necessarily commute with either. If we impose and , the eigenstates will coincide, which guarantees a state transfer without final excitations. In typical applications the Hamiltonians and are given, and set the initial and final configurations of the external parameters. Then we define and its eigenvectors accordingly, so that the commutation relations are obeyed at the boundary times and, finally, is designed via Eq. (6). While the may be taken as fully free time-dependent phases in principle, they may also be constrained by a pre-imposed or assumed structure of . Sections 3–5 present examples of how this works for expansions, transport, and internal state control.
A generalization of this inverse method for non-Hermitian Hamiltonians was considered in Ibáñez et al. (2011). Inverse engineering was applied to accelerate the slow expansion of a classical particle in a time-dependent harmonic oscillator without final excitation. This system may be treated formally as a quantum two-level system with non-Hermitian Hamiltonian (Gao et al., 1991,1992).
Quadratic in momentum invariants. Lewis and Riesenfeld (1969) paid special attention to the time-dependent harmonic oscillator and its invariants quadratic in position and momentum. Later on Lewis and Leach (1982) found, in the framework of classical mechanics, the general form of the Hamiltonian compatible with quadratic-in-momentum invariants, which includes non-harmonic potentials. This work, and the corresponding quantum results of Dhara and Lawande (1984), constitutes the basis of this section.
A one-dimensional Hamiltonian with a quadratic-in-momentum invariant must have the form , 3 with the potential (Lewis and Leach, 1982; Dhara and Lawande, 1984)
(7)
, and are arbitrary functions of time that satisfy the auxiliary equations
(8)
(9)
where is a constant. Their physical interpretation will be explained below and depends on the operation. A quadratic-in- dynamical invariant is given, up to a constant factor, by
(10)
Now in Eq. (4) satisfies (Lewis and Riesenfeld, 1969; Dhara and Lawande, 1984)
(11)
and the function can be written as (Dhara and Lawande, 1984)
(12)
in terms of the solution (normalized in -space) of the auxiliary Schrödinger equation
(13)
The strategy of invariant-based inverse engineering here is to design and first so that and commute at initial and final times, except for launching or stopping atoms as in Torrontegui et al. (2011). Then is deduced from Eq. (7). Applications will be discussed in Sections 3 and 4.
Read full chapter
URL:
https://www.sciencedirect.com/science/article/pii/B9780124080904000025
Higher-Order Differential Equations
Martha L. Abell , James P. Braselton , in Differential Equations with Mathematica (Fourth Edition), 2016
4.6.1 Second-Order Cauchy-Euler Equations
Consider the second-order homogeneous Cauchy-Euler equation
(4.22)
where a≠0. Notice that because the coefficient of y″ is zero if x = 0, we must restrict our domain to either x > 0 or x < 0 in order to ensure that the theory of second-order equations stated in Section 4.1 holds.
Suppose that y = x m , x > 0, for some constant m. Substitution of y = x m with derivatives y′ = mx m−1 and y″ = m(m − 1)x m−2 into equation (4.22) yields
Then, y = x m is a solution of equation (4.22) if m satisfies
(4.23)
which is called the characteristic equation(or auxiliary equation ) associated with the Cauchy-Euler equation of order two. The solutions of the characteristic equation completely determine the general solution of the homogeneous Cauchy-Euler equation of order two. Let m 1 and m 2 denote the two solutions of the characteristic (or auxiliary) equation (4.23):
Hence, we can obtain two real roots, one repeated real root, or a complex conjugate pair depending on the values of a, b, and c. We state a general solution that corresponds to the different types of roots.
Theorem 8
Let m 1 and m 2 be the solutions of equation (4.23).
- 1.
-
If m 1≠m 2 are real and distinct, two linearly independent solutions of equation (4.22) are and ; a general solution of (4.22) is
- 2.
-
If m 1 = m 2, two linearly independent solutions of equation (4.22) are and ; a general solution of (4.22) is
- 3.
-
If m 1, 2 = α ± βi, β≠0, two linearly independent solutions of equation (4.22) are and ; a general solution of (4.22) is
Example 4.6.1
Solve each of the following equations: (a) 3x 2 y″ − 2xy′ + 2y = 0, x > 0; (b) x 2 y″ − xy′ + y = 0, x > 0; (c) x 2 y″ − 5xy′ + 10y = 0, x > 0.
Solution
(a) If y = x m , y′ = mx m−1, and y″ = m(m − 1)x m−2, substitution into the differential equation yields
Hence, the auxiliary equation is
with roots m 1 = 2/3 and m 2 = 1. Therefore, a general solution is y = c 1 x 2/3 + c 2 x. We obtain the same results with DSolve. Entering
Clear[x, y]gensol=DSolve[3x ∧ 2y"[x]−2xy′[x]+2y[x]==0, y[x], x]
{{y[x] → x 2/3 C[1] + xC[2]}}
finds a general solution of the equation, naming the result gensol, and then entering
toplot=Table[gensol[[1, 1, 2]]/.{C[1]→i, C[2]→j}, {i, −2, 2, 2}, {j, −2, 2, 2}];
p1=Plot[toplot, {x, 0, 12}, PlotRange→{−6, 6}, AspectRatio→1, AxesLabel→{x, y}, PlotLabel→"(a)"]
defines toplot to be the list of functions obtained by replacing c[1] in gensol[[1,1,2]] by − 2, 0, and 2 and C[2] in gensol[[1,1,2]] by − 2, 0, and 2, and graphs the set of functions toplot on the interval [0, 12]. See Figure 4-26 (a).
(b) In this case, the auxiliary equation is
with root m = 1 of multiplicity 2. Hence, a general solution is . As in the previous example, we see that we obtain the same results with DSolve. See Figure 4-26 (b).
Clear[x, y]gensol=DSolve[x ∧ 2y"[x]−xy′[x]+y[x]==0, y[x], x]
{{y[x] → xC[1] + xC[2]Log[x]}}
toplot=Table[gensol[[1, 1, 2]]/.{C[1]→i, C[2]→j}, {i, −2, 2, 2}, {j, −2, 2, 2}];
p2=Plot[toplot, {x, 0, 10}, PlotRange→{−5, 5}, AspectRatio→1, AxesLabel→{x, y}, PlotLabel→"(b)"]
(c) The auxiliary (characteristic) equation is given by
with complex conjugate roots . Thus, a general solution is .
Again, we see that we obtain equivalent results with DSolve. First, we find a general solution of the equation, naming the resulting output gensol.
Clear[x, y]gensol=DSolve[x ∧ 2y"[x]−5xy′[x]+10y[x]==0, y[x], x]
{{y[x] → x 3 C[2]Cos[Log[x]] + x 3 C[1]Sin[Log[x]]}}
Now, we define y(x) to be the general solution obtained in gensol. (The same result is obtained with Part by entering y[x_]=gensol[[1,1,2]].)
y[x_]=x 3 C[2]Cos[Log[x]]+x 3 C[1]Sin[Log[x]]
x 3 C[2]Cos[Log[x]] + x 3 C[1]Sin[Log[x]]
To find the values of C[1] and C[2] so that the solution satisfies the initial conditions y(1) = a and y′(1) = b, we use Solve and name the resulting list cvals.
cvals=Solve[{y[1]==a, y′[1]==b}, {C[1], C[2]}]
{{C[1] →−3a + b, C[2] → a}}
The solution to the initial-value problem
is obtained by replacing C[1] and C[2] in y(x) by the values found in cvals.
y[x_]=y[x]/.cvals[[1]]
ax 3Cos[Log[x]] + (−3a + b)x 3Sin[Log[x]]
Note that when you enter the following Plot command, Mathematica may display several error messages because each solution is undefined if x = 0. Nevertheless, the resulting graphs are displayed correctly.
This solution is then graphed for various initial conditions in Figure 4-26 (c).
toplot=Table[y[x], {a, −2, 2, 2}, {b, −2, 2, 2}];
p3=Plot[toplot, {x, 0, 2}, PlotRange→All, AspectRatio→1, AxesLabel→{x, y}, PlotLabel→"(c)"]
Show[GraphicsRow[{p1, p2, p3}]]
Read full chapter
URL:
https://www.sciencedirect.com/science/article/pii/B9780128047767000048
Comment on: The (G′/G)-expansion method for the nonlinear lattice equations [Commun Nonlinear Sci Numer Simulat 17 (2012) 3490–3498]
İsmail Aslan , in Communications in Nonlinear Science and Numerical Simulation, 2012
Now, we make the following observations:
- (i)
-
Eqs. (1) and (2) are special cases of Eq. (4). Namely, taking a = α, b = −α, and c = β in Eq. (4) leads to Eq. (1); taking a = −α, b = 0, and c = 1 in Eq. (4) leads to Eq. (2). Hence, Eq. (4) includes both Eq. (1) and Eq. (2).
- (ii)
-
The basic (G′/G)-expansion method uses the auxiliary equation where λ and μ are arbitrary constants. Recently, Aslan [5] demonstrated the redundancy of the parameter λ. In other words, one can make the assumption λ = 0 without loss of generality. This approach reduces the number of the parameters at the outset without affecting the generality of the results. Also, it should be clear that (5) is a further generalization of (3) in the meaning that the sum goes from l = −m to l = m instead of from l = 0 to l = m where m is a positive integer. This fact indicates that the solutions obtained by (5) include the solutions obtained by (3).
- (iii)
-
It is also worth to mention here that the so-called two-component Volterra lattice equations
(6)
Read full article
URL:
https://www.sciencedirect.com/science/article/pii/S1007570412001773
Analytic solutions of a second order iterative functional differential equation
Pingping Zhang , Lufang Mi , in Applied Mathematics and Computation, 2009
In this section, we prove the existence of analytic solution to Eq. (1.1).
Theorem 3.1
Under one of the conditions in Theorems 2.1–2.3, Eq. (1.5) has an analytic solution of the form in a neighborhood of the number , where is an analytic solution of (1.7) .
Proof
By Theorems 2.1–2.3, we may find an analytic solution of the auxiliary equation (1.7) in the form of (2.1) such that . Clearly the inverse exists and is analytic in a neighborhood of . If we define . Then
as required. The proof is completed. □
Read full article
URL:
https://www.sciencedirect.com/science/article/pii/S009630030800920X
Find a General Solution to the Cauchy-euler Equation Order 3
Source: https://www.sciencedirect.com/topics/mathematics/auxiliary-equation