Gauss Interpolation Formula

-

Gauss Interpolation Formula :


             The general Newton formula or Newton's Divided Difference Interpolation Formula gives 

 y = f(x) = f(x₀)+(x-x₀)f(x₀,x₁)+(x-x₀)                                         (x-x₁)f(x₀,x₁,x₂)+(x-x₀)(x-x₁)                 (x-x₂)f(x₀,x₁,x₂,x₃)+(x-x₀)(x-x₁)(x-x₂)           (x-x₃)f(x₀,x₁,x₂,x₃,x₄)+(x-x₀)(x-x₁)(x-x₂)(x-x₃)(x-x₄)f(x₀,x₁,x₂,x₃,x₄,x₅)+(x-x₀)(x-x₁) (x-x₂)(x-x₃)(x-x₄)                                                              ×(x-x₅)f(x₀,x₁,x₂,x₃,x₄,x₅,x₆)+....

Now putting x₀=x₀ , x₁= x₀+h ,x₂= x₀-h,             x₃= x₀+2h ,x₄=x₀-2h ,x₅=x₀+3h ,x₆=x₀-3h

We immediately get 
f(x) = f(x₀)+(x-x₀)f(x₀,x₀+h)+(x-x₀)(x-x₀-h)

           f(x₀,x₀+h,x₀-h)+(x-x₀)(x-x₀-h)(x-x₀+h)

         f(x₀,x₀+h,x₀-h,x₀+2h)+(x-x₀)(x-x₀-h)                 (x-x₀+2h)f(x₀,x₀+h,x₀-h,x₀+2h,x₀-2h)

     + (x-x₀)(x-x₀-h)(x-x₀+h)(x-x₀-h)(x-x₀+2h)

       f(x₀,x₀+h,x₀-h,x₀+2h,x₀-2h,x₀+3h)+....
                                    ............(1)
Further , let u= x-x₀ /h ⇒x-x₀ = hu

when 

 f(x) = f(x₀)+hu f(x₀,x₀+h)

            + hu(hu-h) f(x₀-h,x₀+h

      +hu(hu-h)(hu+h) f(x₀-h,x₀,x₀+h,x₀+2h)

   +hu(hu-h)(hu+h)(hu-2h) 

                  f(x₀-2h,x₀-h,x₀,x₀+h,x₀+2h)

   +hu(hu-h)(hu+h)(hu-2h)(hu+2h)

        f(x₀-2h,x₀-h,x₀,x₀+h,x₀+2h,x₀+3h)+......

But , we also know that 

   f(x₀,x₀+h) = Δy₀/h 

 f(x₀-h,x₀,x₀+h) =Δ²y₋₁/2h²

f(x₀-h,x₀x₀+h,x₀+2h) = Δ³y₋₁/3!h³

f(x₀-2h,x₀-h,x₀,x₀+h,x₀+2h)= Δ⁴y₋₂/4!h⁴

f(x₀-2h,x₀-h,x₀,x₀+h,x₀+2h,x₀+3h)=

                                               Δ⁵y₋₂/5!h⁵ etc.

Now ,substituting these values in equation(1) , we readily obtain 

y= y₀+hu Δy₀/h + h²u(u-1) Δ²y₋₁/2h² +

      h³u(u-1)(u+1) Δ³y₋₁/3!h³ 

      + h⁴u(u-1(u+1)(u-2) Δ⁴y₋₂/4!h⁴

     + h⁵u(u-1)(u+1)(u-2)(u+2) Δ⁵y₋₂/5!h⁵+.....

⇒y= y₀+uc₁Δy₀ +uc₂Δ²y₋₁+(u+1)c₃Δ³y₋₁

                 +(u+1)c₄ Δ⁴y₋₂ +..........

This result is known as Gauss Forward Formula for equal intervals.


Comments

Post a Comment

If Any Doubt Ask Me

Popular posts from this blog

Accuracy Of Numbers

-
Image
              Accuracy Of Numbers 1.  Approximate Numbers :            There are two types of numbers exact and approximate . Exact numbers are 2,4,9,13,7,2,6.45....etc . But there are numbers such as 4/3 =1.3333....  ,          √(2)= 1.41413 and finite number of digits . These may be approximated by numbers 1.3333, 1.4141 , 3.1416 respectively . Such numbers which represent the given numbers to certain degree of accuracy are called approximate numbers .  2. Significant Figures :                     The digits used to express a number is called significant digits (figures) . Thus each of the numbers 7845, 3.589, 0.4758 contains four significant figures while the numbers 0.00389 , 0.000587 and 0.0000296 contain only three significant figures since zeros only help to fix the position of the decimal point . Similarly the numbers 45000 and 7300.00 have two significant figures only . 3. Rounding Off :                   There are numbers with
Read more

PFAFFIAN Differential Equations And It's Solutions

-
Image
PFAFFIAN Differential Equations And It's Solutions :-      The relation              n             Σ  Fᵢ(x₁,x₂,......xₙ) dxᵢ = 0            i=1   where Fᵢ's (i=1,2,....n) are functions of all or some of the n independent variables x₁,x₂,.....xₙ is called a Pfaffian differential equation . Here we shall study Pfaffian differential equation in three independent variables which is of the form              P dx + Q dy + R dz = 0    ...........(1)    where P, Q , R are functions of x , y and z . Setting -->                        -->               X = (P,Q,R)  and dr = (dx,dy,dz), equation (1) can be written in the vector notation as                   -->  -->                  X . dr    = 0                  ..............(2)  If there exists a function μ(x,y,z) such that μ(P dx + Q dy + R dz) is an exact differential dΦ , the equation is said to be integrable and to possess an integrating factor μ(x,y,z) . In such case(1) becomes
Read more

Derivation Of Composite Trapezoidal Rule

-
Image
Derivation Of Composite Trapezoidal Rule :                                                   b       Let us find the value of ∫ f(x) dx                                                  a numerically by composite trapezoidal rule. Solution :              Let [a,b] be divided into n equal sub intervals of length h each . Let                            a= x₀,x₁,x₂,.......xₙ=b  be the points of sub division .  i.e x₁ = x₀+ih  for i= 0,1,2,......n By Newton's Forward Difference Interpolation  f(x) = f(x₀)+uΔf(x₀) +u(u-1)/2!  Δ²f(x₀)+........        where x = x₀+uh b                 ∫ f(x) dx = a       xₙ          ∫ [f(x₀)+uΔf(x₀)+u(u-1)/2!  Δ²f(x₀)+...]dx           x₀              n        = h∫ [f(x₀)+uΔf(x₀)+u(u-1)/2  Δ²f(x₀)+...]dx             0     =h[uf(x₀) +u²/2 Δf(x₀) +                                                                                                       n                                  (u³/6 -u²/4)Δ²f
Read more