Solving for

Subroutine REGRES. REGRES is the main subroutine responsible for performing the regression. It solves for the parameters in nonlinear models where all the measured variables are subject to error and are related by one or two constraints. It uses subroutines FUNG, FUNDR, SUMSQ, and SYMINV. 

SOLVES FOR THE PARAMETERS IN NON-LINEAR MEASURED VARIABLES ARE SUBJECT TC ERROR ONE OR TWO CONSTRAINTS. 

A step-limited Newton-Raphson iteration, applied to the Rachford-Rice objective function, is used to solve for A, the vapor to feed mole ratio, for an isothermal flash. For an adiabatic flash, an enthalpy balance is included in a two-dimensional Newton-Raphson iteration to yield both A and T. Details are given in Chapter 7. 

To calculate the vapor load for a single column of a sequence, start by assuming a feed condition such that q can be fixed. Initially assume saturated liquid feed (i.e., q = 1). Equation (5.1) can be written for all NC components of the feed and solved for the necessary values of 0. There are (JVC - 1) real positive values of 0 which satisfy Eq. (5.1), and each lies between the a values of the... 

The Hartree approximation is usefid as an illustrative tool, but it is not a very accurate approximation. A significant deficiency of the Hartree wavefiinction is that it does not reflect the anti-synnnetric nature of the electrons as required by the Pauli principle [7], Moreover, the Hartree equation is difficult to solve. The Hamiltonian is orbitally dependent because the siumnation in equation Al.3.11 does not include the th orbital. This means that if there are M electrons, then M Hamiltonians must be considered and equation A1.3.11 solved for each orbital. 

Applying Flartree-Fock wavefiinctions to condensed matter systems is not routine. The resulting Flartree-Fock equations are usually too complex to be solved for extended systems. It has been argried drat many-body wavefunction approaches to the condensed matter or large molecular systems do not represent a reasonable approach to the electronic structure problem of extended systems. 

The alternative simulation approaches are based on molecular dynamics calculations. This is conceptually simpler that the Monte Carlo method the equations of motion are solved for a system of A molecules, and periodic boundary conditions are again imposed. This method pennits both the equilibrium and transport properties of the system to be evaluated, essentially by numerically solvmg the equations of motion... 

Surface waves at an interface between two innniscible fluids involve effects due to gravity (g) and surface tension (a) forces. (In this section, o denotes surface tension and a denotes the stress tensor. The two should not be coiifiised with one another.) In a hydrodynamic approach, the interface is treated as a sharp boundary and the two bulk phases as incompressible. The Navier-Stokes equations for the two bulk phases (balance of macroscopic forces is the mgredient) along with the boundary condition at the interface (surface tension o enters here) are solved for possible hamionic oscillations of the interface of the fomi, exp [-(iu + s)t + i V-.r], where m is the frequency, is the damping coefficient, s tlie 2-d wavevector of the periodic oscillation and. ra 2-d vector parallel to the surface. For a liquid-vapour interface which we consider, away from the critical point, the vapour density is negligible compared to the liquid density and one obtains the hydrodynamic dispersion relation for surface waves + s>tf. The temi gq in the dispersion relation arises from... 

Often in numerical calculations we detennine solutions g (R) that solve the Scln-odinger equations but do not satisfy the asymptotic boundary condition in (A3.11.65). To solve for S, we rewrite equation (A3.11.65) and its derivative with respect to R in the more general fomi ... 

To solve for the Y, we begin by solving a reference problem wherein the coupling matrix is assumed diagonal with constant couplings within each step. (These could be accomplished by diagonalizing U, but it would be better to avoid this work and use the diagonal U matrix elements.) Then, in temis of the reference U (which we call Uj), we have... 

The polarization P is given in tenns of E by the constitutive relation of the material. For the present discussion, we assume that the polarization P r) depends only on the field E evaluated at the same position r. This is the so-called dipole approximation. In later discussions, however, we will consider, in some specific cases, the contribution of a polarization that has a non-local spatial dependence on the optical field. Once we have augmented the system of equation B 1.5.16. equation B 1.5.17. equation B 1.5.18. equation B 1.5.19 and equation B 1.5.20 with the constitutive relation for the dependence of Pon E, we may solve for the radiation fields. This relation is generally characterized tlirough the use of linear and nonlinear susceptibility tensors, the subject to which we now turn. 

Recall that L contains the frequency or (equation (B2.4.8)). To trace out a spectrum, equation (B2.4.11)) is solved for each frequency. In order to obtain the observed signal v, the sum of the two individual magnetizations can be written as the dot product of two vectors, equation (B2.4.12)). 

Assume that there are two distinct potentials (aside from an additive constant that simply shifts the zero of total energy) V r) and V r) which, when used in //and //, respectively, to solve for a ground... 

CISD Yes/No A/ transformed integrals n N to solve for one Cl energy and eigenvector... 

CAS-MCSCF Yes/Yes A/ transformed integrals to solve for Cl energy many iterations also i... 

In this seetion, we briefly review the basie elements of DFT and the EDA. We then foeus on improvements suggested to remedy some of the shorteomings of the EDA (see seetion B3.2.3.1). A wide variety of teelmiques based on DFT have been developed to ealeulate the eleetron density. Many approaehes do not ealeulate the density direetly but rather solve for either a set of single-eleetron orbitals, or the Green s fiinetion, from whieh tire density is derived. 

The conceptually simplest approach to solve for the -matrix elements is to require the wavefimction to have the fonn of equation (B3.4.4). supplemented by a bound function which vanishes in the asymptote [32, 33, 34 and 35] This approach is analogous to the fiill configuration-mteraction (Cl) expansion in electronic structure calculations, except that now one is expanding the nuclear wavefimction. While successfiti for intennediate size problems, the resulting matrices are not very sparse because of the use of multiple coordinate systems, so that this type of method is prohibitively expensive for diatom-diatom reactions at high energies. 

This approach has one key advantage [30]. Although when solving for x one needs to invert (E -ft +iVr)... 

More sophisticated approaches to describe double layer interactions have been developed more recently. Using cell models, the full Poisson-Boltzmann equation can be solved for ordered stmctures. The approach by Alexander et al shows how the effective colloidal particle charge saturates when the bare particle charge is increased [4o]. Using integral equation methods, the behaviour of the primitive model has been studied, in which all the interactions between the colloidal macro-ions and the small ions are addressed (see, for instance, [44, 45]). 

Thus, the equation we are going to solve for the zeroth-order nuclear motion is... 

---