Linear Equations of Higher Order

The singular solution, p = 0, which is dv/dt = 0, so that y/x = constant is a solution. This solution cannot satisfy the two boundary conditions. 

The most general linear differential equation of n order can be written in the standard form 

The most general solution to Eq. 2.167 is called the homogeneous or complementary solution. The notation complementary comes about when f x) is 

We first focus our efforts in solving the unforced or homogeneous equation, and then concentrate on dealing with solutions arising when forcing is applied through fix) (i.e., the particular solution). 

It is clear in the homogeneous Eq. 2.167 that if all coefficients a, ... a f x) were zero, then we could solve the final equation by n successive integrations of 

---

As stated earlier, much is known about linear equations of higher order, but no general technique is available to solve the nonlinear equations that arise frequently in natural and man-made systems. When analysis fails to uncover the analytical solution, the last recourse is to undertake numerical solution methods, as introduced in Chapters 7 and 8. 

Linear Equations of Higher Order 65 dy/dx with r and d y/dx with so that... 

Linear Equations of Higher Order 73 1. Method of Undetermined Coefficients... 

Linear Equations of Higher Order 87 These are now separable, so that within an arbitrary constant ... 

To the same order of approximation of the equations, that is, with only terms linear in / (v) kept, better approximations to the viscosity may be found by considering the equations of higher order than Eqs. (1-86) and (1-87). These new equations will, to this order of approximation, have zero on the left sides (since the higher order coefficients are taken equal to zero) on the right sides appears the factor (p/fii) multiplied by a series of terms like those in Eq. (1-110). Using these equations, and the first order terms of Eq. (1-86) for arbitrary v,... 

As a general rule it is more difficult to solve differential equations of higher orders than the first. Of these, the linear equation is the most important. A linear equation of the nth order is one in which the dependent variable and its n derivatives are all of the first degree and are not multiplied together. If a higher power appears the equation is not linear, and its solution is, in general, more difficult to find. The typical form is... 

The steady-state kinetic treatment of random reactions is complex and gives rise to rate equations of higher order in substrate and product terms. For kinetic treatment of random reactions that display the Michaelis-Menten (i.e. hyperbolic velocity-substrate relationship) or linear (linearly transformed kinetic plots) kinetic behavior, the quasi-equilibrium assumption is commonly made to analyze enzyme kinetic data. 

The scaling rate, i.e. the mass of scale deposited per unit time was also influenced by supersaturation ratio [57]. The higher the supersaturation ratio the more scale forming components are readily available for the scale to form. Hence it was expected that at higher supersaturation ratios the scaling rates were faster. In this single pipe flow experiment, this relationship has been confirmed [1] as a linear equation of first order (R = 0.9932) between mass of scale deposited per unit time (kgs ) and supersaturation level. Figure 7 illustiates this relationship. 

If the water is deeper than 200 m, the linear long wave equation should be applied. For the region shallower than 200 m, the shallow water theory with a term for bottom friction included should be used. This shallow water theory includes the first order approximation of the amplitude dependent dispersion. Under special conditions, the term for frequency dependent dispersion should be included. If the purpose of the simulation is to determine the runup height, the equations of higher order approximations are not necessary. 

If temis of higher order than linear in t are neglected, the transverse magnetization evolves in the presence of the first bipolar gradient pulse according to (equation Bl.14.2 and equation B 1.14.61 ... 

Assuming that at the initial instant the angular velocity dp/dt — 0, we conclude that the mass m, placed at any point around the point of equilibrium, remains at rest. Of course, it is only an approximation, because we preserved in the power series, (Equation (3.146)), only the linear term and discarded terms of higher orders. Formally, this case is characterized by infinitely large period of free vibrations... 

As done previously, in The Newton-Raphson Algorithm (p.48), we neglect all but the first two terms in the expansion. This leaves us with an approximation that is not very accurate but, since it is a linear equation, is easy to deal with. Algorithms that include additional higher terms in the Taylor expansion, often result in fewer iterations but require longer computation times due to the calculation of higher order derivatives. 

A linear PDF, such as the diffusion equation, can only have first-order and zero-order source or sink rates. But what if the source or sink term is of higher order An example would be the generalized reaction... 

The First-Order Linear Inhomogeneous Differential Equation (FOLIDE) First-Order Reaction Including Back Reaction Reaction of Higher Order Catalyzed Reactions... 

Just as a is the linear polarizability, the higher order terms p and y (equation 19) are the first and second hvperpolarizabilities. respectively. If the valence electrons are localized and can be assigned to specific bonds, the second-order coefficient, 6, is referred to as the bond (hyper) polarizability. If the valence electron distribution is delocalized, as in organic aromatic or acetylenic molecules, 6 can be described in terms of molecular (hyper)polarizability. Equation 19 describes polarization at the atomic or molecular level where first-order (a), second-order (6), etc., coefficients are defined in terms of atom, bond, or molecular polarizabilities, p is then the net bond or molecular polarization. 

The equation (386) is not closed again and produces new Green functions of higher order. And so on. These sequence of equations can not be closed in the general case and should be truncated at some point. Below we consider some possible approximations. The other important point is, that average populations and lesser Green functions should be calculated self-consistently. In equilibrium (linear response) these functions are easy related to the spectral functions. But at finite voltage they should be calculated independently. 

The example of this reaction demonstrates another important facet of kinetics. Figure 5.4 shows side by side the experimental data plotted as a first order-first order reversible reaction and as an irreversible reaction of order 1.5. Over a limited conversion range (here about two thirds of the way to equilibrium) the second plot is linear within the scatter of the data points. Although evaluation of the full conversion range leaves no doubt that the reaction is indeed reversible and first order-first order, its rate up to a rather high conversion is approximated surprisingly well by the equation for an irreversible reaction of higher order, in this instance of order 1.5 ... 

We have proved, therefore, that the effect of higher-order corrections is to renormalize the friction coefficient, while the standard Fokker-Planck form is left completely unchanged. In the linear case this is, therefore, an exact equation. It will be shown later on (Ferrario et al.. Chapter VI) that nonlinearity destroys this standard foim. This will lead us to the following conclusions ... 

Although the linear model is the model most commonly encountered in analytical science, not all relationships between a pair of variables can be adequately described by linear regression. A calibration curve does not have to approximate a straight line to be of practical value. The use of higher-order equations to model the association between dependent and independent variables may be more appropriate. The most popular function to model non-linear data and include curvature in the graph is to fit a power-series polynomial of the form... 

Higher order ODEs (of order n) were converted to a system of n coupled linear first order ODEs in section 2.1.4. This system was then solved using the exponential matrix developed earlier. This approach yields analytical solutions for linear ODEs of any order. In section 2.1.5, the given system of coupled linear ODEs was converted to Laplace domain. The resulting linear system of algebraic equations was then solved for the solution in the Laplace domain. The solution obtained in the Laplace domain was then converted to the time domain. 

The analysis now follows that in the previous section, at least qualitatively. First we derive a set of linearized equations of motion and boundary conditions by neglecting all terms that are 0(e2) or higher order. As far as the equations of motion are concerned, it is enough to derive them in one of the two fluids. We choose fluid 1. The linearized O(s)] Navier-Stokes equations are thus... 

---