Spline Integration

Another method of integrating unequally spaced data points is to interpolate the data using a suitable interpolation method, such as cubic splines, and then evaluate the integral from the relevant polynomial. Therefore, the integral of Eq. (4.66) may be calculated by integrating Eq. (3.143) over the interval jcJ and summing up these terms for all the intervals  

Prior to calculating the integral from Eq. (4.108), the values of the second derivative at the base points should be calculated from Eq. (3.147). Note that if a natural spline interpolation is employed, the second derivatives for the first and the last intervals are equal to zero. Eq. (4.108) is basically an improved trapezoidal formula in which the value of the integral by trapezoidal mle [the first term in the bracket of Eq. (4.108)] is corrected for the curvature of the function [the second term in the bracket of Eq. (4.108)]. 

---

The Jacobian matrix defined in (3.41) can be easily computed by the same interpolation technique. The idea is to differentiate (3.60) with respect to the parameters changing the order of differentiation and spline integration. 

In ref 147 the authors apply Adomian decomposition method to develop an efficient algorithm of a special second-order ordinary initial value problems. The Adomian decomposition method is known that does not require discretization, so is computer time efficient. The authors are studied the Adomian decomposition method and the results obtained are compared with previously known results using the Quintic C -spline integration methods. 

Expander-compressor shafts are preferably designed to operate below the first lateral critical speed and torsional resonance. A flame-plated band of aluminum alloy or similarly suitable material is generally applied to the shaft in the area sensed by the vibration probes to preclude erroneous electrical runout readings. This technique has been used on hundreds of expanders, steam turbines, and other turbomachines with complete success. Unless integral with the shaft, expander wheels (disks) are often attached to the shaft on a special tapered profile, with dowel-type keys and keyways. The latter design attempts to avoid the stress concentrations occasionally associated with splines and conventional keyways. It also reduces the cost of manufacture. When used, wheels are sometimes secured to the tapered ends of the shaft by a common center stretch rod which is pre-stressed during assembly. This results in a constant preload on each wheel to ensure proper contact between wheels and shaft at the anticipated extremes of temperature and speed. 

A generalized partial differential equation solver which handles simultaneous parabolic, one dimensional elliptic, ordinary and integral equations and uses B-splines with an adaptive grid was written to solve the model. Further details on the model and solution method can be found in Reference 14. 

POLYMATH. AIChE Cache Corp, P O Box 7939, Austin TX 78713-7939. Polynomial and cubic spline curvefitting, multiple linear regression, simultaneous ODEs, simultaneous linear and nonlinear algebraic equations, matrix manipulations, integration and differentiation of tabular data by way of curve fit of the data. 

R = 0 leads back to the problem of interpolation by spline functions. It should be noted at this point that the condition stated by eq. (4) is not sufficient for the construction of calibration curves and additional considerations have to take effect. A reformulation of the problem stated in Equations (2) and (3) gives us with 8y. = 1 for all i calibration points another look at the problem that clarifies the role of the integral in Equation (2) as balanced against a value of R. Find S (x) to... 

Differentiation of the experimental concentration-time curve would then need interpolation or smoothing, e.g.,by using splines. Parallelization in a typical robotic environment is easy when using the integral method with a few or even only one single well for characterization of one enzyme variant. 

M64 Function value, derivatives and definite integral of a cubic spline at a given point 6400 6450... 

REM J FUNCTION VALUE, DERIVATIVES AND DEFINITE t 4484 REN t INTEGRAL OF A CUBIC SPLINE AT A GIVEN -POINT t 4484 REN lltltllllltmillltlllltllltltllllllltltlllllttltll 4488 REN INPUT ... 

In addition, we are interested in functions that are at least twice continuously differentiable. One can draw several such curves satisfying (4.27), and the "smoothest" of them is the one minimizing the integral (4.19). It can be shown that the solution of this constrained minimization problem is a natural cubic spline (ref. 12). We call it smoothing spline. 

The input is similar to that of the module 1463. No end condition flag is used since only natural splines can be fitted. On the other hand, you should specify the maximum number IM of iterations. The module returns the array S defined in the description of the module M63, and hence the function value, the derivatives and the integral at a specified X can be computed by calling the module 1464. The important additional inputs needed by the module 1465 are the standard errors given in the vector D. With all D(I) = 0, the module... 

Since spline interpolation and integration is mucht faster than solving the sensitivity equations and the original differential equations, the direct method is superior to the indirect one in terms of numerical efficiency, whenever it is feasible. 

In spite of its simplicity the direct integral method has relatively good statistical properties and it may be even superior to the traditional indirect approach in ill-conditioned estimation problems (ref. 18). Good performance, however, can be expected only if the sampling is sufficiently dense and the measurement errors are moderate, since otherwise spline interpolation may lead to severely biased estimates. 

A. Yermakova, S. Vajda and P. Valkd, Direct integral method via spline approximation for estimating rate constants. Applied Catalysis,... 

Data array calculations with cubic spline fitting and numerical integration appears to yield data of much greater accuracy. 

Packages exist that use various discretizations in the spatial direction and an integration routine in the time variable. PDECOL uses B-splines for the spatial direction and various GEAR methods in time (Ref 247). PDEPACK and DSS (Ref 247) use finite differences in the spatial direction and GEARB in time (Ref 66). REACOL (Ref. 106) uses orthogonal collocation in the radial direction and LSODE in the axial direction, while REACFD uses finite difference in the radial direction both codes are restricted to modeling chemical reactors. 

---