Mth603 Formulae
Mth603 Formulae
Mth603 Formulae
x2 =
x0 + x1 2 + 4
=
=3
2
2
1. Muller Method
x=
x0 2c
+b b 2 4ac
in this formula x0,x1, x2 will be and U will just put the values according to the given formulas.
a=
h 2 f 1 ( h1 + h 2) f 0 + h1 f 2
h1h 2(h1 + h 2)
f 1 = f 0 ah12
b=
h1
c = f ( x 0)
h1 = x1 x 0
h 2 = x0 x 2
2. Regula Falsi Method (Method of False position)
x3 = x 2
xn xn 1
x 2 x1
f ( xn )
f ( x 2) xn +1 = xn
f ( xn ) f ( xn 1 )
f ( x 2) f ( x1)
f (x )
f ( x0)
xn +1 = xn ' n
'
f ( xn )
f ( x0)
4. Secant Method
x2 =
x f ( xn ) xn f ( xn 1 )
x0 f ( x1) x1 f ( x0)
xn +1 = n 1
f ( xn ) f ( xn 1 )
f ( x1) f ( x0)
5. Newtons Formula
1
x1 = ( x0 + n )
x0
2
in this formula x0= 2 perfect square near to 12 such like 9 and 16
6. Graffee root squaring method
f ( x) f ( x) = a32 x 6 (a2 2 2a1a3 ) x 4 + (a12 2a0 a2 ) x 2 x02
f ( x) = a32t 3 (a2 2 2a1a3 )t 2 + (a12 2a0 a2 )t a02 a is constant t is x
y x = y0 + py0 +
p ( p 1) 2
p ( p 1)( p 2) 3
y0 +
y0
2!
3!
Where p = (
x x0
)
h
y x = yn + pyn +
p( p + 1) 2
p( p + 1)( p + 2) 3
yn +
yn
2!
3!
p=(
x xn
)
h
y = f ( x) =
( x x1 )( x x2 )L ( x xn )
( x x0 )( x x2 )L ( x xn )
y0 +
y1 + L
( x0 x1 )( x0 x2 )L ( x0 xn )
( x1 x0 )( x1 x2 ) L ( x1 xn )
y[ x0 , x1 ] =
y1 y0
x1 x0
y[ x0 , x1 , x2 ] =
y1[ x1 , x2 ] y0 [ x0 , x1 ]
x2 x0
y[ x1 , x2 , x3 ] =
y2 [ x2 , x3 ] y1[ x1 , x2 ]
x3 x1
y[ x2 , x3 , x4 ] =
y3[ x3 , x4 ] y2 [ x2 , x3 ]
x4 x2
y[ x0 , x1 , x2 , x3 ] =
y1[ x1 , x2 , x3 ] y0 [ x0 , x1 , x2 ]
x3 x0
y[ x1 , x2 , x3 , x4 ] =
y2 [ x2 , x3 , x4 ] y1[ x1 , x2 , x3 ]
x4 x1
y[ x0 , x1 , x2 , x3 , x4 ] =
y1[ x1 , x2 , x3 , x4 ] y0 [ x0 , x1 , x2 , x3 ]
x4 x0
y( x) =
y ( x + h) y ( x h )
2h
y ( x ) = F (h) =
y ( x + h) y ( x h )
+ ET
2h
h
h
4m Fm 1 m Fm 1 m 1
h
2
2 where m = 1, 2,3,K
Fm m =
m
4 1
2
12.The Trapezoidal
xn
x0
ba
h
f ( x) dx = [ y0 + yn + 2( y1 + y2 + L + yn 1 )] + En h =
2
n
h
ba
fx )dx = [ y0 + yn + 4( y1 + y3 + y5 + ...) + 2( y2 + y4 + ....)] h =
3
2n
h
= [ y0 + 4( y1 + y3 + L + y2 N 1 ) + 2( y2 + y4 + L + y2 N 2 ) + y2 N ]
3
ERROR
E=
x2 N x0 4 (iv )
h y ( )
180
3
f
(
x
)
dx
=
h[ y (a) + 3 y1 + 3 y2 + 2 y3 + 3 y4 + 3 y5 + 2 y6 + L + 2 yn 3 + 3 yn 2 + 3 yn1 + y (b)]
a
8
ba
h=
3n
b
ERROR
E=
xn x0 4 ( iv )
h y ( )
80
ym +1 = ym + hf (tm , ym )
Here h =
approx value at x
given steps
yn +1 = yn +W1k1 +W2 k 2
This well-known fourth-order R-K method is described in the following steps
1
yn +1 = yn + (k1 + 2k2 + 2k3 + k4 )
6
where
k1 = hf (tn , yn )
h
k
k2 = hf tn + , yn + 1
2
2
k
h
k3 = hf tn + , yn + 2
2
2
k4 = hf (tn + h, yn + k3 )
h
C : yn(1)+1 = yn + f (tn , yn ) + f (t n+1 , yn(0)+1 )
2
P : yn(0)+1 = yn + hf (tn , yn )
Milnes Method
4h
(2 yn 2 yn 1 + 2 yn )
h
C : yn +1 = yn 1 + ( yn 1 + 4 yn + yn +1 )
3
P : yn +1 = yn 3 +
Adam-Moulton Method
h
( 55 yn 59 yn 1 + 37 yn 2 9 yn 3 )
24
C : yn +1 = yn + ( 9 yn +1 + 19 yn 5 yn 1 + yn 2 )
24
P : yn +1 = yn +
ROMBERGS INTEGRATION
h
h
4m IT ( m 1) m IT ( m1) m 1
2
2
=
m 1
4
h
ITm m
2
(t t0 ) 2
(t t 0 ) 3
(t t0 ) 4 IV
y (t0 ) +
y (t0 ) +
y (t0 ) +L
2!
3!
4!
AREA
b
f (x) dx
a
VOLUME
d b
g(x, y) dx dy
c a
Key concepts:
1. Integration is a summing process. Thus virtually all numerical
approximations can be represented by
n
b
I = f(x)dx = xf(x i ) + E t
i =0
a
where:
x = weights
xi = sampling points
Et = truncation error
ck are the weighting coefficients given by
b
h =
DIFFERENCE OPREATORS:
Applications: ( Primes)
2 y0 3 y0 4 y0
1
Dy0 = y0 = y0
+
+ L
h
2
3
4
2
d y0
1
11
5
D 2 y0 =
= y0 = 2 2 y0 3 y0 + 4 y0 5 y0 + L
2
dx
h
12
6
Backward Primes
2 y n 3 yn 4 y n
d
1
yn = Dyn = yn = yn
+
+
+ L
dx
h
2
3
4
1
11
5
yn = D 2 yn = 2 ( 2 yn + 3 yn + 4 y n + 5 yn + L
h
12
6
Central Primes
d
1
1
3 5
y = y = Dy = y 3 y +
y L
dx
h
24
640
1
1
1
y = D 2 y = 2 2 y 4 y + 6 y L
h
12
90
First derivative
yi =
yi +1 yi
h
Second derivative
yi =
yi + 2 2 yi +1 + yi
h2
Crouts Reduction
[L][U]=[A]