Research Article Numerical Solutions of the Second-Order One-Dimensional Telegraph Equation Based on Reproducing Kernel Hilbert Space Method

Size: px
Start display at page:

Download "Research Article Numerical Solutions of the Second-Order One-Dimensional Telegraph Equation Based on Reproducing Kernel Hilbert Space Method"

Transcription

1 Hndaw Publshng Corporaton Abstract and Appled Analyss Volume 3, Artcle ID , 3 pages Research Artcle Numercal Solutons of the Second-Order One-Dmensonal Telegraph Equaton Based on Reproducng Kernel Hlbert Space Method Mustafa Inc, Al Akgül,,3 and Adem KJlJçman 4 Department of Mathematcs, Scence Faculty, Fırat Unversty, 39 Elazığ, Turkey Department of Mathematcs, Educaton Faculty, Dcle Unversty, 8 Dyarbakır, Turkey 3 Department of Mathematcs and Statstcs, Mssour Unversty of Scence and Technology, Rolla, MO 6549-, USA 4 Department of Mathematcs and Insttute for Mathematcal Research, Unverst Putra Malaysa, 434Serdang,Selangor,Malaysa Correspondence should be addressed to Adem Kılıçman; klcman@yahoo.com Receved 8 December ; Accepted 5 March 3 Academc Edtor: Mustafa Bayram Copyrght 3 Mustafa Inc et al. Ths s an open access artcle dstrbuted under the Creatve Commons Attrbuton Lcense, whch permts unrestrcted use, dstrbuton, and reproducton n any medum, provded the orgnal work s properly cted. We nvestgate the effectveness of reproducng kernel method (RKM) n solvng partal dfferental equatons. We propose a reproducng kernel method for solvng the telegraph equaton wth ntal and boundary condtons based on reproducng kernel theory. Its exact soluton s represented n the form of a seres n reproducng kernel Hlbert space. Some numercal examples are gven n order to demonstrate the accuracy of ths method. The results obtaned from ths method are compared wth the exact solutons and other methods. Results of numercal examples show that ths method s smple, effectve, and easy to use.. Introducton The hyperbolc partal dfferental equatons model the vbratons of structures (e.g., buldngs, beams, and machnes). These equatons are the bass for fundamental equatons of atomc physcs. In ths paper, we consder the telegraph equaton of the form u (x, t) +α u t t (x, t) +β u (x, t) u (x, t) +f(x, t), x x,t,α>β, wth ntal condtons u (x, ) φ (x), u t (x, ) φ (x), () and approprate boundary condtons u (, t) g (t), u(, t) g (t), t (3) by usng reproducng kernel method (RKM). In recent years, much attenton has been gven n the lterature to () the development, analyss, and mplementaton of stable methods for the numercal soluton of () (3)[ 3]. Mohanty carred out a new technque to solve the lnear one-spacedmensonal hyperbolc equaton ()[4]. Hgh-order accurate method for solvng lnear hyperbolc equaton s presented n [5]. A compact fnte dfference approxmaton of fourth order for dscretzng spatal dervatve of lnear hyperbolc equaton and a collocaton method for the tme component are used n[6]. A numercal scheme s developed to solve the one-dmensonal hyperbolc telegraph equaton usng the collocaton ponts and approxmatng the soluton usng thn plate splnes radal bass functon [7]. Several test problems were gven, and the results of numercal experments were compared wth analytcal solutons to confrm the good accuracy of ther scheme. Yao [8] nvestgated a nonlnear hyperbolc telegraph equaton wth an ntegral condton by reproducng kernel space at α β. Yousef presented a numercal method for solvng the one-dmensonal hyperbolc telegraph equaton by usng Legendre multwavelet Galerkn method [9]. Dehghan and Lakestan presented a numercal technque for the soluton of the second-order

