A Numerical Patching Technique for Singularly Perturbed Nonlinear Differential-Difference Equations with a Negative Shift

Size: px
Start display at page:

Download "A Numerical Patching Technique for Singularly Perturbed Nonlinear Differential-Difference Equations with a Negative Shift"

Transcription

1 Aled Matematcs, (: - DOI:.593/j.am..4 A Numercal Patcg Tecque for Sgularl Perturbed Nolear Dfferetal-Dfferece Equatos wt a Negatve Sft R. Nageswar Rao, P. Pramod Cakravart Deartmet of Matematcs, Vsvesvaraa Natoal Isttute of Tecolog, Nagur, 44, Ida Abstract I ts aer, we reset a umercal atcg tecque for solvg sgularl erturbed olear dfferetal-dfferece equato wt a small egatve sft. Te olear roblem s coverted to a sequece of lear roblems b quaslearzato rocess. After learzato, t s dvded to two roblems, amel er rego roblem ad outer rego roblem. Te boudar codto at te cuttg ot s obtaed from te teor of sgular erturbatos. Usg stretcg trasformato, a modfed er rego roblem s costructed ad s solved b usg te uwd fte dfferece sceme. Te outer rego roblem s solved b a Talor olomal aroac. We combe te solutos of bot roblems to obta a aromate soluto of te orgal roblem. Te roosed metod s teratve o te cuttg ot. Te rocess s reeated for varous coces of te cuttg ot, utl te soluto rofles stablze. Some umercal eamles ave bee solved to demostrate te alcablt of te metod. Te metod s aalzed for stablt ad covergece. Kewords Sgular Perturbatos, Numercal Patcg Tecque, Nolear Dfferetal-Dfferece Equato, Uwd Fte Dfferece Sceme, Talor Polomal Aroac. Itroducto A sgularl erturbed dfferetal-dfferece equato s a ordar dfferetal equato wc te gest dervatve s multled b a small arameter ad volvg at least oe dela term. Suc roblems are foud trougout te lterature o edemcs ad oulato damcs were tese small sfts la a mortat role te modelg of varous real lfe eomea[]. Boudar value roblems dfferetal dfferece equatos arse a ver atural wa studg varatoal roblems cotrol teor were te roblem s comlcated b te effect of tme delas sgal trasmsso[]. I te matematcal model for te determato of te eected frst-et tme te geerato of acto otetal erve cells b radom satc uts dedrtes, te sfts are due to te jums te otetal membrae wc are ver small[]. I[5] te autors C.G. Lage ad R.M. Mura gave a asmtotc aroac te stud of a class of boudar value roblems for lear secod order dfferetal-dfferece equatos wc te gest order dervatve s multled b a small arameter. I[4] te autors ave eteded te deas to te boudar value roblems for sgularl erturbed olear Corresodg autor: ramododla@aoo.co. (P. Pramod Cakravart Publsed ole at tt://joural.saub.org/f Corgt Scetfc & Academc Publsg. All Rgts Reserved dfferetal dfferece equatos ad dscussed te estece ad uqueess of ter solutos. I[3], te autors roosed a umercal metod to solve te boudar value roblems for sgularl erturbed lear dfferetal-dfferece equatos, wc worked well we te dela arameter s of O( or o(. To adle te dela argumet, te costructed a secal te of mes so tat te term cotag dela les o odal ots after dscretzato. Mustafa Gulsu ad Memet Sezer[6] ave roosed a Talor olomal aroac for solvg mt order lear dfferetal-dfferece equatos wt med codtos. Ts metod s based o frst takg te trucated Talor s easos of te fuctos te dfferetal-dfferece equatos ad te substtutg ter matr forms to te equato. Hece te result matr equato ca be solved ad te ukow Talor coeffcets ca be foud aromatel. I[], te autors M.K. Kadalbajoo ad K.K. Sarma reseted te umercal stud to solve te sgularl erturbed olear dfferetal dfferece equatos wt egatve sft. We te sft s o(, te term cotag te sft s eaded Talor seres ad we te sft s O(, a secal te of mes s costructed to adle te sft term. I[], te autors M.K.Kadalbajoo, Devedra Kumar roosed a B-Sle collocato metod for solvg a sgularl erturbed olear dfferetal dfferece equato wt egatve sft. A ecewse uform mes s used to gras te better aromato to te eact soluto te boudar laer rego. Talor seres s used to tackle te dela.

2 R. Nageswar Rao et al.: A Numercal Patcg Tecque for Sgularl Perturbed Nolear Dfferetal-Dfferece Equatos wt a Negatve Sft I ts aer, we reset a umercal atcg tecque for solvg sgularl erturbed olear dfferetal dfferece equatos wt egatve sft ad reset te umercal stud to suc roblems. I order to kow te beavour of te soluto of te sgularl erturbed dfferetal dfferece equatos te boudar laer rego, t s alwas suggestve to dvde te orgal roblem to two roblems amel, te er rego roblem ad te outer rego roblem ad solve tem searatel. Frst, we learze te olear boudar value roblem usg quaslearzato[8] ad obta a sequece of lear boudar value roblems. Te t s dvded to two roblems amel er rego roblem ad outer rego roblem. Te boudar codto at te cuttg ot s obtaed from te teor of sgular erturbatos. Usg stretcg trasformato, a modfed er rego roblem s costructed ad s solved b usg te uwd fte dfferece sceme. Te outer rego roblem s solved b a Talor olomal aroac gve b Mustafa Gulsu ad Memet Sezer[6]. We combe te solutos of bot te roblems to obta a aromate soluto to te orgal roblem. Te roosed metod s teratve o te cuttg ot. Te rocess s to be reeated for varous coces of te cuttg ot, utl te soluto rofles stablze. Te estece ad uqueess of te dscretzed roblem alog wt stablt estmates are dscussed. Some umercal eamles ave bee solved to demostrate te alcablt of te metod.. Numercal Patcg Tecque To descrbe te metod, we frst cosder a olear sgularl erturbed dfferetal dfferece equato of te form: "( = f (, (, '( δ, ( subject to te terval ad boudar codtos ( = φ(; δ ( = β ( were < << s a small erturbato arameter ad < δ <<. Te soluto ( of te boudar value roblem (-( s assumed to be cotuous o[,] ad cotuousl dfferetable o (,. It s assumed tat f (,, ' s smoot fucto satsfg te codtos: f f f f > M > ' ', were M s a ostve costat. v Te growt codto f (,, ' = O(' as ' for all [,] ad all real ad '. For δ=, uder te codtos lsted above te roblem (-( as a uque soluto[9]. To develo a umercal sceme for te boudar value roblem (-(, we frst learze te orgal olear roblem b usg quaslearzato rocess[8]. Te olear dfferetal equato s learzed aroud a omal soluto of te olear dfferetal equato wc satsfes te secfed boudar codtos. Let ( ( be te tal guess to te soluto of te roblem ( satsfg te boudar codtos ( ( = φ(, δ ; ( ( = β. Te a sequece of boudar value roblems ca be obtaed as follows: We ave (" (k + = f wt ( = φ( (k + (, (,( (k + ( δ (k +, δ ; ( = β Eadg te rgt ad sde of (3 Talor seres about, we obta (k + ( = f (, (,( ( δ (k f (( + ( δ ( ( δ (" ' (k + f ( ( ( +... B rearragg te terms of (4 we obta te recurrece relato of sgularl erturbed lear dfferetal-dfferece equatos (k + f (k + (" ( (' ( δ ' f (k + ( f (, (,(' ( δ = (5 f f ( (' ( δ ' (k + (k + wt ( = φ( o δ, ( = β For smlct, we deote f a ( ' =, f b ( =, F ( = f (, (,(' ( δ f f ( ( ( δ ' Te (5 ma be wrtte as (k + (" a (( (k + + ( δ + (6 (k + b ( ( = F ( wt (k + ( = φ( o δ, (k ( = β Te soluto of te reduced roblem of (-( s assumed to be te tal aromato ( (, ad te successve (3 (4 (7

3 Aled Matematcs, (: - 3 aromatos { (} k = are determed b (6-(7. Hece, stead of solvg te orgal olear roblem (-(, we solve te sequece of boudar value roblems for sgularl erturbed secod order lear dfferetal equatos wt egatve sft gve b (6-(7 for k=,,, For large values of k te solutos ( k ( coverge to te soluto ( of te orgal olear dfferetal equato (k + wle umercall, we requre tat ( ( < μ,, were µ s te rescrbed tolerace. Te terato ca be termated we te above codto s satsfed, ad ( k + te rofle ( s te umercal soluto of te olear boudar value roblem (-(. Due to te resece of te sgular erturbato arameter, te soluto of te roblem ebts boudar laer beavor. As, te order of te reduced roblem decreases b oe, terefore te soluto ebts boudar laer beavor at eter of te boudar ots.e., te boudar laer wll be o te left sde or te rgt sde of te doma of cosderato deedg o te sg of te coeffcet a ( k ( of te covecto term tat s accordg as a ( k ( M > or a ( k ( M < resectvel, were M s a ostve costat. We assume tat a ( k ( M > trougout te terval [, ]. Ts assumto mles tat te boudar laer wll be te egborood of =. We set δ=τ wt τ=o(. If τ s ot too large, te laer structure s modfed but mataed at te same ed [5]. We cosder = O( te cuttg ot or te tckess of te boudar laer. Now te learzed roblem s dvded to two roblems amel te er rego roblem ad te outer rego roblem. Te er rego roblem s defed te terval ad te outer rego roblem s defed te terval... Boudar Codto at te Cuttg Pot We assume tat { (} k = are suffcetl dfferet- able. B eadg te retarded term (' ( δ usg Talor seres, we obta ( ( δ (' ( δ(" ( Sce δ=τ, we ave ( ( δ (' ( τ(" ( (8 B substtutg (8 (6, we get [ a (τ]( ( + a (( ( + (9 b ( ( = F (. We sall seek a outer soluto as a asmtotc easo te form ( k ( ( = = ( were ( k ( are ukow fuctos to be determed. Substtutg ( (9 ad smlfg, we get a a (( wt ( + b ( (( [ a (τ] ( k+ Te soluto of ( s ( = F ( wt ( = β ad ( ( + b ( b (ξ F ( = e dξ a (ξ a ( ( =, =,,3... (s s b e (s a ( k+ ( = ( (ξ dξds + (ξ β Te fuctos (, ( ca be obtaed b solvg equato ( for =,,3. Tus te easo for ( k ( gve equato ( s obtaed. Hece te boudar codto at te cuttg ot ca be obtaed from ( ad deote ( = ( = γ (sa (3 = Sce te termal ot s commo to bot te er ad outer regos, t defes te er rego roblem as a boudar value roblem ( ( + a (( ( δ +, (4 b ( ( = F ( wt te terval ad boudar codtos ( = φ(; δ (5 ( = γ. ad te outer rego roblem as a boudar value roblem ( ( + a (( ( δ +, (6 b ( ( = F ( wt te boudar codtos ( = γ ad ( = β (7.. Ier Rego Problem From (4-(5 te er rego roblem s gve b ( ( + a (( ( δ +, (8 b ( ( = F ( wt te terval ad boudar codtos ( k ( = φ(; δ ( = γ. (9 We coose te trasformato =t to create a ew er rego roblem. B rescalg te equato (8 wt

4 4 R. Nageswar Rao et al.: A Numercal Patcg Tecque for Sgularl Perturbed Nolear Dfferetal-Dfferece Equatos wt a Negatve Sft ( = (t = Y (t (' ( = ( (t = (Y (t (" ( = ( (t = (Y (t a ( = a (t = A (t b ( = b (t = B (t F ( = F (t = Fˆ(t ( we obta te ew dfferetal-dfferece equato for te er rego soluto as (Y" (t + A B (t(y' (ty (t τ + (t = Fˆ(t ( δ were τ =. Te boudar codtos for te equato ( are determed b (9 ad ( as Y Y (t = Φ(t = φ(t; τ t (t = γ. ( We solve te ew er rego roblem (-( to obta te soluto over te terval t t. We costruct a umercal sceme for solvg (-( based o a uwd fte dfferece sceme. We dvde te terval [, t ] to equal arts. To tackle te dela term, we coose te mes τ arameter as =, were m=q, s a ostve teger ad m q s te matssa of τ. Te dfferece sceme for (-( s gve b D+ D Y + A D+ Y m + B Y = Fˆ (3, =,,3...,. ( k Y =Φ, = m, m+,., (4 ad Y ( k + = γ. were Y Y + Y + Y Y D + D Y =, + D Y + = ad A = A (t, B = B (t, Fˆ = Fˆ (t. Now (3-(4 become Y + ( B Y + Y + = Fˆ, for =,, m-. A ( Φ m+ Φ m Y + ( B Y + Y + + A Y =, for = m. Fˆ + A Φ Y + ( B Y + Y + + (5 A ( Y m+ Y m = Fˆ for = m+, m+, - wt Y = Φ = Φ( ad Y ( k + = γ (6 Te dscrete roblem (5-(6 reduces to a sstem of (+ lear dfferece equatos gve b G Y = H, (7 T were Y = [Y Y... Y, ],...,H [ H, H ] T H = ad G = [gj], te ozero etres of te sstem matr beg gve b g(, = g(, = + B for =,,...,m g(, + = H = Fˆ A [ Φ m Φ m] + ad g(, = g(, = + B g(, + = for = m,..., g(, m + = A g(, m = A H = Fˆ H = Φ, H = γ, g(, = = g(, Te sstem (7 s solved b Gauss elmato metod wt artal votg. I fact, a umercal metod or aaltcal metod ca be used..3. Outer Rego Problem Sce te termal ot s commo to bot te er ad outer regos, t defes te outer rego roblem as a boudar value roblem ( ( + a (( ( δ +, (8 b ( ( = F ( wt te boudar codtos ( = γ ad ( = β (9 We solve te outer rego roblem (8-(9 b emlog te Talor olomal aroac gve b Mustafa Gulsu ad Memet Sezer[6] to obta te soluto over te terval..4. Soluto of te Orgal Problem After gettg te soluto of te er rego roblem ad outer rego roblem, we combe bot to obta te aromate soluto of te orgal roblem (-( over te terval. We reeat te rocess for varous coces of, utl te soluto rofles do ot dffer materall from terato to terato. For comutatoal uroses we use a absolute error crtero, amel + ( Y (t ( Y (t σ; t t were ( Y (t s te t terate of te er rego soluto ad σ s te rescrbed tolerace boud.

5 Aled Matematcs, (: Error Estmate Now we sall fd te error estmate for te modfed er rego roblem (3 uder te boudar codtos (4. Trougout te aalss smlct we relace (k Y +, A, B, Fˆ wt Y, A, B, F resectvel. Case : We B (t θ <, were θ s a ostve costat. Lemma : (Dscrete mmum rcle Let π be a mes fucto tat satsfes π ad π ad L (π te π for all =,,,,. Proof: Let k {,,,..., } be suc tat πk = m π ad assume tat π k <. Clearl k {,}. Te for k m we ave We ave L (πk = D D πk B + + (t k πk πk πk + π = k + + B (t k πk = (πk + πk (πk πk + B (t k πk >, sce π B k + πk ad πk πk, (t k < Tus, we ave L ( πk > for k m wc cotradcts te otess tat L ( πk for =,..., of te dscrete mmum rcle. Ts cotradcto arose sce we assumed tat π k <. Terefore our assumto s wrog. Hece, π k. But k {,,,...,} s a arbtrar ostve teger, so π for all =,,...,. Teorem : If A ( t M > ad B ( t θ <, were M ad θ are ostve costats, te te soluto of te dscrete roblem (3 wt te boudar codtos (4 ests, s uque ad satsfes Y θ F + K( Φ + γ, (3, were K s a ostve costat. Here,. s te dscrete l, orm defed b t ma t = Proof: Suose { u } = ad { v } = be two solutos to te dscrete roblem (3, (4. Te, z = u v s a mes fucto satsfg z =, z = for, we ave L (z = L (u L (v. Sce u ad v satsf (5, terefore L (z =,. Tus, te mes fucto z satsfes te otess of te dscrete mmum rcle ad so b te alcato of t to te mes fucto z we get z = u v, (3 Aga let z = (u v, te z satsfes z = = z ad roceedg as above we get L (z =,.. Tus, te dscrete mmum rcle ca be aled for te mes fucto z wc gves z = (u v..e., u v,. (3 Hece, from (3 & (3 we get u v =, wc rove te uqueess of te soluto to te dscrete roblem (3-(4 ad for lear equatos, te estece s mled b uqueess. Now to rove te boud o { Y } = we cosder two barrer fuctos ψ defed b ψ = θ F + K( Φ + γ Y,, were K s a arbtrar ostve costat. Te, we ave ψ = θ F + K( Φ + γ Y, = θ F + ( K Φ Φ + K γ, sce Y = Φ, s ce Φ Φ ad K ad ψ = θ F + K( Φ + γ Y, = θ F + K Φ + K γ γ, sce Y = γ γ γ ad K, sce For m, we ave L (ψ D D ( = + ψ + B (t ψ = B (t θ F + K( Φ + γ L ( (33, Y Also for m, we ave L (Y Y B Y = D+ D + (t (Φ Φ F(t A (t m+ = m Te from (33 we get B L (ψ = (t θ F + K( Φ + γ (Φ Φ F(t A (t m+ m Sce, B (t θ <, tat s B (t θ, we get L B (ψ F + (t K( Φ + γ (Φ Φ F(t A (t m+ m F F(t B + (t K( Φ + γ (34 (Φ Φ A (t m+ m Sce, te above equalt (34 te frst ad secod terms are egatve, so we coose te costat K suc tat te sum of te modul of te frst ad secod terms domates te modulus of te trd term te above equalt. We te obta

6 6 R. Nageswar Rao et al.: A Numercal Patcg Tecque for Sgularl Perturbed Nolear Dfferetal-Dfferece Equatos wt a Negatve Sft L (ψ, m. (35 For m,we ave L (ψ = D+ D ψ + A (td+ ψ m + B (t ψ = B (t θ F K( Φ γ L (Y + + (36 We kow tat for =m, m+,..,-, L (Y D D Y A (t D Y B (t Y = m + = F(t Te from (36 we get L (ψ B (t θ F K( Φ γ = + + F(t Sce, B (t θ <, tat s B (t θ, we get L (ψ F B (t K( Φ γ + + F(t F F(t B + (t K( Φ + γ (37 Hece, L (ψ, m (38 Combg te results of (35 ad (38 we get L (ψ,. Tus a alcato of Lemma to te mes fucto ψ gves ψ θ F K( Φ γ Y = + +, wc roves te requred boud o te dscrete soluto, were θ s a ostve co- { Y } =. Case : We B (t θ > stat. Lemma : (Dscrete mamum rcle Let π be a mes fucto tat satsfes π ad π ad L (π te π for all =,,,,. Proof: Let k {,,,...,} be suc tat πk = ma π ad assume tat π k <. Clearl k {,}. Te, for k m we ave L (π D D π + B (t π k = + k k k πk πk + π = k + + B (t k πk (πk πk (πk πk = + + B (t k πk ad π <, s ce π π k +, B (t π > k k k Tus, we ave L (πk < for k m wc cotradcts te otess tat L (πk for =,,,- of te dscrete mamum rcle. Ts cotradcto arose sce we assumed tat π k <. Terefore our assumto s wrog. Hece, π k. But k {,,,...,} s a arbtrar ostve teger, so π for all =,,,. Teorem : If A (t M > ad B (t θ >, were M ad θ are ostve costats, te te soluto of te dscrete k roblem (3 wt te boudar codtos (4 ests, s uque ad satsfes Y θ F + K( Φ + γ, (39, were K s a ostve costat. Proof: Let { u } = ad { v } = be two solutos to te dscrete roblem (3, (4. Te, z = u v s a mes fucto satsfg z =, z, we ave Sce = L (z = L (u L (v, u ad v satsf (3, terefore L (z =,. Tus, te mes fucto z satsfes te otess of te dscrete mamum rcle ad so b te alcato of t to te mes fucto z we get z = u v, (4 Aga we cosder z = (u v, te z s a mes fucto satsfg z = = z ad roceedg as above we get L (z =,. Tus, te dscrete mamum rcle s aled o te mes fucto z wc gves z = (u v..e., u v,. (4 Hece, from te (4, (4 we get u v =, for wc roves te uqueess of te soluto to te dscrete roblem (3-(4, ad for lear equatos te estece s mled b uqueess. Now to rove te boud o { Y } = we cosder two barrer fuctos ψ defed b ψ = θ F + K Φ + γ Y,, ( were K s a arbtrar ostve costat. Te, we ave ψ = θ F + K Φ + γ Y ( = θ F ( K Φ Φ + + K γ, sce Y = Φ, sce Φ Φ ad K ad ψ = θ F + K( Φ + γ Y, = θ F + K Φ + ( K γ γ, sce Y γ =, sce γ γ ad K. ad for m, we ave L (ψ = D+ D ψ + B (t ψ = B ( t θ F K( Φ γ L (Y + + (4 We kow tat for =,,..m-,

7 Aled Matematcs, (: - 7 L (Y Y B Y = D+ D + (t = (Φ Φ F(t A (t m+ m Te from (4 we get B L (ψ = (t θ F + K( Φ + γ F(t (Φ Φ A (t m+ m Sce, B (t θ >, tat s B (t θ, we get L (ψ F + B t K Φ + γ F F(t ( ( (Φ A (t m+ Φ F(t B + (tk (Φ Φ A m+ (t m ( Φ + γ m (43 Sce, te above equalt (43 te frst ad secod terms are ostve, so we coose te costat K so tat te sum of te modul of te frst ad secod terms domates te modulus of te trd term te above equalt. We te obta For m, we ave L (ψ L (ψ, m. (44 = D+ D ψ + A (td + ψ m + B (t ψ + + = B (t θ F K( Φ γ L We kow tat for =m, m+,.., -, L (Y D D Y + A (t D Y + B (t Y = + + m = (Y F(t Te from (45 we get L (ψ B (t θ F K Φ γ = + + F(t Sce, B (t θ >, tat s B (t θ, we get L (ψ F + B (t K Φ + γ F(t ( ( K( Φ + γ (45 F F(t B + (t (46 Hece, L (ψ, m (47 Combg te results of (44 ad (47 we get L (ψ,. Tus a alcato of Lemma to te mes fucto gves ψ = θ F + K Φ + γ Y, ( ψ wc roves te requred boud o te dscrete soluto { Y } =. Tus, Teorems ( ad ( ml tat te soluto to te dscrete roblem (3, (4 s uforml bouded, deedetl of te mes arameter ad te arameter, wc roves tat te dfferece sceme s stable for all mes szes. Corollar: Uder te assumto tat A (t M >, te error e = Y (t Y betwee te soluto Y (t of te cotuous roblem (, ( ad te soluto Y of te dscrete roblem (3,(4 satsfes te estmate e θ T, were T,, satsfes ( Y 4 ma (t + 6 t t t+ T ( A Y + Y (t ma "(t (t + Y 4 (t t m t t m+ 3 Proof: Te trucato error T s gve b Y Y + Y + T = Y (t + Y + Y A m m (t Y (t m Now usg Talor s seres ad after smlfcatos, we obta Y (4 ma (t + 6 t t t + T A Y + Y (t ma "(t (t + Y (4 (t t m t t m+ 3 We ave L ( ( Y L ( Y e(t = L (t = T, =,,,- ad e =e =. Te b usg Teorems ( ad ( we obta te requred error estmate. 4. Numercal Results We ave aled te reset metod o two olear sgularl erturbed dfferetal-dfferece equatos wt small egatve sft. Te olear roblems are frst coverted to sequece lear sgularl erturbed dfferetal-dfferece equatos b usg quaslearzato metod. Te soluto of te reduced roblem s take as tal aromato. ( Eamle [,.593]: "( + ' ( δ e =, wt te terval ad boudar codtos ( = ; δ, ( = Te soluto at te cuttg ot s gve b te reduced roblem soluto ( = ( + ( were ( = loge 3 τ 3 ad ( = loge (3 ad ece = τ ( loge 3 3

8 8 R. Nageswar Rao et al.: A Numercal Patcg Tecque for Sgularl Perturbed Nolear Dfferetal-Dfferece Equatos wt a Negatve Sft Te learzed form of te gve olear equato b quaslearzato metod s "( + ' ( δ ( = log e wt te codtos ( = ; δ ad ( =. Te modfed er rego roblem s gve b Y"(t + Y' (t τ Y(t = log 3 t, t t e 3 t 3 t uder te codtos Y ( =, Y(t = γ. Te outer rego roblem s gve b "( + '( δ ( = loge 3 3 3, uder te codtos ( = γ, ( = Te umercal results are gve tables & for τ=.5 ad = 3, = 4 resectvel. For τ=.5 adτ=.5, te er soluto s lotted gras ad sow fg. to fg 4 for = 3, = 4 resectvel. Numercal Soluto Table. Numercal results of Eamle for = 3, τ=.5. (t=3 (t=6 (t= Fgure. Ier soluto of Eamle for = 3, τ=.5. Numercal Soluto Numercal Soluto Numercal Soluto Fgure. Ier soluto of Eamle for = 3, τ= Fgure 3. Ier soluto of Eamle for = 4, τ= Fgure 4. Ier soluto of Eamle for = 4, τ=.5. Table. Numercal results of Eamle for = 4, τ=.5. (t=3 (t=6 (t= Eamle [,.e9]: "( + (.'( δ ( = wt te terval ad boudar codtos

9 Aled Matematcs, (: - 9 ( = ; δ, ( = Te soluto at te cuttg ot s gve b te reduced roblem soluto ( = ( + ( were ( = ad ( = ad ece ( = Te learzed form of te gve olear equato b quaslearzato metod s "( + '( δ = wt te codtos ( = ; δ, ( = Te er rego roblem s gve b Y "(t + ty'(t τ = t, t t uder te codtos Y ( =, Y(t = γ Te outer rego roblem s gve b "( + '( δ =, uder te codtos ( = γ, ( = Te umercal results are gve tables 3 & 4 for τ=.5 ad = 3, = 4 resectvel. For τ=.5 ad τ=.5, te umercal soluto s lotted gras ad sow fg. 5 to fg. 8 for = 3, = 4 resectvel. Table 3. Numercal results of Eamle for = -3, τ=.5. Numercal Soluto Numercal Soluto Fgure 6. Numercal soluto of Eamle for = 3, τ= Fgure 7. Numercal soluto of Eamle for = 4, τ=.5..8 Numercal Soluto (t=3 (t=6 (t= Fgure 5. Numercal soluto of Eamle for = 3, τ=.5. Numercal Soluto Fgure 8. Numercal soluto of Eamle for = 4, τ=.5. Table 4. Numercal results of Eamle for = -4, τ=.5. (t=3 (t=6 (t=

10 R. Nageswar Rao et al.: A Numercal Patcg Tecque for Sgularl Perturbed Nolear Dfferetal-Dfferece Equatos wt a Negatve Sft 5. Coclusos I order to kow te beavor of te soluto of te sgularl erturbed dfferetal-dfferece equatos te boudar laer rego, t s alwas suggestve to dvde te orgal roblem to two roblems amel te er rego roblem ad te outer rego roblem ad solve tem searatel. We ave reseted a umercal atcg tecque for solvg sgularl erturbed olear dfferetal-dfferece equatos wt te boudar laer at oe ed ot. Te orgal olear boudar value roblem s learzed usg quaslearzato. Te t s dvded to two roblems amel er rego roblem ad outer rego roblem. Te boudar codto at te cuttg ot s obtaed from te teor of sgular erturbatos. A ew er rego roblem s costructed ad solved b usg uwd fte dfferece sceme. Te outer rego roblem s solved b Talor olomal aroac. We ave mlemeted te reset metod o two olear eamles ebtg a left ed boudar laer. Te roosed metod s teratve o te cuttg ot. Te rocess s to be reeated for varous coces of te cuttg ot, utl te soluto rofles do ot dffer materall from terato to terato. Numercal results are reseted tables. It ca be observed from te tables tat te reset metod aromates te soluto avalable te lterature as well. ACKNOWLEDGMENTS Te autors ws to tak te Deartmet of Scece & Tecolog, Govermet of Ida, for ter facal suort uder te roject No. SR/S4/MS: 598/9. (5 e99-e94 [] Moa K. Kadalbajoo ad Devedra Kumar, A comutatoal metod for sgularl erturbed olear dfferetal-dfferece equatos wt small sft, Aled Matematcal Modellg 34 ( [3] Moa K. Kadalbajoo ad Kal K. Sarma, A umercal metod based o fte dfferece for boudar value roblems for sgularl erturbed dela dfferetal equatos, Aled Matematcs ad Comutato 97 ( [4] C.G. Lage ad R.M. Mura, Sgular erturbato aalss of boudar-value roblems for dfferetal dfferece equatos, IV, a olear eamle wt laer beavor, Stud. Al. Mat. 84 ( [5] C.G. Lage ad R.M. Mura, Sgular erturbato aalss of boudar-value roblems for dfferetal dfferece equatos. V. Small sfts wt laer beavor, SIAM Joural o Aled Matematcs 54 ( [6] Mustafa Gulsu ad Memet Sezer, A Talor olomal aroac for solvg dfferetal- dfferece equatos, Joural of Comutatoal ad Aled Matematcs 86 ( [7] E.P. Doola, J.J.H. Mller, W.H.A. Sclders, Uform Numercal Metods for Problems wt Ital ad Boudar Laers, Boole Press, Dubl, 98 [8] R.E. Bellma, R.E. Kalaba, Quaslearzato ad Nolear Boudar Value Problems, Amerca Elsever, New York, 965 [9] F.A. Howes, Sgular Perturbatos ad Dfferetal Iequaltes, vol. 68, Memors of te Amerca Matematcal Socet, Provdece, RI, 976 [] Y. Kuag, Dela Dfferetal equatos wt alcatos oulato damcs, Academc Press, 993 [] R.B. Ste, A teoretcal aalss of euroal varablt, Boscal Joural 5 ( REFERENCES [] Moa K. Kadalbajoo ad Kal K. Sarma, Numercal treatmet for sgularl erturbed olear dfferetal dfferece equatos wt egatve sft, Nolear Aalss 63 [] L.E. El sgol ts, Qualtatve Metods Matematcal Aalses, Traslatos of Matematcal Moogras, Amerca matematcal socet, Provdece, RI, 964

Numerical Analysis Formulae Booklet

Numerical 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 information

On the Behavior of Positive Solutions of a Difference. equation system:

On the Behavior of Positive Solutions of a Difference. equation system: Aled Mathematcs -8 htt://d.do.org/.6/am..9a Publshed Ole Setember (htt://www.scr.org/joural/am) O the Behavor of Postve Solutos of a Dfferece Equatos Sstem * Decu Zhag Weqag J # Lg Wag Xaobao L Isttute

More information

BERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS. Aysegul Akyuz Dascioglu and Nese Isler

BERNSTEIN 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 information

Unique Common Fixed Point of Sequences of Mappings in G-Metric Space M. Akram *, Nosheen

Unique Common Fixed Point of Sequences of Mappings in G-Metric Space M. Akram *, Nosheen Vol No : Joural of Facult of Egeerg & echolog JFE Pages 9- Uque Coo Fed Pot of Sequeces of Mags -Metrc Sace M. Ara * Noshee * Deartet of Matheatcs C Uverst Lahore Pasta. Eal: ara7@ahoo.co Deartet of Matheatcs

More information

UNIT 2 SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS

UNIT 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 information

2006 Jamie Trahan, Autar Kaw, Kevin Martin University of South Florida United States of America

2006 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 information

Asymptotic Formulas Composite Numbers II

Asymptotic Formulas Composite Numbers II Iteratoal Matematcal Forum, Vol. 8, 203, o. 34, 65-662 HIKARI Ltd, www.m-kar.com ttp://d.do.org/0.2988/mf.203.3854 Asymptotc Formulas Composte Numbers II Rafael Jakmczuk Dvsó Matemátca, Uversdad Nacoal

More information

An Iterative Solution for Second Order Linear Fredholm Integro-Differential Equations

An Iterative Solution for Second Order Linear Fredholm Integro-Differential Equations Malasa Joural of Matematcal Sceces 8): 57-7 4) MLYSIN JOURNL OF MTHEMTICL SCIENCES Joural omeage: tt://esem.um.edu.m/oural Iteratve Soluto for Secod Order Lear Fredolm Itegro-Dfferetal Equatos * Elaaraa

More information

Multiple Choice Test. Chapter Adequacy of Models for Regression

Multiple 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

Nonlinear Piecewise-Defined Difference Equations with Reciprocal Quadratic Terms

Nonlinear Piecewise-Defined Difference Equations with Reciprocal Quadratic Terms Joural of Matematcs ad Statstcs Orgal Researc Paper Nolear Pecewse-Defed Dfferece Equatos wt Recprocal Quadratc Terms Ramada Sabra ad Saleem Safq Al-Asab Departmet of Matematcs, Faculty of Scece, Jaza

More information

Fourth Order Four-Stage Diagonally Implicit Runge-Kutta Method for Linear Ordinary Differential Equations ABSTRACT INTRODUCTION

Fourth 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 information

Numerical Solutions of Second Order Initial Value Problems of Bratu-Type equation Using Higher Ordered Rungu-Kutta Method

Numerical 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 information

Beam Warming Second-Order Upwind Method

Beam 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 information

FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM

FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM Joural of Appled Matematcs ad Computatoal Mecacs 04, 3(4), 7-34 FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM Ata Cekot, Stasław Kukla Isttute of Matematcs, Czestocowa Uversty of Tecology Częstocowa,

More information

1 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

1 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 information

Toni Alchofino ( ) Physics Department, Bandung Institute of Technology, Indonesia,

Toni Alchofino ( ) Physics Department, Bandung Institute of Technology, Indonesia, Study of Numercal alyss for Coductve System Usg Fte Dfferece lgortms for Smle Oe Dmesoal Heat-Coducto Problems Wt Costat Termal Dffusvtes To lcofo (53 Pyscs Deartmet, Badug Isttute of Tecology, Idoesa,

More information

Mu Sequences/Series Solutions National Convention 2014

Mu 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 information

( ) 2 2. Multi-Layer Refraction Problem Rafael Espericueta, Bakersfield College, November, 2006

( ) 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 information

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

UNIVERSITY 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 information

Initial-Value Problems for ODEs. numerical errors (round-off and truncation errors) Consider a perturbed system: dz dt

Initial-Value Problems for ODEs. numerical errors (round-off and truncation errors) Consider a perturbed system: dz dt Ital-Value Problems or ODEs d GIVEN: t t,, a FIND: t or atb umercal errors (roud-o ad trucato errors) Cosder a perturbed sstem: dz t, z t, at b z a a Does z(t) (t)? () (uqueess) a uque soluto (t) exsts

More information

Uniform asymptotical stability of almost periodic solution of a discrete multispecies Lotka-Volterra competition system

Uniform asymptotical stability of almost periodic solution of a discrete multispecies Lotka-Volterra competition system Iteratoal Joural of Egeerg ad Advaced Research Techology (IJEART) ISSN: 2454-9290, Volume-2, Issue-1, Jauary 2016 Uform asymptotcal stablty of almost perodc soluto of a dscrete multspeces Lotka-Volterra

More information

General Method for Calculating Chemical Equilibrium Composition

General Method for Calculating Chemical Equilibrium Composition AE 6766/Setzma Sprg 004 Geeral Metod for Calculatg Cemcal Equlbrum Composto For gve tal codtos (e.g., for gve reactats, coose te speces to be cluded te products. As a example, for combusto of ydroge wt

More information

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

PGE 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 information

STRONG CONSISTENCY FOR SIMPLE LINEAR EV MODEL WITH v/ -MIXING

STRONG CONSISTENCY FOR SIMPLE LINEAR EV MODEL WITH v/ -MIXING Joural of tatstcs: Advaces Theory ad Alcatos Volume 5, Number, 6, Pages 3- Avalable at htt://scetfcadvaces.co. DOI: htt://d.do.org/.864/jsata_7678 TRONG CONITENCY FOR IMPLE LINEAR EV MODEL WITH v/ -MIXING

More information

On the characteristics of partial differential equations

On the characteristics of partial differential equations Sur les caractérstques des équatos au dérvées artelles Bull Soc Math Frace 5 (897) 8- O the characterstcs of artal dfferetal equatos By JULES BEUDON Traslated by D H Delhech I a ote that was reseted to

More information

ON BIVARIATE GEOMETRIC DISTRIBUTION. K. Jayakumar, D.A. Mundassery 1. INTRODUCTION

ON BIVARIATE GEOMETRIC DISTRIBUTION. K. Jayakumar, D.A. Mundassery 1. INTRODUCTION STATISTICA, ao LXVII, 4, 007 O BIVARIATE GEOMETRIC DISTRIBUTIO ITRODUCTIO Probablty dstrbutos of radom sums of deedetly ad detcally dstrbuted radom varables are maly aled modelg ractcal roblems that deal

More information

Australian Journal of Basic and Applied Sciences. Full-Sweep SOR Iterative Method to Solve Space-Fractional Diffusion Equations

Australian Journal of Basic and Applied Sciences. Full-Sweep SOR Iterative Method to Solve Space-Fractional Diffusion Equations Australa Joural of Basc ad Aled Sceces, 8(4) Secal 14, Paes: 153-158 AENSI Jourals Australa Joural of Basc ad Aled Sceces ISSN:1991-8178 Joural home ae: www.abasweb.com Full-Swee SOR Iteratve Method to

More information

Entropy, Relative Entropy and Mutual Information

Entropy, Relative Entropy and Mutual Information Etro Relatve Etro ad Mutual Iformato rof. Ja-Lg Wu Deartmet of Comuter Scece ad Iformato Egeerg Natoal Tawa Uverst Defto: The Etro of a dscrete radom varable s defed b : base : 0 0 0 as bts 0 : addg terms

More information

STK3100 and STK4100 Autumn 2017

STK3100 and STK4100 Autumn 2017 SK3 ad SK4 Autum 7 Geeralzed lear models Part III Covers the followg materal from chaters 4 ad 5: Sectos 4..5, 4.3.5, 4.3.6, 4.4., 4.4., ad 4.4.3 Sectos 5.., 5.., ad 5.5. Ørulf Borga Deartmet of Mathematcs

More information

A Conventional Approach for the Solution of the Fifth Order Boundary Value Problems Using Sixth Degree Spline Functions

A 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 information

STK3100 and STK4100 Autumn 2018

STK3100 and STK4100 Autumn 2018 SK3 ad SK4 Autum 8 Geeralzed lear models Part III Covers the followg materal from chaters 4 ad 5: Cofdece tervals by vertg tests Cosder a model wth a sgle arameter β We may obta a ( α% cofdece terval for

More information

A stopping criterion for Richardson s extrapolation scheme. under finite digit arithmetic.

A stopping criterion for Richardson s extrapolation scheme. under finite digit arithmetic. A stoppg crtero for cardso s extrapoato sceme uder fte dgt artmetc MAKOO MUOFUSHI ad HIEKO NAGASAKA epartmet of Lbera Arts ad Sceces Poytecc Uversty 4-1-1 Hasmotoda,Sagamara,Kaagawa 229-1196 JAPAN Abstract:

More information

Supervised learning: Linear regression Logistic regression

Supervised learning: Linear regression Logistic regression CS 57 Itroducto to AI Lecture 4 Supervsed learg: Lear regresso Logstc regresso Mlos Hauskrecht mlos@cs.ptt.edu 539 Seott Square CS 57 Itro to AI Data: D { D D.. D D Supervsed learg d a set of eamples s

More information

0/1 INTEGER PROGRAMMING AND SEMIDEFINTE PROGRAMMING

0/1 INTEGER PROGRAMMING AND SEMIDEFINTE PROGRAMMING CONVEX OPIMIZAION AND INERIOR POIN MEHODS FINAL PROJEC / INEGER PROGRAMMING AND SEMIDEFINE PROGRAMMING b Luca Buch ad Natala Vktorova CONENS:.Itroducto.Formulato.Applcato to Kapsack Problem 4.Cuttg Plaes

More information

ELEC 6041 LECTURE NOTES WEEK 1 Dr. Amir G. Aghdam Concordia University

ELEC 6041 LECTURE NOTES WEEK 1 Dr. Amir G. Aghdam Concordia University ELEC 604 LECTURE NOTES WEEK Dr mr G ghdam Cocorda Uverst Itroducto - Large-scale sstems are the mult-ut mult-outut (MIMO) sstems cosstg of geograhcall searated comoets - Eamles of large-scale sstems clude

More information

Bivariate Vieta-Fibonacci and Bivariate Vieta-Lucas Polynomials

Bivariate Vieta-Fibonacci and Bivariate Vieta-Lucas Polynomials IOSR Joural of Mathematcs (IOSR-JM) e-issn: 78-78, p-issn: 19-76X. Volume 1, Issue Ver. II (Jul. - Aug.016), PP -0 www.osrjourals.org Bvarate Veta-Fboacc ad Bvarate Veta-Lucas Polomals E. Gokce KOCER 1

More information

MEASURES OF DISPERSION

MEASURES OF DISPERSION MEASURES OF DISPERSION Measure of Cetral Tedecy: Measures of Cetral Tedecy ad Dsperso ) Mathematcal Average: a) Arthmetc mea (A.M.) b) Geometrc mea (G.M.) c) Harmoc mea (H.M.) ) Averages of Posto: a) Meda

More information

Analysis of Lagrange Interpolation Formula

Analysis 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 information

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

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 information

Strong Convergence of Weighted Averaged Approximants of Asymptotically Nonexpansive Mappings in Banach Spaces without Uniform Convexity

Strong 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 information

Comparing Different Estimators of three Parameters for Transmuted Weibull Distribution

Comparing Different Estimators of three Parameters for Transmuted Weibull Distribution Global Joural of Pure ad Appled Mathematcs. ISSN 0973-768 Volume 3, Number 9 (207), pp. 55-528 Research Ida Publcatos http://www.rpublcato.com Comparg Dfferet Estmators of three Parameters for Trasmuted

More information

Study of Correlation using Bayes Approach under bivariate Distributions

Study of Correlation using Bayes Approach under bivariate Distributions Iteratoal Joural of Scece Egeerg ad Techolog Research IJSETR Volume Issue Februar 4 Stud of Correlato usg Baes Approach uder bvarate Dstrbutos N.S.Padharkar* ad. M.N.Deshpade** *Govt.Vdarbha Isttute of

More information

Approximate Solution of the Singular-Perturbation Problem on Chebyshev-Gauss Grid

Approximate Solution of the Singular-Perturbation Problem on Chebyshev-Gauss Grid Amerca Joural of Computatoal Mathematcs,,, 9-8 do:.436/ajcm..44 Publshed Ole December (http://www.scrp.org/joural/ajcm) Approxmate Soluto of the Sgular-Perturbato Problem o Chebyshev-Gauss Grd Abstract

More information

Application of Generating Functions to the Theory of Success Runs

Application of Generating Functions to the Theory of Success Runs Aled Mathematcal Sceces, Vol. 10, 2016, o. 50, 2491-2495 HIKARI Ltd, www.m-hkar.com htt://dx.do.org/10.12988/ams.2016.66197 Alcato of Geeratg Fuctos to the Theory of Success Rus B.M. Bekker, O.A. Ivaov

More information

Lecture 12 APPROXIMATION OF FIRST ORDER DERIVATIVES

Lecture 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 information

2SLS Estimates ECON In this case, begin with the assumption that E[ i

2SLS Estimates ECON In this case, begin with the assumption that E[ i SLS Estmates ECON 3033 Bll Evas Fall 05 Two-Stage Least Squares (SLS Cosder a stadard lear bvarate regresso model y 0 x. I ths case, beg wth the assumto that E[ x] 0 whch meas that OLS estmates of wll

More information

MMJ 1113 FINITE ELEMENT METHOD Introduction to PART I

MMJ 1113 FINITE ELEMENT METHOD Introduction to PART I MMJ FINITE EEMENT METHOD Cotut requremets Assume that the fuctos appearg uder the tegral the elemet equatos cota up to (r) th order To esure covergece N must satsf Compatblt requremet the fuctos must have

More information

Idea is to sample from a different distribution that picks points in important regions of the sample space. Want ( ) ( ) ( ) E f X = f x g x dx

Idea is to sample from a different distribution that picks points in important regions of the sample space. Want ( ) ( ) ( ) E f X = f x g x dx Importace Samplg Used for a umber of purposes: Varace reducto Allows for dffcult dstrbutos to be sampled from. Sestvty aalyss Reusg samples to reduce computatoal burde. Idea s to sample from a dfferet

More information

Numerical Simulations of the Complex Modied Korteweg-de Vries Equation. Thiab R. Taha. The University of Georgia. Abstract

Numerical Simulations of the Complex Modied Korteweg-de Vries Equation. Thiab R. Taha. The University of Georgia. Abstract Numercal Smulatos of the Complex Moded Korteweg-de Vres Equato Thab R. Taha Computer Scece Departmet The Uversty of Georga Athes, GA 002 USA Tel 0-542-2911 e-mal thab@cs.uga.edu Abstract I ths paper mplemetatos

More information

Test Paper-II. 1. If sin θ + cos θ = m and sec θ + cosec θ = n, then (a) 2n = m (n 2 1) (b) 2m = n (m 2 1) (c) 2n = m (m 2 1) (d) none of these

Test Paper-II. 1. If sin θ + cos θ = m and sec θ + cosec θ = n, then (a) 2n = m (n 2 1) (b) 2m = n (m 2 1) (c) 2n = m (m 2 1) (d) none of these Test Paer-II. If s θ + cos θ = m ad sec θ + cosec θ =, the = m ( ) m = (m ) = m (m ). If a ABC, cos A = s B, the t s C a osceles tragle a eulateral tragle a rght agled tragle. If cos B = cos ( A+ C), the

More information

( ) = ( ) ( ) Chapter 13 Asymptotic Theory and Stochastic Regressors. Stochastic regressors model

( ) = ( ) ( ) Chapter 13 Asymptotic Theory and Stochastic Regressors. Stochastic regressors model Chapter 3 Asmptotc Theor ad Stochastc Regressors The ature of eplaator varable s assumed to be o-stochastc or fed repeated samples a regresso aalss Such a assumpto s approprate for those epermets whch

More information

ANALYSIS ON THE NATURE OF THE BASIC EQUATIONS IN SYNERGETIC INTER-REPRESENTATION NETWORK

ANALYSIS 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 information

1 Lyapunov Stability Theory

1 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 information

Research Article A New Iterative Method for Common Fixed Points of a Finite Family of Nonexpansive Mappings

Research 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 information

Correlation and Regression Analysis

Correlation and Regression Analysis Chapter V Correlato ad Regresso Aalss R. 5.. So far we have cosdered ol uvarate dstrbutos. Ma a tme, however, we come across problems whch volve two or more varables. Ths wll be the subject matter of the

More information

Regression and the LMS Algorithm

Regression and the LMS Algorithm CSE 556: Itroducto to Neural Netorks Regresso ad the LMS Algorthm CSE 556: Regresso 1 Problem statemet CSE 556: Regresso Lear regresso th oe varable Gve a set of N pars of data {, d }, appromate d b a

More information

On A Two Dimensional Finsler Space Whose Geodesics Are Semi- Elipses and Pair of Straight Lines

On A Two Dimensional Finsler Space Whose Geodesics Are Semi- Elipses and Pair of Straight Lines IOSR Joural of Mathematcs (IOSR-JM) e-issn: 78-578 -ISSN:39-765X Volume 0 Issue Ver VII (Mar-Ar 04) PP 43-5 wwwosrjouralsorg O A Two Dmesoal Fsler Sace Whose Geodescs Are Sem- Elses ad Par of Straght es

More information

å 1 13 Practice Final Examination Solutions - = CS109 Dec 5, 2018

å 1 13 Practice Final Examination Solutions - = CS109 Dec 5, 2018 Chrs Pech Fal Practce CS09 Dec 5, 08 Practce Fal Examato Solutos. Aswer: 4/5 8/7. There are multle ways to obta ths aswer; here are two: The frst commo method s to sum over all ossbltes for the rak of

More information

2. Independence and Bernoulli Trials

2. Independence and Bernoulli Trials . Ideedece ad Beroull Trals Ideedece: Evets ad B are deedet f B B. - It s easy to show that, B deedet mles, B;, B are all deedet ars. For examle, ad so that B or B B B B B φ,.e., ad B are deedet evets.,

More information

CS 1675 Introduction to Machine Learning Lecture 12 Support vector machines

CS 1675 Introduction to Machine Learning Lecture 12 Support vector machines CS 675 Itroducto to Mache Learg Lecture Support vector maches Mlos Hauskrecht mlos@cs.ptt.edu 539 Seott Square Mdterm eam October 9, 7 I-class eam Closed book Stud materal: Lecture otes Correspodg chapters

More information

DKA method for single variable holomorphic functions

DKA 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 information

Evaluating Polynomials

Evaluating Polynomials Uverst of Nebraska - Lcol DgtalCommos@Uverst of Nebraska - Lcol MAT Exam Expostor Papers Math the Mddle Isttute Partershp 7-7 Evaluatg Polomals Thomas J. Harrgto Uverst of Nebraska-Lcol Follow ths ad addtoal

More information

9.1 Introduction to the probit and logit models

9.1 Introduction to the probit and logit models EC3000 Ecoometrcs Lecture 9 Probt & Logt Aalss 9. Itroducto to the probt ad logt models 9. The logt model 9.3 The probt model Appedx 9. Itroducto to the probt ad logt models These models are used regressos

More information

d 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

d 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 information

= lim. (x 1 x 2... x n ) 1 n. = log. x i. = M, n

= lim. (x 1 x 2... x n ) 1 n. = log. x i. = M, n .. Soluto of Problem. M s obvously cotuous o ], [ ad ], [. Observe that M x,..., x ) M x,..., x ) )..) We ext show that M s odecreasg o ], [. Of course.) mles that M s odecreasg o ], [ as well. To show

More information

Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations

Derivation 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 information

Lecture 9: Tolerant Testing

Lecture 9: Tolerant Testing Lecture 9: Tolerat Testg Dael Kae Scrbe: Sakeerth Rao Aprl 4, 07 Abstract I ths lecture we prove a quas lear lower boud o the umber of samples eeded to do tolerat testg for L dstace. Tolerat Testg We have

More information

Comparison of Dual to Ratio-Cum-Product Estimators of Population Mean

Comparison of Dual to Ratio-Cum-Product Estimators of Population Mean Research Joural of Mathematcal ad Statstcal Sceces ISS 30 6047 Vol. 1(), 5-1, ovember (013) Res. J. Mathematcal ad Statstcal Sc. Comparso of Dual to Rato-Cum-Product Estmators of Populato Mea Abstract

More information

Introducing Sieve of Eratosthenes as a Theorem

Introducing Sieve of Eratosthenes as a Theorem ISSN(Ole 9-8 ISSN (Prt - Iteratoal Joural of Iovatve Research Scece Egeerg ad echolog (A Hgh Imact Factor & UGC Aroved Joural Webste wwwrsetcom Vol Issue 9 Setember Itroducg Seve of Eratosthees as a heorem

More information

Binary classification: Support Vector Machines

Binary classification: Support Vector Machines CS 57 Itroducto to AI Lecture 6 Bar classfcato: Support Vector Maches Mlos Hauskrecht mlos@cs.ptt.edu 539 Seott Square CS 57 Itro to AI Supervsed learg Data: D { D, D,.., D} a set of eamples D, (,,,,,

More information

CS 2750 Machine Learning Lecture 5. Density estimation. Density estimation

CS 2750 Machine Learning Lecture 5. Density estimation. Density estimation CS 750 Mache Learg Lecture 5 esty estmato Mlos Hausrecht mlos@tt.edu 539 Seott Square esty estmato esty estmato: s a usuervsed learg roblem Goal: Lear a model that rereset the relatos amog attrbutes the

More information

Can we take the Mysticism Out of the Pearson Coefficient of Linear Correlation?

Can we take the Mysticism Out of the Pearson Coefficient of Linear Correlation? Ca we tae the Mstcsm Out of the Pearso Coeffcet of Lear Correlato? Itroducto As the ttle of ths tutoral dcates, our purpose s to egeder a clear uderstadg of the Pearso coeffcet of lear correlato studets

More information

Solving Constrained Flow-Shop Scheduling. Problems with Three Machines

Solving 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 information

Special Instructions / Useful Data

Special Instructions / Useful Data JAM 6 Set of all real umbers P A..d. B, p Posso Specal Istructos / Useful Data x,, :,,, x x Probablty of a evet A Idepedetly ad detcally dstrbuted Bomal dstrbuto wth parameters ad p Posso dstrbuto wth

More information

Chapter 3. Differentiation 3.2 Differentiation Rules for Polynomials, Exponentials, Products and Quotients

Chapter 3. Differentiation 3.2 Differentiation Rules for Polynomials, Exponentials, Products and Quotients 3.2 Dfferetato Rules 1 Capter 3. Dfferetato 3.2 Dfferetato Rules for Polyomals, Expoetals, Proucts a Quotets Rule 1. Dervatve of a Costat Fucto. If f as te costat value f(x) = c, te f x = [c] = 0. x Proof.

More information

Point Estimation: definition of estimators

Point Estimation: definition of estimators Pot Estmato: defto of estmators Pot estmator: ay fucto W (X,..., X ) of a data sample. The exercse of pot estmato s to use partcular fuctos of the data order to estmate certa ukow populato parameters.

More information

Econometric Methods. Review of Estimation

Econometric Methods. Review of Estimation Ecoometrc Methods Revew of Estmato Estmatg the populato mea Radom samplg Pot ad terval estmators Lear estmators Ubased estmators Lear Ubased Estmators (LUEs) Effcecy (mmum varace) ad Best Lear Ubased Estmators

More information

Algorithms behind the Correlation Setting Window

Algorithms behind the Correlation Setting Window Algorths behd the Correlato Settg Wdow Itroducto I ths report detaled forato about the correlato settg pop up wdow s gve. See Fgure. Ths wdow s obtaed b clckg o the rado butto labelled Kow dep the a scree

More information

ECE606: Solid State Devices Lecture 13 Solutions of the Continuity Eqs. Analytical & Numerical

ECE606: Solid State Devices Lecture 13 Solutions of the Continuity Eqs. Analytical & Numerical ECE66: Sold State Devces Lecture 13 Solutos of the Cotuty Eqs. Aalytcal & Numercal Gerhard Klmeck gekco@purdue.edu Outle Aalytcal Solutos to the Cotuty Equatos 1) Example problems ) Summary Numercal Solutos

More information

Introduction to Numerical Differentiation and Interpolation March 10, !=1 1!=1 2!=2 3!=6 4!=24 5!= 120

Introduction to Numerical Differentiation and Interpolation March 10, !=1 1!=1 2!=2 3!=6 4!=24 5!= 120 Itroducto to Numercal Deretato ad Iterpolato Marc, Itroducto to Numercal Deretato ad Iterpolato Larr Caretto Mecacal Egeerg 9 Numercal Aalss o Egeerg stems Marc, Itroducto Iterpolato s te use o a dscrete

More information

Random Variables. ECE 313 Probability with Engineering Applications Lecture 8 Professor Ravi K. Iyer University of Illinois

Random Variables. ECE 313 Probability with Engineering Applications Lecture 8 Professor Ravi K. Iyer University of Illinois Radom Varables ECE 313 Probablty wth Egeerg Alcatos Lecture 8 Professor Rav K. Iyer Uversty of Illos Iyer - Lecture 8 ECE 313 Fall 013 Today s Tocs Revew o Radom Varables Cumulatve Dstrbuto Fucto (CDF

More information

4 Round-Off and Truncation Errors

4 Round-Off and Truncation Errors HK Km Slgtly moded 3//9, /8/6 Frstly wrtte at Marc 5 4 Roud-O ad Trucato Errors Errors Roud-o Errors Trucato Errors Total Numercal Errors Bluders, Model Errors, ad Data Ucertaty Recallg, dv dt Δv v t Δt

More information

Power Flow S + Buses with either or both Generator Load S G1 S G2 S G3 S D3 S D1 S D4 S D5. S Dk. Injection S G1

Power Flow S + Buses with either or both Generator Load S G1 S G2 S G3 S D3 S D1 S D4 S D5. S Dk. Injection S G1 ower Flow uses wth ether or both Geerator Load G G G D D 4 5 D4 D5 ecto G Net Comple ower ecto - D D ecto s egatve sg at load bus = _ G D mlarl Curret ecto = G _ D At each bus coservato of comple power

More information

CHAPTER 6. d. With success = observation greater than 10, x = # of successes = 4, and

CHAPTER 6. d. With success = observation greater than 10, x = # of successes = 4, and CHAPTR 6 Secto 6.. a. We use the samle mea, to estmate the oulato mea µ. Σ 9.80 µ 8.407 7 ~ 7. b. We use the samle meda, 7 (the mddle observato whe arraged ascedg order. c. We use the samle stadard devato,

More information

Objectives of Multiple Regression

Objectives of Multiple Regression Obectves of Multple Regresso Establsh the lear equato that best predcts values of a depedet varable Y usg more tha oe eplaator varable from a large set of potetal predctors {,,... k }. Fd that subset of

More information

Unimodality Tests for Global Optimization of Single Variable Functions Using Statistical Methods

Unimodality 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 information

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Aalyss of Varace ad Desg of Exermets-I MODULE II LECTURE - GENERAL LINEAR HYPOTHESIS AND ANALYSIS OF VARIANCE Dr Shalabh Deartmet of Mathematcs ad Statstcs Ida Isttute of Techology Kaur Tukey s rocedure

More information

Maximum Likelihood Estimation

Maximum 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 information

A Remark on the Uniform Convergence of Some Sequences of Functions

A 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 information

ESS Line Fitting

ESS 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

Fractional Order Finite Difference Scheme For Soil Moisture Diffusion Equation And Its Applications

Fractional Order Finite Difference Scheme For Soil Moisture Diffusion Equation And Its Applications IOS Joural of Mathematcs (IOS-JM e-iss: 78-578. Volume 5, Issue 4 (Ja. - Feb. 3, PP -8 www.osrourals.org Fractoal Order Fte Dfferece Scheme For Sol Mosture Dffuso quato Ad Its Applcatos S.M.Jogdad, K.C.Takale,

More information

BASIC PRINCIPLES OF STATISTICS

BASIC PRINCIPLES OF STATISTICS BASIC PRINCIPLES OF STATISTICS PROBABILITY DENSITY DISTRIBUTIONS DISCRETE VARIABLES BINOMIAL DISTRIBUTION ~ B 0 0 umber of successes trals Pr E [ ] Var[ ] ; BINOMIAL DISTRIBUTION B7 0. B30 0.3 B50 0.5

More information

THE ROYAL STATISTICAL SOCIETY 2016 EXAMINATIONS SOLUTIONS HIGHER CERTIFICATE MODULE 5

THE ROYAL STATISTICAL SOCIETY 2016 EXAMINATIONS SOLUTIONS HIGHER CERTIFICATE MODULE 5 THE ROYAL STATISTICAL SOCIETY 06 EAMINATIONS SOLUTIONS HIGHER CERTIFICATE MODULE 5 The Socety s provdg these solutos to assst cadtes preparg for the examatos 07. The solutos are teded as learg ads ad should

More information

Chapter 3. Differentiation 3.3 Differentiation Rules

Chapter 3. Differentiation 3.3 Differentiation Rules 3.3 Dfferetato Rules 1 Capter 3. Dfferetato 3.3 Dfferetato Rules Dervatve of a Costat Fucto. If f as te costat value f(x) = c, te f x = [c] = 0. x Proof. From te efto: f (x) f(x + ) f(x) o c c 0 = 0. QED

More information

{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:

{ }{ ( )} (, ) = ( ) ( ) ( ) 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 information

G S Power Flow Solution

G 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 information

On a Semi-symmetric Non-metric Connection Satisfying the Schur`s Theorem on a Riemannian Manifold

On a Semi-symmetric Non-metric Connection Satisfying the Schur`s Theorem on a Riemannian Manifold O a Sem-symmetrc No-metrc oecto Satsfy te Scur`s Teorem o a emaa Mafod Ho Ta Yu Facuty of Matematcs, Km I Su versty, D.P..K Abstract: 99, Aace ad ae troduced te cocet of a sem-symmetrc o-metrc coecto[].

More information

8.1 Hashing Algorithms

8.1 Hashing Algorithms CS787: Advaced Algorthms Scrbe: Mayak Maheshwar, Chrs Hrchs Lecturer: Shuch Chawla Topc: Hashg ad NP-Completeess Date: September 21 2007 Prevously we looked at applcatos of radomzed algorthms, ad bega

More information

Chapter Two. An Introduction to Regression ( )

Chapter Two. An Introduction to Regression ( ) ubject: A Itroducto to Regresso Frst tage Chapter Two A Itroducto to Regresso (018-019) 1 pg. ubject: A Itroducto to Regresso Frst tage A Itroducto to Regresso Regresso aalss s a statstcal tool for the

More information

On the convergence of derivatives of Bernstein approximation

On 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 information