Numerical Solution of Linear Second Order Ordinary Differential Equations with Mixed Boundary Conditions by Galerkin Method
|
|
- Phebe Caldwell
- 6 years ago
- Views:
Transcription
1 Mathematcs ad Computer Scece 7; (5: do:.648/.mcs.75. Numercal Soluto of Lear Secod Order Ordary Dfferetal Equatos wth Mxed Boudary Codtos by Galer Method Aalu Abrham Aulo,, Alemayehu Shferaw Kbret, Geaew Gofe Gofa, Ayaa Deressa Negassa Departmet of Mathematcs, Isttute of Techology, Dre Dawa Uversty, Dre Dawa, Ethopa Departmet of Mathematcs, Jmma Uversty, Jmma, Ethopa Emal address: (A. A. Aulo, (A. S. Kbret, (G. G. Gofa, (A. D. Negassa Correspodg author To cte ths artcle: Aalu Abrham Aulo, Alemayehu Shferaw Kbret, Geaew Gofe Gofa, Ayaa Deressa Negassa. Numercal Soluto of Lear Secod Order Ordary Dfferetal Equatos wth Mxed Boudary Codtos by Galer Method. Mathematcs ad Computer Scece. Vol., No. 5, 7, pp do:.648/.mcs.75. Receved: Aprl 8, 7; Accepted: Jue 6, 7; Publshed: September 8, 7 Abstract: I ths paper, the Galer method s appled to secod order ordary dfferetal equato wth mxed boudary after covertg the gve lear secod order ordary dfferetal equato to equvalet boudary value problem by cosderg a vald assumpto for the depedet varable ad also covertg mxed boudary codto to Neuma type by usg secat ad Ruge-Kutta methods. The resultg system of equato s solved by drect method. I order to chec to what extet the method approxmates the exact soluto, a test example wth ow exact soluto s solved ad compared wth the exact soluto graphcally as well as umercally. Keywords: Secod Order Ordary Dfferetal Equato, Mxed Boudary Codtos, Ruge-Kutta, Secat Method, Galer Method, Chebyshev Polyomals. Itroducto The goal of umercal aalyss s to fd the approxmate umercal soluto to some real physcal problems by usg dfferet umercal techques, especally whe aalytcal solutos are ot avalable or very dffcult to obta. Sce most of mathematcal models of physcal pheomea are expressed terms of ordary dfferetal equatos, ad these equatos due to ther ature ad further applcatos to use computers, t eeds to establsh approprate umercal methods correspodg to the type of the dfferetal equato ad codtos that gover the mathematcal model of the physcal pheomea. The codtos may be specfed as a tal Value (IVP or at the boudares of the system, Boudary Value (BVP []. May problems egeerg ad scece ca be formulated as two-pot BVPs, le mechacal vbrato aalyss, vbrato of sprg, electrc crcut aalyss ad may others. Ths shows that the umercal methods used to approxmate the solutos of two-pot boudary value problems play a vtal role all braches of sceces ad egeerg []. Amog dfferet umercal methods used to approxmate two-pot boudary value problems terms of dfferetal equatos are shootg method, fte dfferece methods, fte elemet methods (FEM, Varatoal methods (Weghted resdual methods, Rtz method ad others have bee used to solve the two-pot boudary value problems [3]. Both FEM ad Varatoal methods the ma attempts were to loo a approxmato soluto the form of a lear combato of sutable approxmato fucto ad udetermed coeffcets [4]. For a vector space of fuctos S = φ ( x be bass of V, a set of learly V, f { } = depedet fuctos, ay fucto f ( x V could be uquely wrtte as a lear combato of the bass as:
2 Mathematcs ad Computer Scece 7; (5: f ( x = c φ ( x ( = The weghted resdual methods use a fte umber of φ ( x as tral fucto. learly depedet fuctos { } = Suppose that the approxmato soluto of the dfferetal equato, D( u = L( u( x + f ( x =, o the boudary B( u = [ a, b] s the form: N φ φ ( = u( x U ( x = c ( x + ( x N Where UN ( x s the approxmate soluto, u( x s the exact soluto L s a dfferetal operator, f s a gve fucto, φ ( x' s are fte umber of bass fuctos ad uow coeffcets for =,,..., N. The resdual R( x, c s defed as: R( x, c = D( U ( x ( L( U ( x + f ( x. Now determe N N c by requrg R to vash a weghted-resdual sese: b w ( x R( x, c = ( =,,..., N (3 a Where w ( x are a set of learly depedet fuctos, called weght fuctos, whch geeral ca be dfferet from the approxmato fuctos φ ( x, ths method s ow as the weghted-resdual method.. Galer Method If φ ( x = w ( x equato (3, the the specal ame of the weghted-resdual method s ow as the Galer method. Thus Galer method s oe of the weghted resdual methods whch the approxmato fucto s the same as the weght fucto ad hece t s also used to fd the approxmate soluto of two-pot boudary value problems [4]. The Galer method was veted 95 by Russa mathematca Bors Grgoryevch Galer ad the org of the method s geerally assocated wth a paper publshed by Galer 95 o the elastc equlbrum of rods ad th plates. He publshed hs fte elemet method 95. The Galer method ca be used to approxmate the soluto to ordary dfferetal equatos, partal dfferetal equatos ad tegral equatos [5]. May authors have bee used the Galer method to fd approxmate soluto of ordary dfferetal equatos wth boudary codto. Amog ths, a sple soluto of two pot boudary value problems troduced [6], a method for solutos of olear secod order mult-pot boudary value problems produced [7], [8] lear ad o-lear c dfferetal equatos were solved umercally by Galer method usg a Berste polyomals bass, [9] a umercal method s establshed to solved secod order ordary dfferetal equato wth Neuma ad Cauchy boudary codtos usg Hermte polyomals, [] a parametrc cubc sple soluto of two pot boudary value problems were obtaed, a secod-order Neuma boudary value problem wth sgular olearty for exact three postve solutos were solved [], a Numercal soluto of a sgular boudary-value problem o-newtoa flud mechacs were establshed [], a Fourth Order Boudary Value Problems by Galer Method wth Cubc B-sples were solved by cosderg dfferet cases o the boudary codto [3] ad a specal successve approxmatos method for solvg boudary value problems cludg ordary dfferetal equatos were proposed.[4] I ths paper Galer method wll be appled to the lear secod order ordary dfferetal equato of the form d y dy α( x + β ( x + δ ( x y = g( x ; a x b y( a = µ wth boudary codto ad y '( b = µ d y dy α( x + β( x + δ ( x y = g( x ; a x b wth y '( a = β boudary codto y '( b = β 3. Chebyshev Polyomal The polyomals whose propertes ad applcatos are dscussed ths paper were 'dscovered' almost a cetury ago by the Russa mathematca Chebyshev. It s a fucto defed usg trgoometrc fuctos cosθ ad sθ for x [,]. Chebyshev polyomals of frst d wth degree for x [, ] defed as: T ( x = cos θ, such that cos θ = x, for x ad Thus ( = T x cos ( cos x, T ( x = cos( = ad T ( x = cos(cos x = x From the trgoometrc detty, ( l θ ( l cos + + cos θ = cos θ cosθ s ( θ s( θ + s ( θ s( θ ( ( ( ( ( cos + l θ + cos l θ = cos θ cosθ T x = xt x - T x. Thus usg the recursve + - relato above for =, there s a seres of Chebyshev polyomal T ( x = x
3 68 Aalu Abrham Aulo et al.: Numercal Soluto of Lear Secod Order Ordary Dfferetal Equatos wth Mxed Boudary Codtos by Galer Method T ( x = x - 3 T 3 ( x = 4x - 3x 4 T 4 ( x = 8x -8x + Here the coeffcet of T ( x = 6 x x + 5 x etc. x T ( x s. The fgure below shows the graph of the frst seve Chebyshev polyomals for x [,]. Y-axs.5 The Graph of the frst eght Chebyshev Polyomals T(x T(x T3(x T4(x T5(x T6(x T7(x T8(x X-axs Fgure. The graph of the frst eght Chebyshev polyomals for x [,]. 4. Ruge-Kutta Method for Secod Order ODE Ruge-Kutta method s a umercal method used to fd approxmate soluto for tal value problems. I order to use Ruge-Kutta method to fd a approxmate soluto of secod order ODE, t eeds to covert to a system of two frst order ODEs. For two evaluato of f the method s gve by Where y+ = y + hy + ( K + K y + = y + hy + ( K + 3 K h h K = f ( x, y h 4 K = f ( x + h, y + hy + K Secat Method Ths method approxmates the graph of the fucto y = f ( x the eghborhood of the root by a straght le (secat passg through the pots ( x, f ( x ( x, f ( x, where f f ( x ad = ad tae the pot of tersecto of ths le wth the x-axs as the ext terate. Hece Or x x x x f + = f f x + x f x f = f f, =,,..., =,,... Where x ad x are two cosecutve terates. I ths method there are two tal approxmatos x ad x. Ths method s also called Secat method. 6. Mathematcal Formulato of the Method Cosder a geeral lear secod order dfferetal equato wth two type boudary codtos d y dy Type I: α( x + β( x + δ ( x y = g( x ; a x b wth boudary codto y( a = µ (Mxed type (4 y '( b = µ Type II: d y dy α( x + β( x + δ ( x y = g( x ; a x b wth boudary codto where fuctos for α( x, β ( x, δ ( x, y '( a = β y '( b = β (Neuma type (5 g( x are gve cotuous a x b where β, β, µ ad µ are gve
4 Mathematcs ad Computer Scece 7; (5: costats ad y( x s uow fucto or exact soluto of the boudary value problem whch s to be determed. I a BVP wth mxed boudary codto, the soluto s requred to satsfy a Drchlet or a Neuma boudary codto a mutually exclusve way o dsot parts of the boudary. Now Cosder the BVP of type II (Neuma type. To use a approxmatg polyomal defed for x [,] the gve BVP defed o arbtrary terval [a, b] must be coverted to a equvalet BVP defed o [-, ]. So that the approxmatg polyomal should be defed o [-, ]. Sce Chebyshev polyomal s defed o [-, ], t s possble to use Chebyshev polyomal after covertg the BVP defed o arbtrary terval [a, b] to a equvalet BVP defed o [-, ]. 6.. Coverso of the Doma of the BVP The dfferetal equato (5 together wth the Neuma boudary codto ca be coverted to a equvalet problem o [-, ] by lettg b a b + a x = x +, for x ad a x b The equato (4 wth boudary codto s equvalet to the BVP gve by ( d y dy ɶ α x + ɶ( x ( x y g( x; - x β + ɶ δ = ɶ Subect to the boudary codto, where y '( = d y '( = d 4 b a b + a b a b + a ɶ α( x = α x sce for x x, + = + ( b a Ths mples that dy dy d y 4 d y ad = ( b a = ( b a Now equatg the D. E (4 wth (6, d y dy α( x + β( x + δ ( x y g( x = d y dy ɶ α( x + ɶ β( x + ɶ δ ( x y gɶ ( x d y α x d y ɶ α( x = (, ɶ dy dy β ( x = β( x, ɶ δ ( x = δ ( x (6 (7 ad gɶ ( x = g( x (8 Therefore, the DE (5 wth Neuma boudary codto s a equvalet BVP wth the BVP (6. Up o substtuto of (7 to (4, equato (4 yelds d y d y α x = ɶ α x 4 ( (, ( b a dy dy β ( x = ɶ β( x, ( b a ɶ δ ( x y = δ ( x y ad gɶ ( x = g( x 4 b a b + a ɶ α( x = α( x +, ( b a ɶ b a b + a β ( x = β ( x +, ( b a ɶ b a b + a δ ( x= δ ( x + ad b a b + a gɶ ( x=g( x Applyg Galer Method To apply the techque of Galer method to fd a approxmate soluto of (4, say y( x, wrtte as a lear combato of base fuctos ad uow costats. That s; y( x = ct ( x (9 = where T ( x are pecewse polyomal, amely Chebyshev polyomals of degree ad c 's are uow parameters, to be determed. Now applyg Galer method wth the bass fucto T ( x gves, d y dy [ ɶ α ( x + ɶ β ( x + ɶ δ ( x y] T ( x = gɶ ( x T ( x ( Itegratg the frst term by parts o the left had sde of (, that s d y d y T ( x ɶ α( x, u = T ( ( ad x ɶ α x dv = d dy du = T ( x ɶ α( x ad v = d y ɶ α( x = uv vdu dy dy d = T( x ( x T( x ( x ɶ α α ɶ ( Upo substtuto of ( to (, yelds
5 7 Aalu Abrham Aulo et al.: Numercal Soluto of Lear Secod Order Ordary Dfferetal Equatos wth Mxed Boudary Codtos by Galer Method - - ( ɶ α dy d dy [- ( x T ( x + ɶ β ( x T ( x + ɶ δ ( x y( x T ( x] = gɶ ( x T ( x + ɶ α(- y '(- T (- - ɶ α( y '( T ( ( But from equato (9 the approxmate soluto s gve by y( x = ct ( x = Substtutg ths to equato ( yelds, - = ( ɶ α x T x d [- ct ( x ( ( + ɶ β( x ct ( x T ( x + ɶ δ ( x ct ( x T ( x] = = - = gɶ ( x T ( x + ɶ α(- uɶ (- T (- - ɶ α( uɶ ( T ( ( ɶ α d c [- T ( x ( x T ( x + ɶ β( x T ( x T ( x = - + ɶ δ ( x T ( x T ( x] = gɶ ( x T ( x + ɶ α(- y '(- T (- - - ɶ α( y '( T ( (3 I the left had sde of equato (3 above t eeds to ow the values of y (- ad y ( whch approxmately equal to y (- ad y ( respectvely, where y s the exact soluto of the BVP equato ( The Resultg System of Equato Sce the values of y ( ad y ( are ow from the boudary codto, substtutg these values to (3, ad equato (3 gves a system of equatos to solve the parameters c 's thus equato (3 matrx form becomes: ck = F (4 = ( d = ɶ ( [- T ( x ( x T ( x α ( = ɶ β (3 = ɶ ( δ ( x T ( x T ( x f - ( x T ( x T ( x = g ɶ ( x T ( x ( =ɶ α ɶ α f (- y (- T (-- ( y ( T ( Now, the uow parameters c 's are determed by solvg the system of equato (4 by drect method ad substtutg these values to (9 yelds the approxmate soluto y( x of the DE (4 satsfyg the gve boudary codtos (6. Cosder a BVP of type I. I ths case t s mpossble to use the above method drectly; sce y ( a s ot gve ad hece stead t eeds to covert the BVP to type II. The coverso s made by usg dfferet umercal methods. Cosder to solve the followg boudary-value problem: α d y dy ( x + ( x ( x y g( x ; a x b β + δ = (5 Wth boudary codto y( a = µ y '( b = µ The dea of shootg method for (5 s to solve for y ( a hopg that y ( b = µ. I order to fd y ( a such that y ( b = µ, guess y ( a = z ad solve for y ( b usg Ruge-Kutta method for secod order ODE, after havg a value usg the guess, deote ths approxmate soluto y z ad hope y z ( b = µ. If ot, use aother guess for y ( a, ad try to solve usg the Ruge-Kutta method. Ths process s repeated ad ca be doe systematcally utl ths choce satsfy y ( b. To do ths, follow the steps below. Step:- select z so that y z ( b = µ, letψ ( z = y z( b µ. The guess for z Step :- Now the obectve s smply to solve for ψ ( z =, hece secat method ca be used. Step 3:- How to compute z Suppose that the solutos yz ( b ad y z ( b obtaed from guesses z ad z respectvely. Step 4:- Now usg secat method to fd z gve by; z z z z y, =,,... + = z yz y z Followg ths sequece of terato there exsts z such that Where ( ( (3 ( ( ad K = + + F = f + f such that y ( b = y ( b z Thus, the Neuma boudary codto for the DE (5
6 Mathematcs ad Computer Scece 7; (5: s gve by y ( a = z y ( b = µ (6 Now to solve the DE wth boudary codto (6 t s coveet to use equato (3. Example: - Cosder the lear boudary value problem d y x + y = x e ; x, subect to the boudary codto y ( =.678, y ( = Whose exact soluto s:- y = / 5 s( x(349 s( e 5 / cos( / e ( ( x 3839 / 5 cos( x + / ( x + e Soluto: -The above problem s a mxed boudary codto or (type II; to apply the above method t eeds to covert the gve boudary codto to Neuma boudary codto. Now assume a guess depedg o the value of y ( =, let y ( = be the frst guess ad hopg that y ( =. The ext step s usg Ruge-Kutta method for secod order dfferetal equato, where y ( x = f ( x, y, y. But for ths problem, y ( x = f ( x, y sce f s depedet of y x =, ad x = ad tae step sze h=.5, ed ( y = y( x = for x = y =.678 y ( x = for x = y = y+ = y + hy + ( K + K y + = y + hy + ( K + 3 K h h K = f ( x, y h 4 K = f ( x + h, y + hy + K for =,, Ths gves the result table for the frst terato, where the th step x = x, y = y( x ad y = y ( x Referrg to table, tae y z =.588. But, y ( y (, thus t eeds to guess aother value for z y (. Let y ( = hopg that y ( =. Usg Ruge- Kutta method for x = to x = ad tag h =.5, y = y( x = for x = ad y = y ( x = for x =, ths yelds the followg result, where the y = y( x ad y = y ( x th step x = x, st terato d terato Now Table. Shows the result of y ad y o the frst terato. x y y Table. Shows the result of y ad y o secod terato. x y y y ( =.548 where z = ad sce z ad z gve the fd z z z z ( z = z yz = (.548 =.99 y y z z Usg z =.99 applyg Ruge-Kutta method where h=.5, y = y( x = for x = ad y = y ( x =.99 for x =
7 7 Aalu Abrham Aulo et al.: Numercal Soluto of Lear Secod Order Ordary Dfferetal Equatos wth Mxed Boudary Codtos by Galer Method gves the followg result for the 3 rd terato. 3 rd terato Table 3. Shows the result of y ad y o the thrd terato. x y y Now calculate the ext guess z z.99 z3 = z yz =.99 (.99 =.3 y y z z Ruge-Kutta method for y =.3 yelds the followg result for the 4 th terato. Table 4. Shows the result of y ad y o the fourth terato. 4 th terato x y y As table 4 shows the guess for y (.3. Thus, the Neuma boudary value problem gve by d y x + y = x e ; x y ( =.3, ad y ( = The ext step s covertg the BVP to equvalet BVP defed for x by lettg ( b a ( b + a x = x + = 5x + 5, sce a = ad b =. The equvalet BVP for the above problem o x becomes, d y ( 5x+ 5 + y = ( 5x + 5 e, for x 5 wth boudary codto (7 y'( =.3 y'( = (7 Now, suppose that y be the approxmate soluto of (7, gve by a lear combato of costats c ' s ad a approxmatg polyomal, called Chebyshev polyomal, thus y = ct ( x (8 = Upo substtuto of y, the approxmate soluto, to the dfferetal equato (7 gves a equato called resdue gve by: d y ( R( x, c y 5x 5 e 5 5x+ 5 = + ( +, for x (9 Applyg Galer method; R( c, x T ( x = d y ( 5x+ 5 [ + y ( 5x + 5 e ] T ( 5 x + = + 5 ( 5x+ 5 [ y '' T ( x yt ( x] ( 5x 5 e T ( x ( Usg tegrato by parts to smplfy the frst term the rght had sde, let u = T ( x du = T ( x ad dv = y ( x v = y ( x, hece t gves y T ( x y ( x T ( x y T ( x 5 5 = Therefore equato ( becomes
8 Mathematcs ad Computer Scece 7; (5: ( 5x+ 5 yt ( x y ' T '( x = y '( x T ( x + ( 5x + 5 e T ( x ( Substtutg the approxmate soluto y = ct ( x to ( yelds, = ( 5x+ 5 c ( ( '( '( '( ( T x T x T x T x = y x T x + ( 5x + 5 e T ( x 5 5 = = ( 5x x + 5 e T ( x c T ( x T ( x T ( x T ( x y ( T ( y = ( T ( ( ( Now, usg the gve boudary codto to (, equato ( becomes ( 5x+ 5 c T ( x T ( x T ( x T ( x = ( 5x + 5 e T ( x + (.3 T ( (3 5 5 = For equato (3 there s a system of equato gve by: ck = F (4 = where ( ( ( ( = + ad = + K F f f ( = such that T ( x T ( x ( = 5 T ( x T ( x I order to fd the value of fucto. For = 6 : x x 3 4x 3x T ' = 4 8x 8x x x + 5x 6 4 3x 48x + 8x ( f = (.3 T ( 5 ( ( 5x+ 5 f = 5x + 5 e T ( x ( c ' s tae tral fuctos defed for x [,], usg Chebyshev polyomals as tral, where T = [ T T T3 T4 T5 T6 ] ad T ' s the traspose of T.
9 74 Aalu Abrham Aulo et al.: Numercal Soluto of Lear Secod Order Ordary Dfferetal Equatos wth Mxed Boudary Codtos by Galer Method /3 / 5 / 4 /5 38 /5 6 / 35 / 5 34 / 35 / 63 = T ( x T ( x = 38 /5 6 / / 99 / / / 99 6 / / 99 4 /43 ( /5 /5 /5 3/75 8/375 88/875 /5 38/5 4/75 = T ( x T ( x = 5 8/ / /65 /5 4/75 6/35 88/ /65 564/965 ( To fd the value of the coeffcet matrx equato (4 use ( ( = + K So, by expressg the coeffcet matrx K, ad the uow coeffcet c ' s a system of equato matrx form: ( ( = + K 44/75 / 5 9 / 55 c 38 / / / 7875 c / 5 6 / /575 c3 = [ f ] 846 / / / 8665 c4 9 / / / 3465 c5 34 / / / 55 c6 (5 Where f s 6 a colum vector gve by:- (- - 4 / / 5e (- - / /5e (- 8 /5-478 /5e f = (- -38 / / 65e 8 / / 35e 88 / /565e (- (- Now a 6 6 coeffcet matrx whch s symmetrc, 6 uow colum vector that represet c ' s ad 6 colum vector that represets f. So there are sx equatos wth sx uows. Usg ( values of the sx uows are:, c = K f to solve (5, the c =.3655 c = c 3 = c 4 =.3697 c 5 = c 6 =.9 Now t s possble to express the approxmate soluto as a lear combato of costats c ' s ad a approxmatg polyomal. So substtutg c ' s ad polyomals for =6 the approxmate soluto s: Chebyshev 3 4 y =.3655x.596(x (4x 3 x (8x 8x (6 x x + 5 x.9(3x 48x + 8x The graph of the exact ad approxmate soluto, to loo the covergece of the approxmate soluto to the graph of the exact soluto, loos le fgure below.
10 Mathematcs ad Computer Scece 7; (5: The Graph of Exact ad Approxmate Solutos for =4 ad =6 The exact soluto Approx. soluto for =4 Approx. soluto for =6.5 Y-axs X-axs -<x< Fgure. The graph of the exact ad approxmate soluto for =4 ad =6. The above graph shows that the graph of the approxmate soluto approaches the graph of the exact soluto for the dfferetal equato wth boudary codto problem. Cosderg =8, the the approxmate soluto ad the Chebyshev polyomals are gve by: 8 y = ct ( x (6 x x 3 4x 3x 4 8x 8x + T ' =, where T = [ T 5 3 6x x 5x T T3 T4 T5 T6 T7 T8 ] ad T ' s the traspose of T x 48x + 8x x x + 56x 7x x + 6x 56x + 8x Usg MATLAB code; = ( K = /3 -/5 -/ -/45 4/5-38/5-6/35-34/3465 -/5 34/35 -/63-38/495 ] -38/5 6/63-34/99-58/45 -/ -/63 98/99-46/49-6/35-34/99 4/43 -/65 -/45-38/495-46/49 94/95-34/ /45 -/65 54/55
11 76 Aalu Abrham Aulo et al.: Numercal Soluto of Lear Secod Order Ordary Dfferetal Equatos wth Mxed Boudary Codtos by Galer Method ( K = /5 /5 /5 /5 3/75 8/375 88/875 5/575 /5 38/5 4/75 86/375 8/ / /65 83/8665 /5 4/75 6/35 686/475 88/ /65 564/ / /5 86/ / /6875 5/575 83/ / /65 f = (- - 4 / / 5e (- - / /5e (- 8 /5-478 /5e (- -38 / / 65e (- 8 / / 35e 88 / /565e 3544 / / 785e 544 / / 785e (- (- (- Thus the uow parameters c ' s are c = c = c 3 = c 4 =.8693 c 5 =.4546 c 6 =.79 c 7 =.579 c 8 =.9963 Now substtutg the uow parameters ad eght Chebyshev polyomals to (6, the approxmate soluto s: y = x x + x x + x x (.95956( ( (6x x 5 x.79(3x 48x 8x (64x x 56x 7 x.9963( 3x 6x 56x 8 x The graph of the exact ad approxmate soluto, to loo the covergece of the approxmate soluto to the graph of the exact soluto, loos le fgure 3 below.
12 Mathematcs ad Computer Scece 7; (5: The graph of the exact ad approxmate souluto for =6 ad =8 Exact soluto Approxmate soluto for =6 Approxmate soluto for =8.5 Y-axs X-axs -<=X<= Fgure 3. The graph of the exact soluto ad approxmate soluto for =6 ad =8. From the graph above the approxmate soluto approaches the graph of the exact soluto whe the umber of the tral fuctos creases from 6 to 8. Now compare the absolute error whch s gve by error = y y exact approx Table 5. Shows the computed the absolute error for problem whe =6 ad =8. x Exact soluto Approxmate soluto for =6 Approxmate soluto for =8 Absolute error for =6 Absolute error for = As observed from table 5 the approxmate soluto s approachg to the exact value as the value of creases. Tae =, the graph of the correspodg approxmate soluto together wth the graph of the approxmate soluto for =6 ad =8 the same plae wth the exact soluto s show below.
13 78 Aalu Abrham Aulo et al.: Numercal Soluto of Lear Secod Order Ordary Dfferetal Equatos wth Mxed Boudary Codtos by Galer Method.5.5 The Graph of Exact ad Approxmate solutos for =6,=8 ad = The exact soluto Approx. soluto for =6 Approx. soluto for =8 Approx. soluto for = Y-axs X-axs -<x< Fgure 4. The graph of the exact soluto ad approxmate soluto for =6, =8 ad =. As the umber of Chebyshev polyomal creases the correspodg approxmate soluto of the dfferetal equato wth mxed boudary codto problem approaches the graph of the exact soluto. 7. Cocluso Ths study troduces that, by applyg Galer method to lear secod order ordary dfferetal equatos wth mxed ad Neuma boudary codtos, t s possble to fd ther approxmate solutos. The umercal results that are obtaed usg ths method coverges to the exact soluto as the umber of Chebyshev polyomal creases, that wll be used as a tral fucto; ad also usg small step sze h, whle covertg the gve lear secod order ordary dfferetal equato from mxed type to Neuma boudary codto, creases the accuracy of the approxmate soluto. So that usg ths method better results wll be obtaed as the umber of Chebyshev polyomal creases ad usg small step sze whle usg Rug-Kutta method. 8. Future Scope Ths study has led to a attetveess of several topcs that requre further vestgato, for stace lear ordary dfferetal equatos wth dfferet order, o-lear ordary dfferetal equatos. I order to fll the gap terms of accuracy t s mportat to aalyze the error of the method. Thus aalyzg the error ad creasg the accuracy of ths method s left for future vestgato. Refereces [] Yatteder Rsh Dubey: A approxmate soluto to buclg of plates by the Galer method, (August 5. [] Ta-Ra Hsu: Mechacal Egeerg 3 Appled Egeerg aalyss, Sa Jose State Uversty, (Sept 9. [3] E. Sul: Numercal Soluto of Ordary Dfferetal Equatos, (Aprl 3. [4] J. N. Reddy: A Itroducto to the fte elemet method, 3rd edto, McGraw-Hll, (Ja [5] Marcos Cesar Rugger: Theory of Galer method ad explaato of MATLAB code, (6. [6] Jall Rashda ad Reza Jalla: Sple soluto of two pot boudary value problems, Appl. Comput. Math 9 ( [7] S. Das, Sul Kumar ad O. P. Sgh: Solutos of olear secod order multpot boudary value problems by Homotopy perturbato method, Appl. Appl. Math. 5 ( [8] M. Idress Bhatt ad P. Brace: Solutos of dfferetal equatos a Berste polyomals bass, J. Comput. Appl. Math. 5 ( [9] M. M. Rahma. et.al: Numercal Solutos of Secod Order Boudary Value Problems by Galer Method wth Hermte Polyomals, (. [] Arshad Kha: Parametrc cubc sple soluto of two pot boudary value problems, Appl. Math. Comput. 54 ( [] Yuqag Feg ad Guagu L: Exact three postve solutos to a secod-order Neuma boudary value problem wth sgular olearty, Araba J. Sc. Eg. 35 ( [] P. M. Lma ad M. Carpeter: Numercal soluto of a sgular boudary-value problem o-newtoa flud mechacs, Computer Phys. Commuca. 6( 4-. [3] K. N. S. Kas Vswaadham ad Sreevasulu Ballem: Fourth Order Boudary Value Problems by Galer Method wth Cubc B-sples, (May 3. [4] Jahashah et al.: A specal successve approxmatos method for solvg boudary value problems cludg ordary dfferetal equatos, (August 3. [5] L. Fox ad I. B. Parer: Chebyshev Polyomals Numercal Aalyss, Oxford Uversty Press, ( May, 967.
Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem
Joural of Amerca Scece ;6( Cubc Nopolyomal Sple Approach to the Soluto of a Secod Order Two-Pot Boudary Value Problem W.K. Zahra, F.A. Abd El-Salam, A.A. El-Sabbagh ad Z.A. ZAk * Departmet of Egeerg athematcs
More informationBERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS. Aysegul Akyuz Dascioglu and Nese Isler
Mathematcal ad Computatoal Applcatos, Vol. 8, No. 3, pp. 293-300, 203 BERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS Aysegul Ayuz Dascoglu ad Nese Isler Departmet of Mathematcs,
More informationChapter 9 Jordan Block Matrices
Chapter 9 Jorda Block atrces I ths chapter we wll solve the followg problem. Gve a lear operator T fd a bass R of F such that the matrx R (T) s as smple as possble. f course smple s a matter of taste.
More informationUNIT 2 SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS
Numercal Computg -I UNIT SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS Structure Page Nos..0 Itroducto 6. Objectves 7. Ital Approxmato to a Root 7. Bsecto Method 8.. Error Aalyss 9.4 Regula Fals Method
More informationDerivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations
Dervato of -Pot Block Method Formula for Solvg Frst Order Stff Ordary Dfferetal Equatos Kharul Hamd Kharul Auar, Kharl Iskadar Othma, Zara Bb Ibrahm Abstract Dervato of pot block method formula wth costat
More informationChapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II
CEE49b Chapter - Free Vbrato of Mult-Degree-of-Freedom Systems - II We ca obta a approxmate soluto to the fudametal atural frequecy through a approxmate formula developed usg eergy prcples by Lord Raylegh
More informationC-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory
ROAD MAP... AE301 Aerodyamcs I UNIT C: 2-D Arfols C-1: Aerodyamcs of Arfols 1 C-2: Aerodyamcs of Arfols 2 C-3: Pael Methods C-4: Th Arfol Theory AE301 Aerodyamcs I Ut C-3: Lst of Subects Problem Solutos?
More informationFunctions of Random Variables
Fuctos of Radom Varables Chapter Fve Fuctos of Radom Varables 5. Itroducto A geeral egeerg aalyss model s show Fg. 5.. The model output (respose) cotas the performaces of a system or product, such as weght,
More informationApplication of Legendre Bernstein basis transformations to degree elevation and degree reduction
Computer Aded Geometrc Desg 9 79 78 www.elsever.com/locate/cagd Applcato of Legedre Berste bass trasformatos to degree elevato ad degree reducto Byug-Gook Lee a Yubeom Park b Jaechl Yoo c a Dvso of Iteret
More informationAnalysis of Lagrange Interpolation Formula
P IJISET - Iteratoal Joural of Iovatve Scece, Egeerg & Techology, Vol. Issue, December 4. www.jset.com ISS 348 7968 Aalyss of Lagrage Iterpolato Formula Vjay Dahya PDepartmet of MathematcsMaharaja Surajmal
More informationLecture 07: Poles and Zeros
Lecture 07: Poles ad Zeros Defto of poles ad zeros The trasfer fucto provdes a bass for determg mportat system respose characterstcs wthout solvg the complete dfferetal equato. As defed, the trasfer fucto
More informationNumerical Analysis Formulae Booklet
Numercal Aalyss Formulae Booklet. Iteratve Scemes for Systems of Lear Algebrac Equatos:.... Taylor Seres... 3. Fte Dfferece Approxmatos... 3 4. Egevalues ad Egevectors of Matrces.... 3 5. Vector ad Matrx
More informationMOLECULAR VIBRATIONS
MOLECULAR VIBRATIONS Here we wsh to vestgate molecular vbratos ad draw a smlarty betwee the theory of molecular vbratos ad Hückel theory. 1. Smple Harmoc Oscllator Recall that the eergy of a oe-dmesoal
More informationA NEW NUMERICAL APPROACH FOR SOLVING HIGH-ORDER LINEAR AND NON-LINEAR DIFFERANTIAL EQUATIONS
Secer, A., et al.: A New Numerıcal Approach for Solvıg Hıgh-Order Lıear ad No-Lıear... HERMAL SCIENCE: Year 8, Vol., Suppl., pp. S67-S77 S67 A NEW NUMERICAL APPROACH FOR SOLVING HIGH-ORDER LINEAR AND NON-LINEAR
More informationA Conventional Approach for the Solution of the Fifth Order Boundary Value Problems Using Sixth Degree Spline Functions
Appled Matheatcs, 1, 4, 8-88 http://d.do.org/1.4/a.1.448 Publshed Ole Aprl 1 (http://www.scrp.org/joural/a) A Covetoal Approach for the Soluto of the Ffth Order Boudary Value Probles Usg Sth Degree Sple
More informationFourth Order Four-Stage Diagonally Implicit Runge-Kutta Method for Linear Ordinary Differential Equations ABSTRACT INTRODUCTION
Malasa Joural of Mathematcal Sceces (): 95-05 (00) Fourth Order Four-Stage Dagoall Implct Ruge-Kutta Method for Lear Ordar Dfferetal Equatos Nur Izzat Che Jawas, Fudzah Ismal, Mohamed Sulema, 3 Azm Jaafar
More informationAssignment 5/MATH 247/Winter Due: Friday, February 19 in class (!) (answers will be posted right after class)
Assgmet 5/MATH 7/Wter 00 Due: Frday, February 9 class (!) (aswers wll be posted rght after class) As usual, there are peces of text, before the questos [], [], themselves. Recall: For the quadratc form
More informationPGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation
PGE 30: Formulato ad Soluto Geosystems Egeerg Dr. Balhoff Iterpolato Numercal Methods wth MATLAB, Recktewald, Chapter 0 ad Numercal Methods for Egeers, Chapra ad Caale, 5 th Ed., Part Fve, Chapter 8 ad
More information2006 Jamie Trahan, Autar Kaw, Kevin Martin University of South Florida United States of America
SOLUTION OF SYSTEMS OF SIMULTANEOUS LINEAR EQUATIONS Gauss-Sedel Method 006 Jame Traha, Autar Kaw, Kev Mart Uversty of South Florda Uted States of Amerca kaw@eg.usf.edu Itroducto Ths worksheet demostrates
More informationd dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin
Learzato of the Swg Equato We wll cover sectos.5.-.6 ad begg of Secto 3.3 these otes. 1. Sgle mache-fte bus case Cosder a sgle mache coected to a fte bus, as show Fg. 1 below. E y1 V=1./_ Fg. 1 The admttace
More informationThe Mathematical Appendix
The Mathematcal Appedx Defto A: If ( Λ, Ω, where ( λ λ λ whch the probablty dstrbutos,,..., Defto A. uppose that ( Λ,,..., s a expermet type, the σ-algebra o λ λ λ are defed s deoted by ( (,,...,, σ Ω.
More informationANALYSIS ON THE NATURE OF THE BASIC EQUATIONS IN SYNERGETIC INTER-REPRESENTATION NETWORK
Far East Joural of Appled Mathematcs Volume, Number, 2008, Pages Ths paper s avalable ole at http://www.pphm.com 2008 Pushpa Publshg House ANALYSIS ON THE NATURE OF THE ASI EQUATIONS IN SYNERGETI INTER-REPRESENTATION
More informationOn the convergence of derivatives of Bernstein approximation
O the covergece of dervatves of Berste approxmato Mchael S. Floater Abstract: By dfferetatg a remader formula of Stacu, we derve both a error boud ad a asymptotc formula for the dervatves of Berste approxmato.
More informationL5 Polynomial / Spline Curves
L5 Polyomal / Sple Curves Cotets Coc sectos Polyomal Curves Hermte Curves Bezer Curves B-Sples No-Uform Ratoal B-Sples (NURBS) Mapulato ad Represetato of Curves Types of Curve Equatos Implct: Descrbe a
More informationLecture 12 APPROXIMATION OF FIRST ORDER DERIVATIVES
FDM: Appromato of Frst Order Dervatves Lecture APPROXIMATION OF FIRST ORDER DERIVATIVES. INTRODUCTION Covectve term coservato equatos volve frst order dervatves. The smplest possble approach for dscretzato
More informationSolving Constrained Flow-Shop Scheduling. Problems with Three Machines
It J Cotemp Math Sceces, Vol 5, 2010, o 19, 921-929 Solvg Costraed Flow-Shop Schedulg Problems wth Three Maches P Pada ad P Rajedra Departmet of Mathematcs, School of Advaced Sceces, VIT Uversty, Vellore-632
More information16 Homework lecture 16
Quees College, CUNY, Departmet of Computer Scece Numercal Methods CSCI 361 / 761 Fall 2018 Istructor: Dr. Sateesh Mae c Sateesh R. Mae 2018 16 Homework lecture 16 Please emal your soluto, as a fle attachmet,
More informationA Collocation Method for Solving Abel s Integral Equations of First and Second Kinds
A Collocato Method for Solvg Abel s Itegral Equatos of Frst ad Secod Kds Abbas Saadatmad a ad Mehd Dehgha b a Departmet of Mathematcs, Uversty of Kasha, Kasha, Ira b Departmet of Appled Mathematcs, Faculty
More informationArithmetic Mean and Geometric Mean
Acta Mathematca Ntresa Vol, No, p 43 48 ISSN 453-6083 Arthmetc Mea ad Geometrc Mea Mare Varga a * Peter Mchalča b a Departmet of Mathematcs, Faculty of Natural Sceces, Costate the Phlosopher Uversty Ntra,
More informationh-analogue of Fibonacci Numbers
h-aalogue of Fboacc Numbers arxv:090.0038v [math-ph 30 Sep 009 H.B. Beaoum Prce Mohammad Uversty, Al-Khobar 395, Saud Araba Abstract I ths paper, we troduce the h-aalogue of Fboacc umbers for o-commutatve
More information02/15/04 INTERESTING FINITE AND INFINITE PRODUCTS FROM SIMPLE ALGEBRAIC IDENTITIES
0/5/04 ITERESTIG FIITE AD IFIITE PRODUCTS FROM SIMPLE ALGEBRAIC IDETITIES Thomas J Osler Mathematcs Departmet Rowa Uversty Glassboro J 0808 Osler@rowaedu Itroducto The dfferece of two squares, y = + y
More informationResearch Article A New Iterative Method for Common Fixed Points of a Finite Family of Nonexpansive Mappings
Hdaw Publshg Corporato Iteratoal Joural of Mathematcs ad Mathematcal Sceces Volume 009, Artcle ID 391839, 9 pages do:10.1155/009/391839 Research Artcle A New Iteratve Method for Commo Fxed Pots of a Fte
More informationMultivariate Transformation of Variables and Maximum Likelihood Estimation
Marquette Uversty Multvarate Trasformato of Varables ad Maxmum Lkelhood Estmato Dael B. Rowe, Ph.D. Assocate Professor Departmet of Mathematcs, Statstcs, ad Computer Scece Copyrght 03 by Marquette Uversty
More informationSolution of General Dual Fuzzy Linear Systems. Using ABS Algorithm
Appled Mathematcal Sceces, Vol 6, 0, o 4, 63-7 Soluto of Geeral Dual Fuzzy Lear Systems Usg ABS Algorthm M A Farborz Aragh * ad M M ossezadeh Departmet of Mathematcs, Islamc Azad Uversty Cetral ehra Brach,
More informationEVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM
EVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM Jose Javer Garca Moreta Ph. D research studet at the UPV/EHU (Uversty of Basque coutry) Departmet of Theoretcal
More information1 Lyapunov Stability Theory
Lyapuov Stablty heory I ths secto we cosder proofs of stablty of equlbra of autoomous systems. hs s stadard theory for olear systems, ad oe of the most mportat tools the aalyss of olear systems. It may
More informationLINEAR RECURRENT SEQUENCES AND POWERS OF A SQUARE MATRIX
INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 6 2006, #A12 LINEAR RECURRENT SEQUENCES AND POWERS OF A SQUARE MATRIX Hacèe Belbachr 1 USTHB, Departmet of Mathematcs, POBox 32 El Ala, 16111,
More informationbest estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best
Error Aalyss Preamble Wheever a measuremet s made, the result followg from that measuremet s always subject to ucertaty The ucertaty ca be reduced by makg several measuremets of the same quatty or by mprovg
More informationNon-uniform Turán-type problems
Joural of Combatoral Theory, Seres A 111 2005 106 110 wwwelsevercomlocatecta No-uform Turá-type problems DhruvMubay 1, Y Zhao 2 Departmet of Mathematcs, Statstcs, ad Computer Scece, Uversty of Illos at
More information. The set of these sums. be a partition of [ ab, ]. Consider the sum f( x) f( x 1)
Chapter 7 Fuctos o Bouded Varato. Subject: Real Aalyss Level: M.Sc. Source: Syed Gul Shah (Charma, Departmet o Mathematcs, US Sargodha Collected & Composed by: Atq ur Rehma (atq@mathcty.org, http://www.mathcty.org
More informationOverview of the weighting constants and the points where we evaluate the function for The Gaussian quadrature Project two
Overvew of the weghtg costats ad the pots where we evaluate the fucto for The Gaussa quadrature Project two By Ashraf Marzouk ChE 505 Fall 005 Departmet of Mechacal Egeerg Uversty of Teessee Koxvlle, TN
More informationX ε ) = 0, or equivalently, lim
Revew for the prevous lecture Cocepts: order statstcs Theorems: Dstrbutos of order statstcs Examples: How to get the dstrbuto of order statstcs Chapter 5 Propertes of a Radom Sample Secto 55 Covergece
More informationDynamic Analysis of Axially Beam on Visco - Elastic Foundation with Elastic Supports under Moving Load
Dyamc Aalyss of Axally Beam o Vsco - Elastc Foudato wth Elastc Supports uder Movg oad Saeed Mohammadzadeh, Seyed Al Mosayeb * Abstract: For dyamc aalyses of ralway track structures, the algorthm of soluto
More informationComplete Convergence and Some Maximal Inequalities for Weighted Sums of Random Variables
Joural of Sceces, Islamc Republc of Ira 8(4): -6 (007) Uversty of Tehra, ISSN 06-04 http://sceces.ut.ac.r Complete Covergece ad Some Maxmal Iequaltes for Weghted Sums of Radom Varables M. Am,,* H.R. Nl
More informationChapter 5 Properties of a Random Sample
Lecture 6 o BST 63: Statstcal Theory I Ku Zhag, /0/008 Revew for the prevous lecture Cocepts: t-dstrbuto, F-dstrbuto Theorems: Dstrbutos of sample mea ad sample varace, relatoshp betwee sample mea ad sample
More informationBeam Warming Second-Order Upwind Method
Beam Warmg Secod-Order Upwd Method Petr Valeta Jauary 6, 015 Ths documet s a part of the assessmet work for the subject 1DRP Dfferetal Equatos o Computer lectured o FNSPE CTU Prague. Abstract Ths documet
More informationESS Line Fitting
ESS 5 014 17. Le Fttg A very commo problem data aalyss s lookg for relatoshpetwee dfferet parameters ad fttg les or surfaces to data. The smplest example s fttg a straght le ad we wll dscuss that here
More information( ) 2 2. Multi-Layer Refraction Problem Rafael Espericueta, Bakersfield College, November, 2006
Mult-Layer Refracto Problem Rafael Espercueta, Bakersfeld College, November, 006 Lght travels at dfferet speeds through dfferet meda, but refracts at layer boudares order to traverse the least-tme path.
More informationGeneralized One-Step Third Derivative Implicit Hybrid Block Method for the Direct Solution of Second Order Ordinary Differential Equation
Appled Mathematcal Sceces, Vol. 1, 16, o. 9, 417-4 HIKARI Ltd, www.m-hkar.com http://dx.do.org/1.1988/ams.16.51667 Geeralzed Oe-Step Thrd Dervatve Implct Hybrd Block Method for the Drect Soluto of Secod
More informationQuantization in Dynamic Smarandache Multi-Space
Quatzato Dyamc Smaradache Mult-Space Fu Yuhua Cha Offshore Ol Research Ceter, Beg, 7, Cha (E-mal: fuyh@cooc.com.c ) Abstract: Dscussg the applcatos of Dyamc Smaradache Mult-Space (DSMS) Theory. Supposg
More information4 Inner Product Spaces
11.MH1 LINEAR ALGEBRA Summary Notes 4 Ier Product Spaces Ier product s the abstracto to geeral vector spaces of the famlar dea of the scalar product of two vectors or 3. I what follows, keep these key
More informationx y exp λ'. x exp λ 2. x exp 1.
egecosmcd Egevalue-egevector of the secod dervatve operator d /d hs leads to Fourer seres (se, cose, Legedre, Bessel, Chebyshev, etc hs s a eample of a systematc way of geeratg a set of mutually orthogoal
More informationStrong Convergence of Weighted Averaged Approximants of Asymptotically Nonexpansive Mappings in Banach Spaces without Uniform Convexity
BULLETIN of the MALAYSIAN MATHEMATICAL SCIENCES SOCIETY Bull. Malays. Math. Sc. Soc. () 7 (004), 5 35 Strog Covergece of Weghted Averaged Appromats of Asymptotcally Noepasve Mappgs Baach Spaces wthout
More informationA unified matrix representation for degree reduction of Bézier curves
Computer Aded Geometrc Desg 21 2004 151 164 wwwelsevercom/locate/cagd A ufed matrx represetato for degree reducto of Bézer curves Hask Suwoo a,,1, Namyog Lee b a Departmet of Mathematcs, Kokuk Uversty,
More informationIntroduction to local (nonparametric) density estimation. methods
Itroducto to local (oparametrc) desty estmato methods A slecture by Yu Lu for ECE 66 Sprg 014 1. Itroducto Ths slecture troduces two local desty estmato methods whch are Parze desty estmato ad k-earest
More informationON THE LOGARITHMIC INTEGRAL
Hacettepe Joural of Mathematcs ad Statstcs Volume 39(3) (21), 393 41 ON THE LOGARITHMIC INTEGRAL Bra Fsher ad Bljaa Jolevska-Tueska Receved 29:9 :29 : Accepted 2 :3 :21 Abstract The logarthmc tegral l(x)
More informationLecture Note to Rice Chapter 8
ECON 430 HG revsed Nov 06 Lecture Note to Rce Chapter 8 Radom matrces Let Y, =,,, m, =,,, be radom varables (r.v. s). The matrx Y Y Y Y Y Y Y Y Y Y = m m m s called a radom matrx ( wth a ot m-dmesoal dstrbuto,
More informationLINEAR REGRESSION ANALYSIS
LINEAR REGRESSION ANALYSIS MODULE V Lecture - Correctg Model Iadequaces Through Trasformato ad Weghtg Dr. Shalabh Departmet of Mathematcs ad Statstcs Ida Isttute of Techology Kapur Aalytcal methods for
More information1 0, x? x x. 1 Root finding. 1.1 Introduction. Solve[x^2-1 0,x] {{x -1},{x 1}} Plot[x^2-1,{x,-2,2}] 3
Adrew Powuk - http://www.powuk.com- Math 49 (Numercal Aalyss) Root fdg. Itroducto f ( ),?,? Solve[^-,] {{-},{}} Plot[^-,{,-,}] Cubc equato https://e.wkpeda.org/wk/cubc_fucto Quartc equato https://e.wkpeda.org/wk/quartc_fucto
More informationThird handout: On the Gini Index
Thrd hadout: O the dex Corrado, a tala statstca, proposed (, 9, 96) to measure absolute equalt va the mea dfferece whch s defed as ( / ) where refers to the total umber of dvduals socet. Assume that. The
More informationMaximum Likelihood Estimation
Marquette Uverst Maxmum Lkelhood Estmato Dael B. Rowe, Ph.D. Professor Departmet of Mathematcs, Statstcs, ad Computer Scece Coprght 08 b Marquette Uverst Maxmum Lkelhood Estmato We have bee sag that ~
More informationJournal of Mathematical Analysis and Applications
J. Math. Aal. Appl. 365 200) 358 362 Cotets lsts avalable at SceceDrect Joural of Mathematcal Aalyss ad Applcatos www.elsever.com/locate/maa Asymptotc behavor of termedate pots the dfferetal mea value
More informationCHAPTER 4 RADICAL EXPRESSIONS
6 CHAPTER RADICAL EXPRESSIONS. The th Root of a Real Number A real umber a s called the th root of a real umber b f Thus, for example: s a square root of sce. s also a square root of sce ( ). s a cube
More informationEstimation of Stress- Strength Reliability model using finite mixture of exponential distributions
Iteratoal Joural of Computatoal Egeerg Research Vol, 0 Issue, Estmato of Stress- Stregth Relablty model usg fte mxture of expoetal dstrbutos K.Sadhya, T.S.Umamaheswar Departmet of Mathematcs, Lal Bhadur
More informationarxiv: v4 [math.nt] 14 Aug 2015
arxv:52.799v4 [math.nt] 4 Aug 25 O the propertes of terated bomal trasforms for the Padova ad Perr matrx sequeces Nazmye Ylmaz ad Necat Tasara Departmet of Mathematcs, Faculty of Scece, Selcu Uversty,
More informationA Remark on the Uniform Convergence of Some Sequences of Functions
Advaces Pure Mathematcs 05 5 57-533 Publshed Ole July 05 ScRes. http://www.scrp.org/joural/apm http://dx.do.org/0.436/apm.05.59048 A Remark o the Uform Covergece of Some Sequeces of Fuctos Guy Degla Isttut
More informationn -dimensional vectors follow naturally from the one
B. Vectors ad sets B. Vectors Ecoomsts study ecoomc pheomea by buldg hghly stylzed models. Uderstadg ad makg use of almost all such models requres a hgh comfort level wth some key mathematcal sklls. I
More informationLecture 3 Probability review (cont d)
STATS 00: Itroducto to Statstcal Iferece Autum 06 Lecture 3 Probablty revew (cot d) 3. Jot dstrbutos If radom varables X,..., X k are depedet, the ther dstrbuto may be specfed by specfyg the dvdual dstrbuto
More information{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:
Chapter 4 Exercses Samplg Theory Exercse (Smple radom samplg: Let there be two correlated radom varables X ad A sample of sze s draw from a populato by smple radom samplg wthout replacemet The observed
More informationMu Sequences/Series Solutions National Convention 2014
Mu Sequeces/Seres Solutos Natoal Coveto 04 C 6 E A 6C A 6 B B 7 A D 7 D C 7 A B 8 A B 8 A C 8 E 4 B 9 B 4 E 9 B 4 C 9 E C 0 A A 0 D B 0 C C Usg basc propertes of arthmetc sequeces, we fd a ad bm m We eed
More informationDKA method for single variable holomorphic functions
DKA method for sgle varable holomorphc fuctos TOSHIAKI ITOH Itegrated Arts ad Natural Sceces The Uversty of Toushma -, Mamhosama, Toushma, 770-8502 JAPAN Abstract: - Durad-Kerer-Aberth (DKA method for
More informationTESTS BASED ON MAXIMUM LIKELIHOOD
ESE 5 Toy E. Smth. The Basc Example. TESTS BASED ON MAXIMUM LIKELIHOOD To llustrate the propertes of maxmum lkelhood estmates ad tests, we cosder the smplest possble case of estmatg the mea of the ormal
More informationSimple Linear Regression
Statstcal Methods I (EST 75) Page 139 Smple Lear Regresso Smple regresso applcatos are used to ft a model descrbg a lear relatoshp betwee two varables. The aspects of least squares regresso ad correlato
More informationNewton s Power Flow algorithm
Power Egeerg - Egll Beedt Hresso ewto s Power Flow algorthm Power Egeerg - Egll Beedt Hresso The ewto s Method of Power Flow 2 Calculatos. For the referece bus #, we set : V = p.u. ad δ = 0 For all other
More informationarxiv:math/ v1 [math.gm] 8 Dec 2005
arxv:math/05272v [math.gm] 8 Dec 2005 A GENERALIZATION OF AN INEQUALITY FROM IMO 2005 NIKOLAI NIKOLOV The preset paper was spred by the thrd problem from the IMO 2005. A specal award was gve to Yure Boreko
More informationTaylor s Series and Interpolation. Interpolation & Curve-fitting. CIS Interpolation. Basic Scenario. Taylor Series interpolates at a specific
CIS 54 - Iterpolato Roger Crawfs Basc Scearo We are able to prod some fucto, but do ot kow what t really s. Ths gves us a lst of data pots: [x,f ] f(x) f f + x x + August 2, 25 OSU/CIS 54 3 Taylor s Seres
More informationCS5620 Intro to Computer Graphics
CS56 Itro to Computer Graphcs Geometrc Modelg art II Geometrc Modelg II hyscal Sples Curve desg pre-computers Cubc Sples Stadard sple put set of pots { } =, No dervatves specfed as put Iterpolate by cubc
More informationInvestigating Cellular Automata
Researcher: Taylor Dupuy Advsor: Aaro Wootto Semester: Fall 4 Ivestgatg Cellular Automata A Overvew of Cellular Automata: Cellular Automata are smple computer programs that geerate rows of black ad whte
More informationFibonacci Identities as Binomial Sums
It. J. Cotemp. Math. Sceces, Vol. 7, 1, o. 38, 1871-1876 Fboacc Idettes as Bomal Sums Mohammad K. Azara Departmet of Mathematcs, Uversty of Evasvlle 18 Lcol Aveue, Evasvlle, IN 477, USA E-mal: azara@evasvlle.edu
More informationDIFFERENTIAL GEOMETRIC APPROACH TO HAMILTONIAN MECHANICS
DIFFERENTIAL GEOMETRIC APPROACH TO HAMILTONIAN MECHANICS Course Project: Classcal Mechacs (PHY 40) Suja Dabholkar (Y430) Sul Yeshwath (Y444). Itroducto Hamltoa mechacs s geometry phase space. It deals
More informationLecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model
Lecture 7. Cofdece Itervals ad Hypothess Tests the Smple CLR Model I lecture 6 we troduced the Classcal Lear Regresso (CLR) model that s the radom expermet of whch the data Y,,, K, are the outcomes. The
More informationG S Power Flow Solution
G S Power Flow Soluto P Q I y y * 0 1, Y y Y 0 y Y Y 1, P Q ( k) ( k) * ( k 1) 1, Y Y PQ buses * 1 P Q Y ( k1) *( k) ( k) Q Im[ Y ] 1 P buses & Slack bus ( k 1) *( k) ( k) Y 1 P Re[ ] Slack bus 17 Calculato
More informationNumerical Solutions of Second Order Initial Value Problems of Bratu-Type equation Using Higher Ordered Rungu-Kutta Method
Iteratoal Joural of Scetfc ad Research Publcatos Volume 7 Issue October 7 87 ISSN 5-353 Numercal Solutos of Secod Order Ital Value Problems of Bratu-Tpe equato Usg Hgher Ordered Rugu-Kutta Method Habtamu
More informationThe internal structure of natural numbers, one method for the definition of large prime numbers, and a factorization test
Fal verso The teral structure of atural umbers oe method for the defto of large prme umbers ad a factorzato test Emmaul Maousos APM Isttute for the Advacemet of Physcs ad Mathematcs 3 Poulou str. 53 Athes
More informationResearch Article Gauss-Lobatto Formulae and Extremal Problems
Hdaw Publshg Corporato Joural of Iequaltes ad Applcatos Volume 2008 Artcle ID 624989 0 pages do:055/2008/624989 Research Artcle Gauss-Lobatto Formulae ad Extremal Problems wth Polyomals Aa Mara Acu ad
More informationSummary of the lecture in Biostatistics
Summary of the lecture Bostatstcs Probablty Desty Fucto For a cotuos radom varable, a probablty desty fucto s a fucto such that: 0 dx a b) b a dx A probablty desty fucto provdes a smple descrpto of the
More informationONE GENERALIZED INEQUALITY FOR CONVEX FUNCTIONS ON THE TRIANGLE
Joural of Pure ad Appled Mathematcs: Advaces ad Applcatos Volume 4 Number 205 Pages 77-87 Avalable at http://scetfcadvaces.co. DOI: http://.do.org/0.8642/jpamaa_7002534 ONE GENERALIZED INEQUALITY FOR CONVEX
More information1. A real number x is represented approximately by , and we are told that the relative error is 0.1 %. What is x? Note: There are two answers.
PROBLEMS A real umber s represeted appromately by 63, ad we are told that the relatve error s % What s? Note: There are two aswers Ht : Recall that % relatve error s What s the relatve error volved roudg
More informationOn the construction of symmetric nonnegative matrix with prescribed Ritz values
Joural of Lear ad Topologcal Algebra Vol. 3, No., 14, 61-66 O the costructo of symmetrc oegatve matrx wth prescrbed Rtz values A. M. Nazar a, E. Afshar b a Departmet of Mathematcs, Arak Uversty, P.O. Box
More informationOnline Publication Date: 12 December, 2011 Publisher: Asian Economic and Social Society
Ole Publcato Date: December, Publsher: Asa Ecoomc ad Socal Socety Soluto Of A System Of Two Partal Dfferetal Equatos Of The Secod Order Usg Two Seres Hayder Jabbar Abood (Departmet of Mathematcs, College
More information18.413: Error Correcting Codes Lab March 2, Lecture 8
18.413: Error Correctg Codes Lab March 2, 2004 Lecturer: Dael A. Spelma Lecture 8 8.1 Vector Spaces A set C {0, 1} s a vector space f for x all C ad y C, x + y C, where we take addto to be compoet wse
More informationIntroduction to Probability
Itroducto to Probablty Nader H Bshouty Departmet of Computer Scece Techo 32000 Israel e-mal: bshouty@cstechoacl 1 Combatorcs 11 Smple Rules I Combatorcs The rule of sum says that the umber of ways to choose
More informationPTAS for Bin-Packing
CS 663: Patter Matchg Algorthms Scrbe: Che Jag /9/00. Itroducto PTAS for B-Packg The B-Packg problem s NP-hard. If we use approxmato algorthms, the B-Packg problem could be solved polyomal tme. For example,
More informationUnimodality Tests for Global Optimization of Single Variable Functions Using Statistical Methods
Malaysa Umodalty Joural Tests of Mathematcal for Global Optmzato Sceces (): of 05 Sgle - 5 Varable (007) Fuctos Usg Statstcal Methods Umodalty Tests for Global Optmzato of Sgle Varable Fuctos Usg Statstcal
More informationA Study on Generalized Generalized Quasi hyperbolic Kac Moody algebra QHGGH of rank 10
Global Joural of Mathematcal Sceces: Theory ad Practcal. ISSN 974-3 Volume 9, Number 3 (7), pp. 43-4 Iteratoal Research Publcato House http://www.rphouse.com A Study o Geeralzed Geeralzed Quas (9) hyperbolc
More informationA Markov Chain Competition Model
Academc Forum 3 5-6 A Marov Cha Competto Model Mchael Lloyd, Ph.D. Mathematcs ad Computer Scece Abstract A brth ad death cha for two or more speces s examed aalytcally ad umercally. Descrpto of the Model
More informationGeneralized Convex Functions on Fractal Sets and Two Related Inequalities
Geeralzed Covex Fuctos o Fractal Sets ad Two Related Iequaltes Huxa Mo, X Su ad Dogya Yu 3,,3School of Scece, Bejg Uversty of Posts ad Telecommucatos, Bejg,00876, Cha, Correspodece should be addressed
More informationUNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS
UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS Postpoed exam: ECON430 Statstcs Date of exam: Jauary 0, 0 Tme for exam: 09:00 a.m. :00 oo The problem set covers 5 pages Resources allowed: All wrtte ad prted
More information5 Short Proofs of Simplified Stirling s Approximation
5 Short Proofs of Smplfed Strlg s Approxmato Ofr Gorodetsky, drtymaths.wordpress.com Jue, 20 0 Itroducto Strlg s approxmato s the followg (somewhat surprsg) approxmato of the factoral,, usg elemetary fuctos:
More informationMultiple Choice Test. Chapter Adequacy of Models for Regression
Multple Choce Test Chapter 06.0 Adequac of Models for Regresso. For a lear regresso model to be cosdered adequate, the percetage of scaled resduals that eed to be the rage [-,] s greater tha or equal to
More information