2 Abstract and Appled Analyss one-dmensonal lnear hyperbolc equaton []. Lakestan and Saray used nterpolatng scalng functons for solvng () (3) []. Dehghan provded a soluton of the second-order one-dmensonal hyperbolc telegraph equaton by usng the dualrecproctyboundaryntegralequaton(drbie)method [].Theproblemhasexplctsolutonthatcanbeobtanedby the method of separaton of varables n [3]. In ths paper, the problem s solved easly and elegantly by usng RKM. The technque has many advantages over the classcal technques. It also avods dscretzaton and provdes an effcent numercal soluton wth hgh accuracy, mnmal calculaton, and avodance of physcally unrealstc assumptons. In the next secton, we wll descrbe the procedure of ths method. The theory of reproducng kernels was used for the frst tme at the begnnng of the th century by Zaremba n hs work on boundary value problems for harmonc and bharmonc functons [4].Reproducng kernel theory has mportant applcaton n numercal analyss, dfferental equatons, probablty, and statstcs. Recently, usng the RKM, some authors dscussed telegraph equaton [5], Troesch s porblem [6], MHD Jeffery-Hamel flow [7], Bratu s problem [8], KdV equaton [9], fractonal dfferental equaton [], nonlnear oscllator wth dscontnuty [], nonlnear twopont boundary value problems [], ntegral equatons [3], and nonlnear partal dfferental equatons [4]. The paper s organzed as follows. Secton ntroduces several reproducng kernel spaces. The representaton n W(Ω) and a lnear operator are presented n Secton 3. Secton 4 provdes the man results. The exact and approxmate solutons of () (3) and an teratve method are developed for the knd of problems n the reproducng kernel space. We have proved that the approxmate soluton converges to the exact soluton unformly. Numercal experments are llustrated n Secton 5. Some conclusons are gven n Secton 6.. Reproducng Kernel Spaces In ths secton, some useful reproducng kernel spaces are defned. Defnton (reproducng kernel functon). Let E. A functon K : E E C s called areproducngkernel functon of the Hlbert space H f and only f (a) K(, t) H for all t E, (b) φ, K(, t) φ(t) for all t Eand all φ H. The last condton s called the reproducng property as the value of the functon φ at the pont t s reproduced by the nner product of φ wth K(, t). Defnton. Hlbert functon space H s a reproducng kernel space f and only f for any fxed x X,thelnear functonal I(f) f(x) s bounded [5,page5]. Defnton 3. We defne the space H [, ] by H [, ] {u AC[, ] :u L [, ]}. (4) The nner product and the norm n H [, ] are defned by u, g H u() g () + u (t) g (t) dt, u, g H [, ], u H u, u H, u H [, ]. (5) Lemma 4. The space H [, ] s a reproducng kernel space, and ts reproducng kernel functon q s s gven by [5,page3] q s (t) { +t, t s, + s, t > s. Defnton 5. We defne the space F 3 [, T] by F 3 [, T] {u AC[, T] :u,u AC[, T], u (3) L [, T],u() u () }. The nner product and the norm n F 3 [, T] are defned by u, g F 3 u () () g () () T + u (3) (t) g (3) (t) dt, u, g F 3 [, T], u F 3 u, u F 3, u F 3 [, T]. Lemma 6. The space F 3 [, T] s a reproducng kernel space, and ts reproducng kernel functon r s s gven by [5,page48] 4 s t + s t 3 4 st4 + t5, { t s, r s (t) { 4 s t + s3 t 4 ts4 + s5, { t>s. Defnton 7. We defne the space W 3 [, ] by W 3 [, ] {u AC[, ] :u,u AC[, ], u (3) L [, ],u() u() }. The nner product and the norm n W 3 [, ] are defned by u, g W 3 u () () g () () (6) (7) (8) (9) () + u (3) (x) g (3) (x) dx, u, g W 3 [, ], u W 3 u, u W 3, u W 3 [, ]. () The space W 3 [, ] s a reproducng kernel space, and ts reproducng kernel functon R y s gven by the followng theorem.

3 Abstract and Appled Analyss 3 Theorem 8. The space W 3 [, ] s a reproducng kernel space, and ts reproducng kernel functon R y s gven by c (y) x, { R y (x) 5 { d (y) x, { 5 x y, x>y, c (y) 5 56 y4 56 y5 5 6 y 5 78 y y, c (y) 5 64 y4 64 y5 + 4 y 5 3 y3 5 6 y, c 3 (y) 5 87 y4 87 y y y y, () c 4 (y) y y y y3 4 y, 5 c 5 (y) 3744 y4 87 y5 64 y 87 y3 56 y+, d (y) y5, d (y) 4 y4 56 y5 5 6 y 5 78 y y, d (y) 7 4 y y4 64 y5 + 4 y 5 6 y, d 3 (y) 5 87 y4 87 y5 5 3 y y y, d 4 (y) y y y y y, d 5 (y) 3744 y4 87 y5 64 y 87 y3 56 y. (3) Proof. Let u W 3 [, ] and let y.bydefnton 7 and ntegratng by parts two tmes, we obtan that u, R y W 3 u () () R () y () + u (3) (x) R (3) y (x) dx u() R y () +u () R y () +u () R y () +u () R (3) y () u () R (3) y () u () R (4) y () +u () R (4) y () + u (x) R (5) y (x) dx. (4) After substtutng the values of R y (), R y (), R y (), R(3) y (), R (4) y (), R(3) y (),andr(4) y () nto the above equaton, we get u, R y W 3 u() +u () ( 3 3 y 5 78 y3 5 6 y y4 56 y5 ) +u () ( 5 3 y 5 56 y3 + 5 y y4 3 y5 ) +u () u () ( 5 3 y 5 56 y3 + 5 y y4 3 y5 ) u () +u () ( 3 3 y y y 5 56 y y5 ) + u (x) R (5) y (x) dx; (5) thus we obtan that u, R y W 3 u (x) R (5) y (x) dx y u (x) R (5) y (x) dx+ y u (x) R (5) y (x) dx y u (x) ( 3 y 5 78 y3 5 6 y y4 56 y5 ) dx + u (x) ( y 3 y 5 78 y3 5 6 y y4 56 y5 ) dx (u(y) u()) ( 3 y 5 78 y3 5 6 y y4 56 y5 )

4 4 Abstract and Appled Analyss +(u() u (y)) ( 3 y 5 78 y3 5 6 y y4 56 y5 ) u(y) u() ( 3 y 5 78 y3 5 6 y y4 56 y5 ) +u() ( 3 y 5 78 y3 5 6 y y4 56 y5 ). By Defnton 7,wehaveu() u().so Ths completes the proof. (6) u, R y W 3 u(y). (7) Defnton 9. We defne the bnary space W(Ω) by W (Ω) {u : 4 u x s completely contnuous n t Ω[, ] [, ], u (x, ), u (x, ) t 6 u x 3 t 3 L (Ω),, u (, t),u(, t) }. The nner product and the norm n W(Ω) are defned by u, g W + [ 3 t 3 j j u (, t) x 3 t 3 g (, t)] dt x t j u (x, ), j t j g (x, ) W 3 + [ 3 3 x 3 u (x, t) t3 3 3 x 3 g (x, t)] dt dx, t3 (8) (9) u W u, u W, u W(Ω). () Lemma. W(Ω) s a reproducng kernel space, and ts reproducng kernel functon K (y,s) s gven by [5,page48] K (y,s) R y r s. () Defnton. We defne the bnary space W(Ω) by W (Ω) {u : u s completely contnuous n Ω[, ] [, ], u x t L (Ω)}. The nner product and the norm n W(Ω) are defned by u, g W [ t u (, t) g (, t)] dt t u W u, u W, u W (Ω). () + u(x, ),g(x, ) H + [ (3) u (x, t) g (x, t)] dt dx, x t x t Lemma. W(Ω) s a reproducng kernel space, and ts reproducng kernel functon G (y,s) s gven as [5,page48] G (y,s) Q y q s. (4) Remark 3. Hlbert spaces can be completely classfed: there s a unque Hlbert space up to somorphsm for every cardnalty of the base. Snce fnte-dmensonal Hlbert spaces are fully understood n lnear algebra and snce morphsms of Hlbert spaces can always be dvded nto morphsms of spaces wth Aleph-null (χ ) dmensonalty, functonal analyss of Hlbert spaces mostly deals wth the unque Hlbert space of dmensonalty Aleph-null and ts morphsms. One of the open problems n functonal analyss stoprovethateveryboundedlnearoperatoronahlbert space has a proper nvarant subspace. Many specal cases of ths nvarant subspace problem have already been proven [6]. 3. Soluton Representaton n W(Ω) In ths secton, the soluton of () s gven n the reproducng kernel space W(Ω).On defnng the lnear operatorl : W(Ω) W(Ω) by LV V t V x +α V t +β V (x, t), (5) after homogenzng the ntal and boundary condtons, model problem () (3) changes to the problem LV M(x, t), (x, t) [, ] [, ], V (x, ) V t (x, ) V (, t) V (, t), (6)

5 Abstract and Appled Analyss 5 M (x, t) Z t (x, t) V (x, t) x +α Z t (x, t) +β Z (x, t) +f(x, t) ; for convenence, we agan wrte u nstead of V n (6) Lemma 4. L s a bounded lnear operator. (7) Now, choose a countable dense subset {(x,t ), (x, t ),...}n Ω and defne φ G (x,t ), Ψ L φ, (3) L s the adjont operator of L.Theorthonormalsystem { Ψ } of W(Ω) canbedervedfromtheprocessofgram- Schmdt orthogonalzaton of {Ψ } as Proof. Let u W(Ω) and let (x, t) Ω. ByLemma, we have and thus u (x, t) u,k (x,t) W, (8) Lu (x, t) u,lk (x,t) W, Ψ k β k Ψ k. (33) Theorem 5. Suppose that {(x,t )} s dense n Ω. Then {Ψ } s a complete system n W(Ω),and x Lu (x, t) u, x LK, (x,t) W t Lu (x, t) u, t LK, (x,t) W t x Lu (x, t) u, t x LK. (x,t) W Hence there exst a,b,a,b >such that Therefore, Lu (x, t) a u W, t Lu (x, t) b u W, x Lu (x, t) a u W, t x Lu (x, t) b u W. Lu W [ dt t Lu (, t)] + Lu (x, ),Lu(x, ) H + [ x Ths completes the proof. t Lu (x, t)] dt dx [ t Lu (, t)] dt+[lu (, )] + [ x Lu (x, )] dx + [ x t Lu (x, t)] dt dx (a +a +b +b ) u W. (9) (3) (3) Ψ LK (x,t ) (x, t). (34) Proof. We have Ψ L φ L φ,k (x,t) W φ,lk (x,t) W LK (x,t),g (x,t (35) ) W LK (x,t) (x,t ) LK (x,t ) (x, t). Clearly, Ψ W(Ω).Foreachfxedu W(Ω),f u, Ψ W,,,..., (36) then u,ψ W u,l φ W (37) Lu,φ W Lu(x,t ),,,... Note that {(x,t )} s dense n Ω.Hence,Lu.Fromthe exstence of L,tfollowsthatu.Theproofscompleted. Theorem 6. If {(x,t )} s dense n Ω, then the soluton of (6) s gven by u β k M(x k,t k ) Ψ. (38) k

6 6 Abstract and Appled Analyss Proof. By Theorem 5, {Ψ (x, t)} s a complete system n W(Ω).Thus, u u, Ψ W Ψ β k u, Ψ k W Ψ k u(x) 5 5 β k u, L φ k W Ψ k β k Lu, φ k W Ψ k β k Lu, G (xk,t k ) W Ψ k (39) 4 RKM ES 4 Fgure : Graph of numercal results for Example (α, β, t.5). x β k Lu (x k,t k ) Ψ k β k M(x k,t k ) Ψ. k Ths completes the proof. Now the approxmate soluton u n canbeobtanedfrom the n-term ntercept of the exact soluton u and Obvously n u n β k M(x k,t k ) Ψ. (4) k u n u W, n. (4) Theorem 7. If u W(Ω),then u n u W, n. (4) Moreover, a sequence u n u W s monotoncally decreasng n n. Proof. From (38) and(4), t follows that Thus, u n u W n+ β k M(x k,t k ) Ψ (x,. (43) t) W k u n (x, t) u(x, t) W, n. (44) In addton u n u W n+ n+ ( β k M(x k,t k ) Ψ (x, t) k k W β k M(x k,t k ) Ψ (x, t)). (45) Clearly, u n (x, t) u(x, t) W s monotoncally decreasng n n. 4. The Method Implementaton () If (6) s lnear, then the analytcal soluton of (6) can be obtaned drectly by (38). () If (6) s nonlnear, then the soluton of (6) canbe obtaned by the followng teratve method. We construct an teratve sequence u n,puttng A A n any fxed u W 3 [, ], u n n A Ψ, A β M(x k,t k,u(x k,t k )) k n k β k M(x k,t k,u k (x k,t k )),. β nk M(x k,t k,u k (x k,t k )). (46) (47) Next we wll prove that u n gven by the teratve formula (46) converges to the exact soluton. Theorem 8. Suppose that u n defned by (46) s bounded and (6) has a unque soluton. If {(x,t )} s dense n Ω, then u n converges to the analytcal soluton u of (6),and A s gven by (47). u A Ψ, (48)

7 Abstract and Appled Analyss 7 Table : Numercal results for Example for t.5. RE RE RE x CPU tme (s) α α α α α 5 α 5 β β β5 β5 β β Table : Numercal results for Example for t.. RE RE RE x CPU tme (s) α α α α α 5 α 5 β β β5 β5 β β Table 3: A comparson between nterpolatng scalng functon method [] and RKM for dfferent values of α, β, and t for Example. x CPU tme (s) [] α t.5 RKM α t.5 [] β t RKM β t Table 4: A comparson between nterpolatng scalng functon method [] and RKM for dfferent values of α, β, and t for Example. x [] α t.5 RKM α t.5 [] β5 t RKM β5 t Proof. Frst,weprovetheconvergenceofu n.from(46) and the orthonormalty of { Ψ },wenferthat u n+ n+ A n A +A n+ u n +A n+ u n. (49) By (49), u n s nondecreasng, and by the assumpton, u n s bounded. Thus u n s convergent. By (49), there exsts a constant c such that A c. (5) Ths mples that If m>n,then u m u n k+ u k knu {A } l. (5) m m u k+ u k kn m A k+ kn, m, n. (5)

8 8 Abstract and Appled Analyss Table 5: Numercal results for Example for t.5. RE RE RE x CPU tme (s) α α α α α5 α5 β β β5 β5 β β Table 6: Numercal results for Example for t.. RE RE RE x CPU tme (s) α α α α α5 α5 β β β5 β5 β β u(x) RKM ES 4 6 Fgure : Graph of numercal results for Example (α, β, t.5). The completeness of W(Ω) shows that there exsts u W(Ω) such that u n u as n.now,weprovethat u solves (6). Takng lmts n (4), we get Note that u L u (L u) (x k,t k ) x A Ψ. (53) A L Ψ, (54) A L Ψ (x k,t k ) A L Ψ,G (xk,t k ) W In vew of (47), we have A L Ψ,φ k W A Ψ,L φ k W A Ψ,Ψ k W. (55) L u (x l,t l ) M(x l,t l,u l (x l,t l )). (56) Snce {(x,t )} s dense n Ω,foreach(y, s) Ω, there exsts asubsequence{(x nj,t nj )} j such that We know that (x nj,t nj ) (y, s), j. (57) L u(x nj,t nj )M(x nj,t nj,u nj (x nj,t nj )). (58) Let j.bythecontnutyoff,wehave (L u) (y, s) M(y, s, u (y, s)), (59) whch ndcates that u satsfes (6). Remark 9. Inthesamemanner,tcanbeprovedthat u n x u, n, (6) x u x Ψ A x, and A s gven by (47). n u n x Ψ A x, (6)

9 Abstract and Appled Analyss 9 Table 7: RMS errors for Example. [] RKM [] RKM [] RKM N α α α α α α β β β β β5 β5 t.5 t.5 t. t. t.5 t Table 8: Numercal results for Example for t.5. RE RE RE x α α CPU tme (s) α α α 5 α 5 β β β5 β5 β β Numercal Results To test the accuracy of the present method, some numercal experments are presented n ths secton. Usng our method, we chose 36 ponts n Ω and obtaned the approxmate soluton u 36. The comparson between nterpolatng scalng functon method [] and RKM for dfferent values of α, β, and t sgvenntables4 and 5. We solve these examples for a set of ponts {x,...,x ( ) h,...,x N }, h N. (6) In Tables 7 and we calculate the RMS error by the followng formula: RMS error N+ (u u 36). (63) N+ ItcanbeseenfromTables4 and 5 and 7 that the results obtaned by the RKM are more accurate than those obtaned by the methods n [, ]. Ths ndcates that RKM s a relable method. The CPU tme (s) s gven n Tables,, 3, 4, 5, 6, 7, 8, 9,and. Numercal solutons are descrbed n the extended doman [ 4, 4] [ 3, 3]. ThecomparsonofRMSerrors gven for our method and Chebyshev method. Example. Consder the followng telegraph equaton wth ntal and boundary condtons: u (x, t) +α u t t (x, t) +β u (x, t) u +f(x, t), x u (x, ) φ (x) snh (x), u t (x, ) φ (x) snh (x), u (, t) g (t), t, u (, t) g (t) exp ( t) snh (), t, (64) f (x, t) (3 4α+β ) exp ( t) snh (x). (65) The exact soluton of (64)s gven by[] If we apply u (x, t) exp ( t) snh (x). (66) V (x, t) u(x, t) xsnh () (exp ( t) +t) + snh (x)(t ) to (64), then the followng equaton (68)sobtaned: V t V x +α V t +β V (x, t) M(x, t), V (x, ) V (x, ) V (, t) V (, t), t M (x, t) (3 4α+β ) exp ( t) snh (x) snh (x)(t ) αsnh () ( exp ( t) +)x +4αsnh (x) β snh () (exp ( t) +t)x (67) (68) +β snh (x)(t ) 4xexp ( t) snh (). (69)

10 Abstract and Appled Analyss Table 9: Numercal results for Example for t.. RE RE RE x CPU tme (s) α α α α α 5 α 5 β β β5 β5 β β Table : RMS errors for Example. [] RKM [] RKM [] RKM N α α α5 α5 α5 α5 β5 β5 β β β β t. t. t.5 t.5 t. t After homogenzng the ntal and boundary condtons and usng the above method, we obtan Tables 4 and Fgure. Example. Consder the followng telegraph equaton wth ntal and boundary condtons: u t +α u t +β u (x, t) u +f(x, t), x u (x, ) φ (x) sn (x), u t (x, ) φ (x), u (, t) g (t), t, u (, t) g (t) cos (t) sn (), t, (7) f (x, t) αsn (t) sn (x) +β cos (t) sn (x). (7) The exact soluton of (7)s gven by[] If we apply u (x, t) cos (t) sn (x). (7) V (x, t) u(x, t) xsn ()(cos (t) ) sn (x) (73) to (7), thenthefollowng (74)sobtaned: V t (x, t) V (x, t) +α V x t +β V (x, t) M(x, t), V (x, ) V (x, ) V (, t) V (, t), t (74) M (x, t) αsn (t) sn (x) +β cos (t) sn (x) sn (x) αsn (t) snh () x+xcos (t) sn () β x sn ()(cos (t) ) β sn (x). (75) After homogenzng the ntal and boundary condtons and usng the above method, we obtan Tables 5 7 and Fgures and 3. Example. Consder the followng telegraph equaton wth ntal and boundary condtons: u t +α u t +β u (x, t) u +f(x, t), x u (x, ) φ (x) tan ( x ), u t (x, ) φ (x) u (, t) g (t) tan ( t ), ( + tan (x/)), t, u (, t) g (t) tan ( +t ), t, (76) f (x, t) α(+tan ( x+t )) +β tan ( x+t ). (77) The exact soluton of (76)sgvenby[] u (x, t) tan ( x+t ). (78)

11 Abstract and Appled Analyss (a) (b) Fgure 3: Plots of numercal results for Example (α 6,β) If we apply V (x, t) u(x, t) + (x ) tan ( t ) xtan ( +t ) tan (x ) (+( +xtan ( t tan (x/) )) tan (/) )(+(t )) to (76), then the followng equaton (8)sobtaned V t (x, t) V (x, t) +α V x t +β V (x, t) M(x, t), V (x, ) V (x, ) V (, t) V (, t), t M (x, t) α(+tan ( x+t )) + β tan ( x+t ) tan (x/) +( + t 4 + 3t 4 tan ( x )) (+tan ( x )) + (x ) α(+tan ( t )) αx(+tan ( t+ )) α tan ( x ) +αxtan ( )+β (x ) tan ( t ) β x tan ( +t ) β tan ( x ) (79) (8) (+ +β x tan ( x x t tan (x/) ) t tan (/) )(+ ) tan ( t )(+tan ( t )) tan (+t )(+tan ( +t )). (8) After homogenzng the ntal and boundary condtons and usng the above method, we obtan Tables 8 and Fgure 4. Remark 3. In Tables 9, weabbrevatetheexactsoluton and the approxmate soluton by AS and ES, respectvely. stands for the absolute error, that s, the absolute value of the dfference of the exact soluton and the approxmate soluton, whlerendcatestherelatveerror,thats,theabsoluteerror dvded by the absolute value of the exact soluton. 6. Concluson In ths study, a second-order one-dmensonal telegraph equaton wth ntal and boundary condtons was solved by reproducng kernel Hlbert space method. We descrbed the method and used t n some test examples n order to show ts applcablty and valdty n comparson wth exact and other numercal solutons. The obtaned results show that ths approach can solve the problem effectvely and need few computatons. The satsfactory results that we obtaned were compared wth the results that were obtaned by [, ]. Numercal experments on test examples show that our proposed schemes are of hgh accuracy and support thetheoretcalresults.asshownntables7 and our results are better than the results that were obtaned by []. Accordng to these results, t s possble to apply RKM to lnear and nonlnear dfferental equatons wth ntal and

12 Abstract and Appled Analyss (a) (b) Fgure 4: Plots of numercal results for Example (α, β 5). boundary condtons. It has been shown that the obtaned results are unformly convergent and the operator that was used s a bounded lnear operator. From the results, RKM can be appled to hgh dmensonal partal dfferental equatons, ntegral equatons, and fractonal dfferental equatons wthout any transformaton or dscretzaton, and good results can be obtaned. References [] M. Dehghan, On the soluton of an ntal-boundary value problem that combnes Neumann and ntegral condton for the wave equaton, Numercal Methods for Partal Dfferental Equatons,vol.,no.,pp.4 4,5. []R.K.Mohanty,M.K.Jan,andK.George, Ontheuseof hgh order dfference methods for the system of one space second order nonlnear hyperbolc equatons wth varable coeffcents, Journal of Computatonal and Appled Mathematcs,vol.7,no.,pp.4 43,996. [3] E. H. Twzell, An explct dfference method for the wave equaton wth extended stablty range, BIT, vol.9,no.3,pp , 979. [4] R. K. Mohanty, An uncondtonally stable dfference scheme for the one-space-dmensonal lnear hyperbolc equaton, Appled Mathematcs Letters,vol.7,no.,pp. 5,4. [5] A. Mohebb and M. Dehghan, Hgh order compact soluton of the one-space-dmensonal lnear hyperbolc equaton, Numercal Methods for Partal Dfferental Equatons,vol.4,no.5,pp. 35, 8. [6] M. Dehghan, Fnte dfference procedures for solvng a problem arsng n modelng and desgn of certan optoelectronc devces, Mathematcs and Computers n Smulaton,vol.7,no., pp. 6 3, 6. [7] M. Dehghan and A. Shokr, A numercal method for solvng the hyperbolc telegraph equaton, Numercal Methods for Partal Dfferental Equatons,vol.4,no.4,pp.8 93,8. [8] H. Yao, Reproducng kernel method for the soluton of nonlnear hyperbolc telegraph equaton wth an ntegral condton, Numercal Methods for Partal Dfferental Equatons,vol.7,no. 4, pp ,. [9] S. A. Yousef, Legendre multwavelet Galerkn method for solvng the hyperbolc telegraph equaton, Numercal Methods for Partal Dfferental Equatons, vol.6,no.3,pp ,. [] M. Dehghan and M. Lakestan, The use of Chebyshev cardnal functons for soluton of the second-order one-dmensonal telegraph equaton, Numercal Methods for Partal Dfferental Equatons,vol.5,no.4,pp ,9. [] M. Lakestan and B. N. Saray, Numercal soluton of telegraph equaton usng nterpolatng scalng functons, Computers & Mathematcs wth Applcatons, vol.6,no.7,pp ,. [] M. Dehghan and A. Ghesmat, Soluton of the second-order one-dmensonal hyperbolc telegraph equaton by usng the dual recprocty boundary ntegral equaton (DRBIE) method, Engneerng Analyss wth Boundary Elements,vol.34,no.,pp. 5 59,. [3] A. N. Tkhonov and A. A. Samarskĭ, Equatons of Mathematcal Physcs,Dover,NewYork,NY,USA,99. [4] N. Aronszajn, Theory of reproducng kernels, Transactons of the Amercan Mathematcal Socety,vol.68,pp ,95. [5] M. Inc, A. Akgül, and A. Kılıçman, Explct soluton of telegraph equaton based on reproducng kernel method, Journal of Functon Spaces and Applcatons,vol.,p.3,. [6] M. Inc and A. Akgül, The reproducng kernel Hlbert space method for solvng Troesch s problem, Journal of the Assocaton of Arab Unverstes for Basc and Appled Scences.Inpress. [7] M. Inc, A. Akgül, and A. Kılıçman, A new applcaton of the reproducng kernel Hlbert space method to solve MHD Jeffery- Hamel flows problem n nonparallel walls, Abstract and Appled Analyss, vol. 3, Artcle ID 39454, pages, 3. [8] M. Inc, A. Akgül, and F. Geng, Reproducng kernel Hlbert space method for solvng Bratu s problem, Bulletn of the Malaysan Mathematcal Scences Socety.Inpress. [9] M. Inc, A. Akgül, and A. Kılıçman, On solvng KdV equaton usng reproducng kernel Hlbert space method, Abstract and Appled Analyss, vol. 3, Artcle ID 57894, pages, 3. [] W. Jang and Y. Ln, Representaton of exact soluton for the tme-fractonal telegraph equaton n the reproducng kernel

13 Abstract and Appled Analyss 3 space, Communcatons n Nonlnear Scence and Numercal Smulaton,vol.6,no.9,pp ,. [] Y. Wang, L. Su, X. Cao, and X. L, Usng reproducng kernel for solvng a class of sngularly perturbed problems, Computers & Mathematcs wth Applcatons, vol. 6, no., pp. 4 43,. [] F. Geng and M. Cu, A novel method for nonlnear two-pont boundary value problems: combnaton of ADM and RKM, Appled Mathematcs and Computaton,vol.7,no.9,pp ,. [3] F. Geng and F. Shen, Solvng a Volterra ntegral equaton wth weakly sngular kernel n the reproducng kernel space, Mathematcal Scences Quarterly Journal, vol.4,no.,pp.59 7,. [4] F. Geng and M. Cu, Homotopy perturbaton-reproducng kernel method for nonlnear systems of second order boundary value problems, Journal of Computatonal and Appled Mathematcs, vol. 35, no. 8, pp. 45 4,. [5] M. Cu and Y. Ln, Nonlnear Numercal Analyss n the Reproducng Kernel Space, Nova Scence, New York, NY, USA, 9. [6] A. M. Krall, Hlbert Space, Boundary Value Problems and Orthogonal Polynomals, vol.33ofoperator Theory: Advances and Applcatons, Brkhäuser, Basel, Swtzerland,.

The Analytical Solution of a System of Nonlinear Differential Equations

The Analytical Solution of a System of Nonlinear Differential Equations Int. Journal of Math. Analyss, Vol. 1, 007, no. 10, 451-46 The Analytcal Soluton of a System of Nonlnear Dfferental Equatons Yunhu L a, Fazhan Geng b and Mnggen Cu b1 a Dept. of Math., Harbn Unversty Harbn,

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation Nonl. Analyss and Dfferental Equatons, ol., 4, no., 5 - HIKARI Ltd, www.m-har.com http://dx.do.org/.988/nade.4.456 Asymptotcs of the Soluton of a Boundary alue Problem for One-Characterstc Dfferental Equaton

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

Research Article Cubic B-Spline Collocation Method for One-Dimensional Heat and Advection-Diffusion Equations

Research Article Cubic B-Spline Collocation Method for One-Dimensional Heat and Advection-Diffusion Equations Appled Mathematcs Volume 22, Artcle ID 4587, 8 pages do:.55/22/4587 Research Artcle Cubc B-Splne Collocaton Method for One-Dmensonal Heat and Advecton-Dffuson Equatons Joan Goh, Ahmad Abd. Majd, and Ahmad

More information

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

Research Article Green s Theorem for Sign Data

Research Article Green s Theorem for Sign Data Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

Sharp integral inequalities involving high-order partial derivatives. Journal Of Inequalities And Applications, 2008, v. 2008, article no.

Sharp integral inequalities involving high-order partial derivatives. Journal Of Inequalities And Applications, 2008, v. 2008, article no. Ttle Sharp ntegral nequaltes nvolvng hgh-order partal dervatves Authors Zhao, CJ; Cheung, WS Ctaton Journal Of Inequaltes And Applcatons, 008, v. 008, artcle no. 5747 Issued Date 008 URL http://hdl.handle.net/07/569

More information

Fixed point method and its improvement for the system of Volterra-Fredholm integral equations of the second kind

Fixed point method and its improvement for the system of Volterra-Fredholm integral equations of the second kind MATEMATIKA, 217, Volume 33, Number 2, 191 26 c Penerbt UTM Press. All rghts reserved Fxed pont method and ts mprovement for the system of Volterra-Fredholm ntegral equatons of the second knd 1 Talaat I.

More information

General viscosity iterative method for a sequence of quasi-nonexpansive mappings

General viscosity iterative method for a sequence of quasi-nonexpansive mappings Avalable onlne at www.tjnsa.com J. Nonlnear Sc. Appl. 9 (2016), 5672 5682 Research Artcle General vscosty teratve method for a sequence of quas-nonexpansve mappngs Cuje Zhang, Ynan Wang College of Scence,

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

Numerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method

Numerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method Appled Mathematcs, 6, 7, 5-4 Publshed Onlne Jul 6 n ScRes. http://www.scrp.org/journal/am http://.do.org/.436/am.6.77 umercal Solutons of a Generalzed th Order Boundar Value Problems Usng Power Seres Approxmaton

More information

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k) ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of

More information

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b Int J Contemp Math Scences, Vol 3, 28, no 17, 819-827 A New Refnement of Jacob Method for Soluton of Lnear System Equatons AX=b F Naem Dafchah Department of Mathematcs, Faculty of Scences Unversty of Gulan,

More information

Convexity preserving interpolation by splines of arbitrary degree

Convexity preserving interpolation by splines of arbitrary degree Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY. Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

Research Article The Solution of Two-Point Boundary Value Problem of a Class of Duffing-Type Systems with Non-C 1 Perturbation Term

Research Article The Solution of Two-Point Boundary Value Problem of a Class of Duffing-Type Systems with Non-C 1 Perturbation Term Hndaw Publshng Corporaton Boundary Value Problems Volume 9, Artcle ID 87834, 1 pages do:1.1155/9/87834 Research Artcle The Soluton of Two-Pont Boundary Value Problem of a Class of Duffng-Type Systems wth

More information

Case Study of Markov Chains Ray-Knight Compactification

Case Study of Markov Chains Ray-Knight Compactification Internatonal Journal of Contemporary Mathematcal Scences Vol. 9, 24, no. 6, 753-76 HIKAI Ltd, www.m-har.com http://dx.do.org/.2988/cms.24.46 Case Study of Marov Chans ay-knght Compactfcaton HaXa Du and

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Solving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions

Solving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions ISSN 746-7659 England UK Journal of Informaton and Computng Scence Vol. 9 No. 3 4 pp. 69-8 Solvng Fractonal Nonlnear Fredholm Integro-dfferental Equatons va Hybrd of Ratonalzed Haar Functons Yadollah Ordokhan

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LIX, 013, f.1 DOI: 10.478/v10157-01-00-y ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION BY ION

More information

Report on Image warping

Report on Image warping Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.

More information

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices Internatonal Mathematcal Forum, Vol 11, 2016, no 11, 513-520 HIKARI Ltd, wwwm-hkarcom http://dxdoorg/1012988/mf20166442 The Jacobsthal and Jacobsthal-Lucas Numbers va Square Roots of Matrces Saadet Arslan

More information

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity Int. Journal of Math. Analyss, Vol. 6, 212, no. 22, 195-114 Unqueness of Weak Solutons to the 3D Gnzburg- Landau Model for Superconductvty Jshan Fan Department of Appled Mathematcs Nanjng Forestry Unversty

More information

Haar wavelet collocation method to solve problems arising in induction motor

Haar wavelet collocation method to solve problems arising in induction motor ISSN 746-7659, England, UK Journal of Informaton and Computng Scence Vol., No., 07, pp.096-06 Haar wavelet collocaton method to solve problems arsng n nducton motor A. Padmanabha Reddy *, C. Sateesha,

More information

6.3.4 Modified Euler s method of integration

6.3.4 Modified Euler s method of integration 6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from

More information

The Finite Element Method: A Short Introduction

The Finite Element Method: A Short Introduction Te Fnte Element Metod: A Sort ntroducton Wat s FEM? Te Fnte Element Metod (FEM) ntroduced by engneers n late 50 s and 60 s s a numercal tecnque for solvng problems wc are descrbed by Ordnary Dfferental

More information

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space.

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space. Lnear, affne, and convex sets and hulls In the sequel, unless otherwse specfed, X wll denote a real vector space. Lnes and segments. Gven two ponts x, y X, we defne xy = {x + t(y x) : t R} = {(1 t)x +

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

Quantum Particle Motion in Physical Space

Quantum Particle Motion in Physical Space Adv. Studes Theor. Phys., Vol. 8, 014, no. 1, 7-34 HIKARI Ltd, www.-hkar.co http://dx.do.org/10.1988/astp.014.311136 Quantu Partcle Moton n Physcal Space A. Yu. Saarn Dept. of Physcs, Saara State Techncal

More information

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method

More information

A Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods

A Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods Appled Mathematcal Scences, Vol. 11, 2017, no. 52, 2579-2586 HIKARI Ltd, www.m-hkar.com https://do.org/10.12988/ams.2017.79280 A Soluton of the Harry-Dym Equaton Usng Lattce-Boltzmannn and a Soltary Wave

More information

Chapter 4 The Wave Equation

Chapter 4 The Wave Equation Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as

More information

2 Finite difference basics

2 Finite difference basics Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Group Analysis of Ordinary Differential Equations of the Order n>2

Group Analysis of Ordinary Differential Equations of the Order n>2 Symmetry n Nonlnear Mathematcal Physcs 997, V., 64 7. Group Analyss of Ordnary Dfferental Equatons of the Order n> L.M. BERKOVICH and S.Y. POPOV Samara State Unversty, 4430, Samara, Russa E-mal: berk@nfo.ssu.samara.ru

More information

A Solution of Porous Media Equation

A Solution of Porous Media Equation Internatonal Mathematcal Forum, Vol. 11, 016, no. 15, 71-733 HIKARI Ltd, www.m-hkar.com http://dx.do.org/10.1988/mf.016.6669 A Soluton of Porous Meda Equaton F. Fonseca Unversdad Naconal de Colomba Grupo

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

829. An adaptive method for inertia force identification in cantilever under moving mass

829. An adaptive method for inertia force identification in cantilever under moving mass 89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,

More information

Some modelling aspects for the Matlab implementation of MMA

Some modelling aspects for the Matlab implementation of MMA Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton

More information

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS

More information

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

Modelli Clamfim Equazione del Calore Lezione ottobre 2014 CLAMFIM Bologna Modell 1 @ Clamfm Equazone del Calore Lezone 17 15 ottobre 2014 professor Danele Rtell danele.rtell@unbo.t 1/24? Convoluton The convoluton of two functons g(t) and f(t) s the functon (g

More information

Existence of Two Conjugate Classes of A 5 within S 6. by Use of Character Table of S 6

Existence of Two Conjugate Classes of A 5 within S 6. by Use of Character Table of S 6 Internatonal Mathematcal Forum, Vol. 8, 2013, no. 32, 1591-159 HIKARI Ltd, www.m-hkar.com http://dx.do.org/10.12988/mf.2013.3359 Exstence of Two Conjugate Classes of A 5 wthn S by Use of Character Table

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

New numerical scheme for solving Troesch s Problem

New numerical scheme for solving Troesch s Problem Mathematcal Theory and Modelng ISSN 4-584 (Paper) ISSN 5-5 (Onlne) Vol4 No 4 wwwsteorg Abstract New numercal scheme for solvng Troesch s Problem Abdelmajd El hajaj * Khald Hlal El ghordaf Jalla and El

More information

Research Article Relative Smooth Topological Spaces

Research Article Relative Smooth Topological Spaces Advances n Fuzzy Systems Volume 2009, Artcle ID 172917, 5 pages do:10.1155/2009/172917 Research Artcle Relatve Smooth Topologcal Spaces B. Ghazanfar Department of Mathematcs, Faculty of Scence, Lorestan

More information

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning Journal of Machne Learnng Research 00-9 Submtted /0; Publshed 7/ Erratum: A Generalzed Path Integral Control Approach to Renforcement Learnng Evangelos ATheodorou Jonas Buchl Stefan Schaal Department of

More information

On the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros

On the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros Appled Mathematcal Scences, Vol. 5, 2011, no. 75, 3693-3706 On the Interval Zoro Symmetrc Sngle-step Procedure for Smultaneous Fndng of Polynomal Zeros S. F. M. Rusl, M. Mons, M. A. Hassan and W. J. Leong

More information

The internal structure of natural numbers and one method for the definition of large prime numbers

The internal structure of natural numbers and one method for the definition of large prime numbers The nternal structure of natural numbers and one method for the defnton of large prme numbers Emmanul Manousos APM Insttute for the Advancement of Physcs and Mathematcs 3 Poulou str. 53 Athens Greece Abstract

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

Integrals and Invariants of Euler-Lagrange Equations

Integrals and Invariants of Euler-Lagrange Equations Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

Canonical transformations

Canonical transformations Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,

More information

2nd International Conference on Electronics, Network and Computer Engineering (ICENCE 2016)

2nd International Conference on Electronics, Network and Computer Engineering (ICENCE 2016) nd Internatonal Conference on Electroncs, Network and Computer Engneerng (ICENCE 6) Postve solutons of the fourth-order boundary value problem wth dependence on the frst order dervatve YuanJan Ln, a, Fe

More information

On Finite Rank Perturbation of Diagonalizable Operators

On Finite Rank Perturbation of Diagonalizable Operators Functonal Analyss, Approxmaton and Computaton 6 (1) (2014), 49 53 Publshed by Faculty of Scences and Mathematcs, Unversty of Nš, Serba Avalable at: http://wwwpmfnacrs/faac On Fnte Rank Perturbaton of Dagonalzable

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

Power law and dimension of the maximum value for belief distribution with the max Deng entropy

Power law and dimension of the maximum value for belief distribution with the max Deng entropy Power law and dmenson of the maxmum value for belef dstrbuton wth the max Deng entropy Bngy Kang a, a College of Informaton Engneerng, Northwest A&F Unversty, Yanglng, Shaanx, 712100, Chna. Abstract Deng

More information

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

arxiv: v1 [math.co] 12 Sep 2014

arxiv: v1 [math.co] 12 Sep 2014 arxv:1409.3707v1 [math.co] 12 Sep 2014 On the bnomal sums of Horadam sequence Nazmye Ylmaz and Necat Taskara Department of Mathematcs, Scence Faculty, Selcuk Unversty, 42075, Campus, Konya, Turkey March

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009 College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:

More information

Research Article A Generalized Sum-Difference Inequality and Applications to Partial Difference Equations

Research Article A Generalized Sum-Difference Inequality and Applications to Partial Difference Equations Hndaw Publshng Corporaton Advances n Dfference Equatons Volume 008, Artcle ID 695495, pages do:0.55/008/695495 Research Artcle A Generalzed Sum-Dfference Inequalty and Applcatons to Partal Dfference Equatons

More information

Cubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two

Cubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two World Academy of Scence Engneerng and echnology Internatonal Journal of Mathematcal and omputatonal Scences Vol: No:0 00 ubc rgonometrc B-Splne Appled to Lnear wo-pont Boundary Value Problems of Order

More information

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 ISSN 1549-644 1 Scence Publcatons do:1.844/jmssp.1.4.8 Publshed Onlne 9 (1) 1 (http://www.thescpub.com/jmss.toc) A MODIFIED METHOD FOR SOLVING SYSTEM OF

More information

Maejo International Journal of Science and Technology

Maejo International Journal of Science and Technology Maejo Int. J. Sc. Technol. () - Full Paper Maejo Internatonal Journal of Scence and Technology ISSN - Avalable onlne at www.mjst.mju.ac.th Fourth-order method for sngularly perturbed sngular boundary value

More information

Inexact Newton Methods for Inverse Eigenvalue Problems

Inexact Newton Methods for Inverse Eigenvalue Problems Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

Representation theory and quantum mechanics tutorial Representation theory and quantum conservation laws

Representation theory and quantum mechanics tutorial Representation theory and quantum conservation laws Representaton theory and quantum mechancs tutoral Representaton theory and quantum conservaton laws Justn Campbell August 1, 2017 1 Generaltes on representaton theory 1.1 Let G GL m (R) be a real algebrac

More information

Modelli Clamfim Equazioni differenziali 7 ottobre 2013

Modelli Clamfim Equazioni differenziali 7 ottobre 2013 CLAMFIM Bologna Modell 1 @ Clamfm Equazon dfferenzal 7 ottobre 2013 professor Danele Rtell danele.rtell@unbo.t 1/18? Ordnary Dfferental Equatons A dfferental equaton s an equaton that defnes a relatonshp

More information

New Method for Solving Poisson Equation. on Irregular Domains

New Method for Solving Poisson Equation. on Irregular Domains Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad

More information

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl RECURSIVE SPLINE INTERPOLATION METHOD FOR REAL TIME ENGINE CONTROL APPLICATIONS A. Stotsky Volvo Car Corporaton Engne Desgn and Development Dept. 97542, HA1N, SE- 405 31 Gothenburg Sweden. Emal: astotsky@volvocars.com

More information

The Quadratic Trigonometric Bézier Curve with Single Shape Parameter

The Quadratic Trigonometric Bézier Curve with Single Shape Parameter J. Basc. Appl. Sc. Res., (3541-546, 01 01, TextRoad Publcaton ISSN 090-4304 Journal of Basc and Appled Scentfc Research www.textroad.com The Quadratc Trgonometrc Bézer Curve wth Sngle Shape Parameter Uzma

More information

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law: CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and

More information

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016 ME140 - Lnear rcuts - Wnter 16 Fnal, March 16, 2016 Instructons () The exam s open book. You may use your class notes and textbook. You may use a hand calculator wth no communcaton capabltes. () You have

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 21: Numerical methods for pricing American type derivatives Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)

More information

Beyond Zudilin s Conjectured q-analog of Schmidt s problem

Beyond Zudilin s Conjectured q-analog of Schmidt s problem Beyond Zudln s Conectured q-analog of Schmdt s problem Thotsaporn Ae Thanatpanonda thotsaporn@gmalcom Mathematcs Subect Classfcaton: 11B65 33B99 Abstract Usng the methodology of (rgorous expermental mathematcs

More information

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

Homogenization of reaction-diffusion processes in a two-component porous medium with a non-linear flux-condition on the interface

Homogenization of reaction-diffusion processes in a two-component porous medium with a non-linear flux-condition on the interface Homogenzaton of reacton-dffuson processes n a two-component porous medum wth a non-lnear flux-condton on the nterface Internatonal Conference on Numercal and Mathematcal Modelng of Flow and Transport n

More information

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP C O L L O Q U I U M M A T H E M A T I C U M VOL. 80 1999 NO. 1 FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP BY FLORIAN K A I N R A T H (GRAZ) Abstract. Let H be a Krull monod wth nfnte class

More information