Interpolation forward Gregory-Newton

An exothermic, relatively slow reaction takes place in a reactor under your supervision. Yesterday, after you left the plant, the temperature of the reactor went out of control, for a yet unknown reason, until the operator put it under control by changing the cooling water flow rate. Your supervisor has asked you to prepare a report regarding this incident. As the first step, you must know when the reactor reached its maximum temperature and what was the value of this maximum temperature. A computer was recording the temperature of the reactor at one-hour intervals. These time-temperature data are given in Table E3.1. Write a general MATLAB function for -order one-dimensional interpolation by Gregory-Newton forward interpolation formula to solve this problem. 

In this section, we will develop two interpolation methods for equally spaced data (I) the Gregory-Newton formulas, which are based on forward and backward differences, and (2) Stirling s interpolation formula, based on central differences. 

The Gregory-Newton forward interpolation formula can be derived using the forward finite difference relations derived in Secs. 3.2 and 3.4. Eq. (3.17), written for the function/... 

This is the Gregory-Newton forward interpolation formula. The general formula of the above... 

Method of Solution The function uses the general formula of the Gregory-Newton forward interpolation [Eq. (3.123)] to perform the n-order interpolation. The input to the function specifying the number of base points must be at least (n + 1). 

Solution to Example 3.1. It interpolates the time-temperature data % given in Table E3.1 by Gregory-Newton forward interpolation % formula and finds the maximum temperature and the time this % maximum happened. 

YI = GregoryNewton(X,Y,XI,N) applies the Nth-order % Gregory-Newton forward interpolation to find YI, the % values of the underlying function Y at the points in the vector XI. The vector X specifies the points at which the data Y is given. 

This method is accomplished by first replacing the function y = f x) with a polynomial approximation, such as the Gregory-Newton forward interpolation formula [Eq. (3.122)], In practice, the interval [x, x ] is being divided into several segments, each of width h, and the Gregory-Newton forward interpolation formula becomes note that =x + h) ... 

Gregory-Newton forward interpolation method. Lagrange polynomial interpolation method. Cubic splines interpolation method. 

These values give the best interpolation at t = 575°F rather than the Newton-Gregory s forward or backward interpolation formulae. 

THE PROGRAM USES THE NEWTON-GREGORY FORWARD AND BACKWARD INTERPOLATIONS AND STIRLING S CENTRAL DIFFERENCE METHOD. 

THIS PROGRAM PERFORMS THE NEWTON-GREGORY FORWARD INTERPOLATION OF POLYNOMIALS... 

---