Stabilized Approximations of Strongly Continuous Semigroups

Size: px
Start display at page:

Download "Stabilized Approximations of Strongly Continuous Semigroups"

Transcription

1 Sabilized Approxiaios of Srogly Coiuous Seigroups Sarah McAlliser ad Frak Neubrader Sarah McAlliser, IBM T.J. Waso Research Ceer, Yorkow Heighs, NY 0 ad Frak Neubrader, Depare of Maheaics, Louisiaa Sae Uiversiy, Bao Rouge, LA 00 Absrac This paper iroduces sabilizaio echiques for irisically usable, high accuracy raioal approxiaio ehods for srogly coiuous seigroup. The ehods o oly sabilize he approxiaios, bu iprove heir speed of covergece by a agiude of up o /. Key words: srogly coiuous seigroups, Cheroff produc forula, raioal approxiaios of he expoeial, sabilizaio, ie discreizaio 000 MSC: A0, A, A0, L0, M0, M. Preliiaries The Lax-Richyer Equivalece Theore ad he Cheroff Produc Forula play a ipora role i he aalysis of approxiaio ehods for seigroups. Oe of he ai igredies i boh resuls is he sabiliy of he approxiaio ehod uder cosideraio. This paper deals wih varias of hese heores ha cover sabilizaio echiques for irisically usable approxiaio ehods. To pu he resuls io perspecive, recall he saee of he Lax-Richyer Equivalece Theore. Suppose ha A geeraes a srogly coiuous seigroup T ( ) oabaachspacex ad le {V (); [0,τ]} be a approxiaio schee of bouded liear operaors wih V (0) = I ha saisfies he cosisecy codiio V ()x x li = Ax 0 Eail address: sccalli@us.ib.co ad eubrad@ah.lsu.edu (Sarah McAlliser ad Frak Neubrader). Prepri subied o Elsevier March 00

2 for all x i a se D D(A) (he doai of A) haisdeseix. Theiwasshow by Lax ad Richyer i [] (wih a sroger cosisecy codiio) ad i fial for by Cheroff i [] ha he followig saees are equivale (for also [0]): (i) V (0) = I is sable; i.e, here exis ω 0adM,τ > 0 such ha V ( ) Me ω for each [0,τ]ad /τ. (ii) li V ( ) x = T ()x for all 0adx X. If ore iforaio o he approxiaio schees is available, error esiaes ca be give. Based o earlier work by Hersh ad Kao [], his was doe by Breer ad Thoée [] for approxiaio schees V () :=r(a) defied by raioal fucios r = P Q (P, Q polyoials) wih he followig wo properies. (a) The Maclauri series of r coicides wih he expoeial series for he firs ers; i.e., r(z) e z = O( z + )as z 0. (b) The fucio r is A-sable; i.e., r(z) forre(z) 0. The uber is called he approxiaio order of r = P Q.Iisawellkowresulof Padé [] ha deg(p )+deg(q) for all raioal approxiaios of he expoeial fucio. Raioal approxiaios of he expoeial for which =deg(p )+deg(q) are called Padé approxiaios. They are of he for r = P Q,where P (z) = deg(p ) deg(q) ( j)! deg(p )!!j!(deg(p ) j)! zj ad Q(z) = ( j)! deg(q)!!j!(deg(q) j)! ( z)j. As show i [], Padé schees are A-sable if ad oly if deg(q) deg(p ) deg(q). If r is a A-sable approxiaio of he expoeial of approxiaio order ad A geeraes a srogly coiuous seigroup, he he operaors V () := r(a) are cosise (see, for exaple, [], Th.). Exaples are (a) r(z) = z =+z + z + z + (Backward Euler, =, Padé ype), (b) r(z) = +z z = + z z =+z + + z + (Crak-Nicolso, =, Padé ype), (c) r(z) =c + c bz + c ( bz), where b = (+ ), c =,c = c, ad c = (Calaha, =, o Padé ype). If r is a A-sable raioal approxiaio of he expoeial ad V () := r(a) where A geeraes a srogly coiuous seigroup T (), he he error V ( ) x T ()x ca be esiaed (see Theore. below). The esiaes ca be iproved if r saisfies he codiio ( ); i.e., if (i) r(is) < fors R\{0} ad r( ) <, (ii) r(is) =e is+ψ(s) wih ψ(s) =O( s q+ ) for soe posiive ieger q as s 0 (iii) Reψ(s) γs p for s ad soe eve ieger p q +. If r is a A-sable raioal approxiaio of he expoeial ad if eiher (a) A geeraes a bouded aalyic seigroup or (b) A geeraes a srogly coiuous seigroup T ( ) o a Hilber space saisfyig T () e ω, he he approxiaio schee V () := r(a) is sable (see, [] ad []). I geeral, he followig resul is he bes possible. Theore. (Breer-Thoée) Le A geerae a srogly coiuous seigroup T ( ) wih T () Me ω.ifv () :=r(a) for soe A-sable raioal approxiaio r of he expoeial, he V ( ) is cosise bu ay o be sable. However, here are cosas C, κ such ha

3 V ( ) CM e ωκ (.) for all 0. Moreover, if r saisfies codiio ( ), he V ( ) CM (/ (q+)/p) e ωκ, (.) where p ad q are as i he defiiio of codiio ( ). The esiae (.) is sharp for he Crak-Nicolso schee r(z) = +z z ad A = d/dx o C 0 (R) (see []) or o he closure of D(A) il (R)(see []). Exaples of A-sable raioal approxiaios r = P Q saisfyig he codiio ( ) are Padé approxiaios wih deg(p )=deg(q) ordeg(p )=deg(q). For hese cases, he codiio ( ) holdswih q =deg(p )+deg(q) adp =deg(q) (.) (see [], [], []). Le r = P Q be of Padéypewihdeg(P )=deg(q). The i was show by A. Ashyralyev ad P. Sobolevskii ha V () = r(a) is sable if ad oly if A geeraes a srogly coiuous seigroup (see [], Theore II...) Observe ha he codiio deg(p )=deg(q) iplies ha he Padé approxiaio order =deg(p )+deg(q) us be odd; i.e., i geeral, approxiaio schees V () = r(a) ay o sable for Padé approxiaios of eve order (see also Table below). Exaples of raioal Padé approxiaios of odd order are he Backward Euler schee ( = ), he hird-order schee r(z) = +z z+z = C b + C z b wih C = + z.i ad b, =± i; ad he fifh-order Radau IIA schee r(z) = wih 0+z+z 0 z+z z poles a.,. ±.00i. For geeraors A of aalyic seigroups T ( ), he error aalysis of V ( ) x T ()x for raioal approxiaio schees V () = r(a) is well udersood. For such seigroups, he followig key esiaes are due o M. Crouzeix, S. Larsso, S. Piskarev, V. Thoée ad L.B. Wahlbi (see [] ad []). For ay A-sable raioal approxiaio schee r ad V () =r(a) oe has ha V ( ) is sable ad V ( ) x T ()x M A s x for s all x D(A s )(0 s ). Moreover, if deg(p ) < deg(q), he V ( ) x T ()x M x for all x X. These esiaes are based, o a large degree, o he Duford- Riesz fucioal calculus r(a)x := πi Γ r(λ)r(λ, A)xdλ,whereR(λ, A) :=(λi A). I he oaalyic case, he ai idea i he Breer-Thoée ad Hersh-Kao papers is o replace he Duford-Riesz fucioal calculus by he Hille-Phillips fucioal calculus r(a)x := T (s)x dµ(s), where T ( ) is a bouded, srogly coiuous seigroup 0 geeraed by A ad r(z) = ˆ 0 e sz dµ(s) (.) for Re(z) 0 ad soe oralized fucio µ of oal bouded variaio; see also [], []. To see how his fucioal calculus leads o error esiaes for V ( ) x T ()x wih V () := r(a), observe firs ha usig parial fracios, oe obais ha all A- sable raioal fucios r ca be represeed as i (.). Now ake a A-sable raioal approxiaio r of he expoeial of order such ha r(z) e z M z + for z sufficiely sall. If Re(z) 0, he bioial forula yields

4 r( z) e z = r( z) (e z ) = r( z)j (e z ) j r( z) e z M z + 0 as.observehar( z) = e sz dµ 0, (s), where µ, (s) ishe-h covoluio of µ( s ) wih iself ad ha ez = 0 e sz dh (s), where H is he oralized Heaviside fucio wih jup a ie 0. Sice ˆ 0 e sz d[µ, (s) H (s)] = r( z) e z 0 as, i is o surprisig ha µ, H i he L -or (wih precise error esiaes). I has o be poied ou ha he proof of hese facs is o-rivial ad a he hear of he Breer-Thoée ad Hersh-Kao papers; for Laplace rasfor proofs of hese saees, see [] or [0]. By he Hille-Phillips fucioal calculus, V ( ) x T ()x = = ˆ 0 ˆ 0 T (s)xd[µ, (s) H (s)] [µ, (s) H (s)]t (s)ax ds for all x D(A). Therefore, for bouded seigroups, V ( ) x T ()x C µ, H Ax 0 as for all x D(A), Togeher wih error-esiaes for he L -or of µ, H ad is aiderivaives, his is he ai srucure of he proof of he Breer-Thoée resul below (which appears i ore geeral for as Theore.. of [] ad is a refiee of groud-breakig resuls of Hersh ad Kao []; see also [0]). Theore. (Breer-Thoée) Le A be he geeraor of a srogly coiuous seigroup T ( ) wih T () Me ω.ler be a A-sable raioal approxiaio of he expoeial of approxiaio order ad defie V () :=r(a). Ifs =0,,,..., + wih s +, he here are posiive cosas c ad C (depedig oly o r) such ha V ( ) x T ()x CMe cω s ( ) β(s) A s x (.) for every 0, N, ad x D(A s ),where s β(s) := s + If s = +,he V ( ) x T ()x CMe cω + for every 0, N, ad x D(A + ). if 0 s< +, if + <s +. ( ) / l( +) A + x (.)

5 For sable raioal approxiaio schees he facor l( +) ca be reoved fro (.) above (his was show for = i [] ad, for arbirary, i [], Corollary.. or [], Corollary.). Moreover, if r saisfies codiio ( ), he Theore. holds wih β(s) replaced by β (s) :=s ( q q + +i 0, (s ( (q +)) q + )), (.) p where p ad q are as i he defiiio of codiio ( ). The followig able liss soe values of β(s); he sybol idicaes ha β(s) is udefied. Table. Values of he fucio β(s) D(A s ) = = = = = = = = s =0 s = s = s = s = s = s = s = s = For iiial daa x D(A s )wih0 s + ad wihou furher assupios o he A-sable raioal approxiaio r of he expoeial fucio defiig he approxiaio schee V () =r(a), he Baach space X, or o he seigroup T ( ), Table coais he followig egaive essages, predicig for V ( ) x T ()x (a) poeial divergece of order for soe x X (s =0), (b) a axial order of covergece of for x D(A) (s = ) idepede of he approxiaio order of he defiig uerical approxiaio schee r, (c) a axial order of covergece of s+ for x D(A s )(0 s +) o aer how large he approxiaio order. If r = P Q is of Padé ypewihdeg(p )=deg(q) ordeg(p )=deg(q), he he order of covergece of V ( ) x T ()x is β (s),whereβ (s)isasi(.)wihq = ad p =deg(q) (see (.)). Sice deg(p )=deg(q) ordeg(p )=deg(q), i follows ha for each here is exacly oe possible value for deg(q) (ad herefore also for deg(p )); aely, if =or if = +, he deg(q) = +. Thus, if is eve, he q = ad p =deg(q) =+; if is odd, he q = ad p =deg(q) =+. I paricular, if is odd ad r is of Padé ypewihdeg(p )=deg(q), he 0 β (s) =s +. If is eve ad if r is of Padé ypewihdeg(p )=deg(q), he 0 0

6 s( +) β (s) = ( +) s + if 0 s +, if + <s +. Table. Values of β (s) forpadé approxias wih deg(p )=deg(q) if is odd or deg(p )=deg(q) if is eve x D(A s ) = = = = = = = = s = s = s = s = s = s = s = s = 0 0 s = I suary, wihou furher assupios o he seigroup T ( ) geeraed by A (like aalyiciy) or o he space X (like X Hilber space ad T () ), approxiaio schees V ( ) = r( A) defied by A-sable raioal approxiaios r of he expoeial of approxiaio order are (a) guaraeed o be sable oly if r = P Q is of Padé ype, has a odd approxiaio order, ad deg(p )=deg(q) ; (b) covergig o T ()x for x D(A s )(0 s + ) axially wih order s if r saisfies codiio ( ) ad of order s+ if o. The purpose of his paper is o show ha for geeraors A of srogly coiuous seigroups all approxiaio schees V () := r(a) defied by A-sable approxiaios r of he expoeial of approxiaio order ca be sabilized by +power-scaled Backward Euler seps ad ha he sabilizaio ehods o oly sabilize he approxiaios, bu iprove heir speed of covergece by a agiude of up o / ifr does o saisfy codiio ( ). A his poi we recoed a glace a Corollary. o ge a feel for he ypes of resuls o be discussed i he followig secio.. Sabilizaio If A geeraes a aalyic seigroup, he all approxiaio ehods V () = r(a) defied by A-sable raioal approxiaios r are sable ad V ( A) x e A x for all x X. However, if r( ) = (e.g., Crak-Nicolso), he he covergece ca be arbirarily slow for o-sooh x X. R. Raacher [], ad i fial for A. Hasbo [], cobied he high accuracy of he Crak-Nicolso schee wih he soohig properies of he Backward Euler schee o provide a reedy (for exesios, see []). More precisely, le r, r s be A-sable raioal approxiaio schees such ha he high accuracy approxiaio r has order, he sabilizig approxiaio r s order, ad

7 r s ( ) =0.The,ifA geeraes a aalyic seigroup ad R(λ, A) exiss for all λ 0, he sabilized approxiaio schee V (A) :=r( A) r s ( A) saisfies V (A)x T ()x C() x for all x X. I his secio we cosider geeraors A of srogly coiuous seigroups T ( ) oa Baach space X ha are o ecessarily aalyic. As explaied above, a approxiaio schee V () = r(a) defied via a A-sable raioal approxiaio r of he expoeial fucio is always cosise bu ay o be sable. The, by he Lax-Richyer Equivalece Theore, here exiss x X such ha V ( ) x does o coverge o T ()x. However, he followig resul shows ha all such approxiaio schees ca be sabilized by akig firs + odified Backward Euler seps (I α A) for soe appropriaely α chose 0 <α<. Moreover, i will be show ha by sabilizig V () =r(a) defied by a A-sable raioal approxiaio schee r ha does o saisfy he codiio ( ), he speed of covergece ca be iproved fro s+ s o s++ if x D(A s )ad s. Before providig a proof for arbirary s, we sar wih he case s =. Theore. Le A be he geeraor of a srogly coiuous seigroup T ( ) o a Baach space X. Defie V () :=r(a), wherer is a A-sable raioal approxiaio of he expoeial of approxiaio order ad a sabilizig schee W () := α R( α,a)=(i α A), (.) where α = +.The,forallτ>0here exiss a cosa M τ such ha V ( ) W ( ( ) + )+ x T ()x M τ ( x + Ax ) (.) for all [0,τ], allx D(A), ad all sufficiely large +.Furherore, li V ( ) W ( )+ x = T ()x for all x X. Proof. Choose M adω 0 such ha T () Me ω for 0. Sice AR(λ, A) = λr(λ, A) I, i follows ha AW () = ( α R( α,a) I) (a his poi, α > 0is α arbirary). By he Hille-Yosida Theore, λr(λ, A) Mλ λ ω for all λ>ω.moreover,if λ>ω 0 >ωhe λr(λ, A) Mλ λ ω M 0 for all λ>ω 0.Thus, W( ) M 0 ad AW ( α ) (M α 0 + ) for all >ω /α 0.Cosider ( ) V ( ) W ( )+ x T ()x (V ( ) T ())W ( )+ x + T ()(W ( )+ x x), ad observe ha li W ( )x = li α R( α α,a)x = x for all x X sice A α geeraes a srogly coiuous seigroup (see [], Proposiio..). The secod er of ( ) ca be esiaed by usig he bioial forula; i.e., T ()(W ( )+ x x) T () W ( )j W ( )x x.

8 Sice W ( ) M 0 for >ω /α 0 i follows ha T ()(W ( )+ x x) 0as. By (.) wih s = +,hereexisc, c > 0 such ha (V ( ) T ())W ( ( ) )+ x CMe cω (AW ( ))+ x ) ( ) α + (M 0 +) + x CMe cω ( α = M e cω ( ) α(+) x. If 0 <α< + he, for every τ>0, here exiss K τ > 0 such ha (V ( ) T ())W ( )+ x K τ x (.0) α(+) for all [0,τ]. This shows ha li V ( ) W ( )+ x = T ()x for all x X. Now suppose ha x D(A). The, for τ>0hereexissm > 0 such ha for α α T ()(W ( )+ x x) T () W ( )j W ( )x x = T () W ( ) j R( α M,A)Ax T () α ω Ax W ( ) j M e ω α ω Ax (M) j α M α Ax α >ω 0 ad [0,τ]. By ( ) ad (.0), for τ>0hereexissm τ such ha V ( ) W ( ( )+ x T ()x M τ + α(+) α α α ) ( x + Ax ) for [0,τ]adx D(A). Choosig α = + yields + α(+) =. α + The followig able copares he resuls of he Breer-Thoée Theore. o he covergece raes obaied by sabilizig A-sable raioal approxiaio schees. Noice ha for iiial daa x D(A), he rae of covergece + (for a sabilized schee) approaches for large. The leer c idicaes covergece of ukow speed. The colus arked Padé give he covergece raes for he Padé schees r = P Q,where deg(p )= ad deg(q) = + if he approxiaio rae is eve, ad deg(p )= + ad deg(q) = if he approxiaio rae is odd.

9 Table. Copariso of Theore. o Theore. x X, β(s) x X, β (s) x X x D(A),β(s) x D(A),β (s) x D(A) usabilized Padé sabilized usabilized Padé sabilized = c c = c = c c = c = c c = c = c c = 0 c As Theore. shows, approxiaios schees V () = r(a) defied via A-sable approxiaios of he expoeial of approxiaio order ca be sabilized by akig firs + odified Backward Euler seps (I α A). This resuls i a approxiaio α speed of order o D(A) (hereby iprovig he axial usabilized rae of /+ covergece of ), bu ay slow dow he order of approxiaio o D(A k )for / k +. To sabilize such ha he approxiaio order for iiial daa i D(A k )is approxiaely like for sufficiely large, we cosider sabilizers of he for k W () = where 0 <b <b < <b k <b k are arbirarily chose, α := k i= a i 0 α R( b i,a), (.) α +k+,ad a i := ( )k+i b k i k j= b j b i. (.) j i Wih his choice of he coefficies a i ad soe basic algebra, he sabilizer W () cabe represeed i he for W () :=( ) k+ kα A k k i= (b i I α A) + I. (.) Theore. Le A geerae a srogly coiuous seigroup T ( ) o a Baach space X. Defie V () :=r(a), where r is a A-sable raioal approxiaio of he expoeial of approxiaio order ad defie a sabilizig schee W ( ) as i (.). The li V ( ) W ( )+ x = T ()x (.) for all x X. Furherore,forallτ>0here exiss a cosa M τ such ha V ( ) W ( ( ) k +k+ )+ x T ()x M τ ( x + A k x ) (.) for all [0,τ], allx D(A k ), ad all sufficiely large +.

10 Proof. Fix k +adlew ( ) beasi(.). The, for all x X, W ()x x = kα A k k i= = kα A k k i= (b i I α A) x α R( b i α,a)x = k i= AR( b i α,a)x. By he Hille-Yosida Theore, here exiss ω>0 such ha λr(λ, A) M λ ω M for λ > ω. Thus, he operaors λr(λ, A) adar(λ, A) = λr(λ, A) I are uiforly bouded for λ>ω ad λr(λ, A)x x as for all x D(A) ad, herefore, for all x X = D(A). I follows ha li W ( )x = x for all x X. Moreover,if x D(A k )ad bi > ω, he here exiss C α > 0 such ha ( k ) W ()x x R( b i α,a) A k x (.) i= ( k ) (M) k α A k x = C kα A k x. b α i= i λ Also, if bi > ω, he here exiss M α > 0 such ha k a i AW () = α AR( b i α,a) k a i b i α α R( b i,a) I] α M α i= i= We proceed o esiae k a i. i= ( ) V ( ) W ( )+ x T ()x (V ( ) T ())W ( )+ x + T ()(W ( )+ x x). Sice he faily W ( ) is uiforly bouded ad li W ( )x = x for all x X, i follows ha he secod er i ( ) coverges o zero; i.e., T ()(W ( )+ x x) T () W ( )j W ( )x x 0 as for all x X. Moreover,ifx D(A k ), he i follows fro (.) ha here exiss C τ > 0 such ha T ()(W ( )+ x x) T () W ( )j W ( )x x C τ kα Ak x (.) for all 0 τ. To discuss he firs er i ( ), recall fro he Breer-Thoée esiae (.) wih s = + ha here are C, c > 0 such ha 0

11 (V ( ) T ()) W ( ( ) )+ x CMe cω (AW ( ))+ x ( ) ( ) ( CMe cω α + k + α M + a i ) x. Therefore, for all τ>0hereexissk τ > 0 such ha V ( ) W ( )+ x T ()W ( )+ x K τ x (.) α(+) for all [0,τ]. This shows ha li V ( ) W ( )+ x = T ()x for all x X if 0 <α< +. Furherore, cobiig (.) wih (.) yields ha for each τ > 0 here exiss M τ > 0 such ha, for [0,τ]adx D(A k ), V ( ) W ( )+ x T ()x M τ ( + α(+) kα )( x + Ak x ). Choosig α = +k+ yields + =.Thus, α(+) kα +k+ k V ( ) W ( ( ) k k++ )+ x T ()x M τ ( x + A k x ). (.) As Theore. shows, by codiioig raioal approxiaio schees of order by firs applyig k-sabilizers of he for (.) yields (a) covergece for all x X ad (b) k covergece of order k++ for x D(Ak ). If k +k, he k-sabilizig iproves upo he rae of covergece o D(A k ) prediced by he Breer-Thoée Theore.. Ideed, if k + k, he + >kad β(k) =k k k++ k if. Exaple. (-sabilizaio) Le A geerae a srogly coiuous coracio seigroup ad V () :=r(a), wherer is a A-sable raioal approxiaio of he expoeial of approxiaio order. Cosider he -sabilizer W () := b (I α A) + b (I α A), b b b b b b where α = + ad 0 < b < b (see (.) ad (.)). By keepig rack of he cosas, i follows fro he proof of he previous heore for k =ha V ( ) W ( ( τ ) )+ + x T ()x M τ ( x + A x ), for all [0,τ], wherem τ := ax(k τ,n τ ) for N τ = i= b ( +b b b ) + b b b +b b b ad K τ = Cτ + ( b +b b b ) + (wih C depedig oly o r). Therefore he choice of bi affecs he size of M τ. Oe ca show ha M τ + ax(cτ,) for appropriaely chose b i (wih 0 <b <b 0 as ). The sabilizaios of raioal approxiaio schees for srogly coiuous seigroups (Theores. ad.) ad aalyic seigroups (see[], [], []) have a absrac exesio ha ca be viewed as a sabilized versio of he Lax-Richyer Equivalece

12 Theore (i.e., a sabilized Cheroff Produc Forula). The proof of he sabilized Cheroff Produc Forula requires a sabilized versio of he Troer-Kao Theore; he sabilized Troer-Kao ad Lax-Cheroff Theores were firs proved (for he case j =) by Y. Zhuag i []. Proofs of he sabilized versios of he Troer-Kao Theore ad he Cheroff Produc Forula ay be foud i []. Theore. (Sabilized Cheroff Produc Forula.) Le A geerae a srogly coiuous seigroup T ( ) wih T () Me ω for all 0 ad le {V () : [0,τ]} ad {W () : [0,τ]} be srogly coiuous failies of bouded liear operaors o X saisfyig W (0) = V (0) = I. IfV is cosise ad V ( ) j W ( )j Me ω ( j, [0,τ]), he li W ( )j V ( ) j x = T ()x for all x X uiforly o [0,τ]. To see how his resul ca be used here, recall fro he Breer-Thoée Theore. ha if V () := r(a) for a A-sable raioal approxiaio r of he expoeial of approxiaio order, he V ( ) x T ()x Ce c Ax for x D(A). Thus V () x Ce c Ax + Me ω x ad here exi M 0,ω 0 such ha, for all [0,τ], x D(A), ad > V ( ) c x Ce ( ) Ax + Me ω x M 0 e ω0 ( Ax + x ). Defie W () :=(I A). The, as show i he begiig of he proof of Theore., here exiss M>0such ha W ( ) M ad AW ( ) M for all N ad >0. Thus, here exiss M τ > 0 such ha V ( ) W ( ) M τe ω0 for ad [0,τ]. Sice V ad W coue, Theore. yields he followig corollary. Corollary. Le A be he geeraor of a srogly coiuous seigroup T ( ) ad le V () :=r(a), where r is a A-sable raioal approxiaio of he expoeial of approxiaio order. Defie W () := R(,A)=(I A). The li V ( ) W ( )x = T ()x for all x X uiforly o copac iervals. Refereces [] W. Ared, C. Bay, M. Hieber ad F. Neubrader, Vecor-Valued Laplace Trasfors ad Cauchy Probles. Birkhäuser, 00. [] A. Ashyralyev, P. Sobolevskii, Well-posedess of Parabolic Differece Equaios. Birkhäuser,. [] O. Axelsso, A class of A-sable ehods, BIT (), -. [] P. Breer ad V. Thoée, O raioal approxiaios of seigroups. SIAM J. Nuer. Aal. (), -.

13 [] P. Breer ad V. Thoée, Sabiliy ad covergece raes i L p for cerai differece schees. Mah. Scad (0), -. [] P. R. Cheroff, Produc Forulas, Noliear Seigroups, ad Addiio of Ubouded Operaors. Meoirs of he Aerica Mah. Soc. 0,. [] A. Chori, T. Hughes, M. McCracke ad J. Marsde, Produc forula ad uerical algorihs. Co. Pure Appl. Mah. (), 0-. [] M. Crouzeix, S. Larsso, S. Piskarev ad V. Thoeée, The Sabiliy of Raioal Approxiaios of Seigroups. BIT (), -. [] B.L.Ehle,O Padé approxiaios o he expoeial fucio ad A-sable ehods for he uerical soluio of iiial value probles, SIAM J. Mah. Aal., (), -0. [0] K.-J. Egel, R. Nagel, Oe-Paraeer Seigroups for Liear Evoluio Equaios. Spriger,. [] S. Flory, O he Sabilizaio ad Regularizaio of Approxiaio Schees for Seigroups, Disseraio, Louisiaa Sae Uiversiy, 00. [] S. Flory, F. Neubrader, ad L. Weis, Cosisecy ad Sabilizaio of Padé Approxiaio Schees for C 0 -seigroups, I: Evoluio Equaios: Applicaios o Physics, Idusry, Life Scieces ad Ecooics, Birkhäuser, 00. [] S. Flory, F. Neubrader ad Y. Zhuag, O he regularizaio ad sabilizaio of approxiaio schees for C 0 -seigroups. Parial Differeial Equaios ad Specral Theory, edied by M. Deuh ad B.-W. Schulze, OT. Birkhäuser, 00. [] E. Hairer ad G. Waer, Solvig Ordiary Differeial Equaios II - Siff ad Differeial- Algebraic Probles. Spriger,. [] A. Hasbo, Nosooh daa error esiaes for daped sigle sep ehods for parabolic equaios i Baach spaces. Calcolo (), -0. [] R. Hersh ad T. Kao, High-accuracy sable differece schees for well-posed iiial value probles. SIAM J. Nuer. Aal. (), o., 0-. [] M. Kovács, O he Qualiaive Properies ad Covergece of Tie-discreizaio Mehods for Seigroups, Disseraio, Louisiaa Sae Uiversiy, 00. [] M. Kovács, O he Covergece of Raioal Approxiaios of Seigroups o Ierediae Spaces, Maheaics of Copuaio (00), -. [] M. Kovács, O posiiviy, shape, ad or-boud preservaio of ie-seppig ehods for seigroups, J. Mah. Aal. Appl. 0 (00), -. [0] M. Kovács ad F. Neubrader, O he iverse Laplace-Sieljes rasfor of A-sable raioal fucios, Prepri, 00. [] P. Lax ad R. Richyer, Survey of sabiliy of liear fiie differece equaios. Co. Pure Appl. Mah. (), -. [] S. Larsso, V. Thoée ad L. B. Wahlbi, Fiie-elee ehods for a srogly daped wave equaio. IMAJ.Nuer.Aal.II (), -. [] M. Luski ad R. Raacher, O he soohig propery of he Crak-Nicolso schee. Applicable Aalysis, (), -. [] S. McAlliser, F. Neubrader ad Y. Zhuag, Sabilized Versios of he Troer-Kao Theore ad he Cheroff Produc Forula, i preparaio. [] M. H. Padé, Sur répreseaio approchée d ue focio par des fracioelles. A. de l Ecole Norale Superieure (). [] R. Raacher, Fiie elee soluio of diffusio probles wih irregular daa. Nuer. Mah. (), 0. [] V. Thoée, Galerki Fiie Elee Mehods for Parabolic Probles. Spriger,. [] Y. Zhuag, Classically Usable Approxiaios for Liear Evoluio Equaios ad Applicaios. Disseraio, Louisiaa Sae Uiversiy, 000.

EXISTENCE THEORY OF RANDOM DIFFERENTIAL EQUATIONS D. S. Palimkar

EXISTENCE THEORY OF RANDOM DIFFERENTIAL EQUATIONS D. S. Palimkar Ieraioal Joural of Scieific ad Research Publicaios, Volue 2, Issue 7, July 22 ISSN 225-353 EXISTENCE THEORY OF RANDOM DIFFERENTIAL EQUATIONS D S Palikar Depare of Maheaics, Vasarao Naik College, Naded

More information

BE.430 Tutorial: Linear Operator Theory and Eigenfunction Expansion

BE.430 Tutorial: Linear Operator Theory and Eigenfunction Expansion BE.43 Tuorial: Liear Operaor Theory ad Eigefucio Expasio (adaped fro Douglas Lauffeburger) 9//4 Moivaig proble I class, we ecouered parial differeial equaios describig rasie syses wih cheical diffusio.

More information

CHARACTERIZATIONS OF THE NON-UNIFORM IN TIME ISS PROPERTY AND APPLICATIONS

CHARACTERIZATIONS OF THE NON-UNIFORM IN TIME ISS PROPERTY AND APPLICATIONS CHARACTERIZATIONS OF THE NON-UNIFORM IN TIME ISS PROPERTY AND APPLICATIONS I. Karafyllis ad J. Tsiias Depare of Maheaics, Naioal Techical Uiversiy of Ahes, Zografou Capus 578, Ahes, Greece Eail: jsi@ceral.ua.gr.

More information

In this section we will study periodic signals in terms of their frequency f t is said to be periodic if (4.1)

In this section we will study periodic signals in terms of their frequency f t is said to be periodic if (4.1) Fourier Series Iroducio I his secio we will sudy periodic sigals i ers o heir requecy is said o be periodic i coe Reid ha a sigal ( ) ( ) ( ) () or every, where is a uber Fro his deiiio i ollows ha ( )

More information

A quadratic convergence method for the management equilibrium model

A quadratic convergence method for the management equilibrium model (IJACSA Ieraioal Joural of Advaced Copuer Sciece ad Applicaios Vol No 9 03 A quadraic covergece ehod for he aagee equilibriu odel Jiayi Zhag Feixia School Liyi Uiversiy Feixia Shadog PRChia Absrac i his

More information

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4)

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4) 7 Differeial equaios Review Solve by he mehod of udeermied coefficies ad by he mehod of variaio of parameers (4) y y = si Soluio; we firs solve he homogeeous equaio (4) y y = 4 The correspodig characerisic

More information

Comparison between Fourier and Corrected Fourier Series Methods

Comparison between Fourier and Corrected Fourier Series Methods Malaysia Joural of Mahemaical Scieces 7(): 73-8 (13) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Joural homepage: hp://eispem.upm.edu.my/oural Compariso bewee Fourier ad Correced Fourier Series Mehods 1

More information

A TAUBERIAN THEOREM FOR THE WEIGHTED MEAN METHOD OF SUMMABILITY

A TAUBERIAN THEOREM FOR THE WEIGHTED MEAN METHOD OF SUMMABILITY U.P.B. Sci. Bull., Series A, Vol. 78, Iss. 2, 206 ISSN 223-7027 A TAUBERIAN THEOREM FOR THE WEIGHTED MEAN METHOD OF SUMMABILITY İbrahim Çaak I his paper we obai a Tauberia codiio i erms of he weighed classical

More information

L-functions and Class Numbers

L-functions and Class Numbers L-fucios ad Class Numbers Sude Number Theory Semiar S. M.-C. 4 Sepember 05 We follow Romyar Sharifi s Noes o Iwasawa Theory, wih some help from Neukirch s Algebraic Number Theory. L-fucios of Dirichle

More information

1 Notes on Little s Law (l = λw)

1 Notes on Little s Law (l = λw) Copyrigh c 26 by Karl Sigma Noes o Lile s Law (l λw) We cosider here a famous ad very useful law i queueig heory called Lile s Law, also kow as l λw, which assers ha he ime average umber of cusomers i

More information

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory Liear Time-Ivaria Sysems (LTI Sysems) Oulie Basic Sysem Properies Memoryless ad sysems wih memory (saic or dyamic) Causal ad o-causal sysems (Causaliy) Liear ad o-liear sysems (Lieariy) Sable ad o-sable

More information

Numerical Solution of Parabolic Volterra Integro-Differential Equations via Backward-Euler Scheme

Numerical Solution of Parabolic Volterra Integro-Differential Equations via Backward-Euler Scheme America Joural of Compuaioal ad Applied Maemaics, (6): 77-8 DOI:.59/.acam.6. Numerical Soluio of Parabolic Volerra Iegro-Differeial Equaios via Bacward-Euler Sceme Ali Filiz Deparme of Maemaics, Ada Mederes

More information

Mean Square Convergent Finite Difference Scheme for Stochastic Parabolic PDEs

Mean Square Convergent Finite Difference Scheme for Stochastic Parabolic PDEs America Joural of Compuaioal Mahemaics, 04, 4, 80-88 Published Olie Sepember 04 i SciRes. hp://www.scirp.org/joural/ajcm hp://dx.doi.org/0.436/ajcm.04.4404 Mea Square Coverge Fiie Differece Scheme for

More information

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws.

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws. Limi of a fucio.. Oe-Sided..... Ifiie limis Verical Asympoes... Calculaig Usig he Limi Laws.5 The Squeeze Theorem.6 The Precise Defiiio of a Limi......7 Coiuiy.8 Iermediae Value Theorem..9 Refereces..

More information

Available online at J. Math. Comput. Sci. 4 (2014), No. 4, ISSN:

Available online at   J. Math. Comput. Sci. 4 (2014), No. 4, ISSN: Available olie a hp://sci.org J. Mah. Compu. Sci. 4 (2014), No. 4, 716-727 ISSN: 1927-5307 ON ITERATIVE TECHNIQUES FOR NUMERICAL SOLUTIONS OF LINEAR AND NONLINEAR DIFFERENTIAL EQUATIONS S.O. EDEKI *, A.A.

More information

TAKA KUSANO. laculty of Science Hrosh tlnlersty 1982) (n-l) + + Pn(t)x 0, (n-l) + + Pn(t)Y f(t,y), XR R are continuous functions.

TAKA KUSANO. laculty of Science Hrosh tlnlersty 1982) (n-l) + + Pn(t)x 0, (n-l) + + Pn(t)Y f(t,y), XR R are continuous functions. Iera. J. Mah. & Mah. Si. Vol. 6 No. 3 (1983) 559-566 559 ASYMPTOTIC RELATIOHIPS BETWEEN TWO HIGHER ORDER ORDINARY DIFFERENTIAL EQUATIONS TAKA KUSANO laculy of Sciece Hrosh llersy 1982) ABSTRACT. Some asympoic

More information

E will be denoted by n

E will be denoted by n JASEM ISSN 9-8362 All rigs reserved Full-ex Available Olie a p:// wwwbiolieorgbr/ja J Appl Sci Eviro Mg 25 Vol 9 3) 3-36 Corollabiliy ad Null Corollabiliy of Liear Syses * DAVIES, I; 2 JACKREECE, P Depare

More information

Ergodic Control of Semilinear Stochastic Equations and the Hamilton Jacobi Equation*

Ergodic Control of Semilinear Stochastic Equations and the Hamilton Jacobi Equation* Joural of Maheaical Aalysis ad Applicaios 34, 5963 999 Aricle ID jaa9996387, available olie a hp:wwwidealibraryco o Ergodic Corol of Seiliear Sochasic Equaios ad he ailojacobi Equaio* Beiai Goldys School

More information

VISCOSITY APPROXIMATION TO COMMON FIXED POINTS OF kn- LIPSCHITZIAN NONEXPANSIVE MAPPINGS IN BANACH SPACES

VISCOSITY APPROXIMATION TO COMMON FIXED POINTS OF kn- LIPSCHITZIAN NONEXPANSIVE MAPPINGS IN BANACH SPACES Joral o Maheaical Scieces: Advaces ad Alicaios Vole Nber 9 Pages -35 VISCOSIY APPROXIMAION O COMMON FIXED POINS OF - LIPSCHIZIAN NONEXPANSIVE MAPPINGS IN BANACH SPACES HONGLIANG ZUO ad MIN YANG Deare o

More information

N! AND THE GAMMA FUNCTION

N! AND THE GAMMA FUNCTION N! AND THE GAMMA FUNCTION Cosider he produc of he firs posiive iegers- 3 4 5 6 (-) =! Oe calls his produc he facorial ad has ha produc of he firs five iegers equals 5!=0. Direcly relaed o he discree! fucio

More information

Energy Density / Energy Flux / Total Energy in 1D. Key Mathematics: density, flux, and the continuity equation.

Energy Density / Energy Flux / Total Energy in 1D. Key Mathematics: density, flux, and the continuity equation. ecure Phys 375 Eergy Desiy / Eergy Flu / oal Eergy i D Overview ad Moivaio: Fro your sudy of waves i iroducory physics you should be aware ha waves ca raspor eergy fro oe place o aoher cosider he geeraio

More information

Online Supplement to Reactive Tabu Search in a Team-Learning Problem

Online Supplement to Reactive Tabu Search in a Team-Learning Problem Olie Suppleme o Reacive abu Search i a eam-learig Problem Yueli She School of Ieraioal Busiess Admiisraio, Shaghai Uiversiy of Fiace ad Ecoomics, Shaghai 00433, People s Republic of Chia, she.yueli@mail.shufe.edu.c

More information

ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003

ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003 ODEs II, Suppleme o Lecures 6 & 7: The Jorda Normal Form: Solvig Auoomous, Homogeeous Liear Sysems April 2, 23 I his oe, we describe he Jorda ormal form of a marix ad use i o solve a geeral homogeeous

More information

FIXED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE

FIXED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE Mohia & Samaa, Vol. 1, No. II, December, 016, pp 34-49. ORIGINAL RESEARCH ARTICLE OPEN ACCESS FIED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE 1 Mohia S. *, Samaa T. K. 1 Deparme of Mahemaics, Sudhir Memorial

More information

Math 6710, Fall 2016 Final Exam Solutions

Math 6710, Fall 2016 Final Exam Solutions Mah 67, Fall 6 Fial Exam Soluios. Firs, a sude poied ou a suble hig: if P (X i p >, he X + + X (X + + X / ( evaluaes o / wih probabiliy p >. This is roublesome because a radom variable is supposed o be

More information

Optimum design of complementary transient experiments for estimating thermal properties

Optimum design of complementary transient experiments for estimating thermal properties Opiu desig of copleeary rasie experies for esiaig heral properies Jaes V. Beck*, Filippo de Moe, Doald E. Aos *Depare of Mechaical Egieerig, Michiga Sae Uiversiy, USA Depare of Idusrial ad Iforaio Egieerig

More information

Extremal graph theory II: K t and K t,t

Extremal graph theory II: K t and K t,t Exremal graph heory II: K ad K, Lecure Graph Theory 06 EPFL Frak de Zeeuw I his lecure, we geeralize he wo mai heorems from he las lecure, from riagles K 3 o complee graphs K, ad from squares K, o complee

More information

Lecture 15 First Properties of the Brownian Motion

Lecture 15 First Properties of the Brownian Motion Lecure 15: Firs Properies 1 of 8 Course: Theory of Probabiliy II Term: Sprig 2015 Isrucor: Gorda Zikovic Lecure 15 Firs Properies of he Browia Moio This lecure deals wih some of he more immediae properies

More information

Bernstein Direct Method for Solving. Variational Problems

Bernstein Direct Method for Solving. Variational Problems Ieraioal Maheaical Foru, 5,, o. 48, 35-37 Bersei Direc Mehod for Solvig Variaioal Probles Sadeep Dixi*, Viee K. Sigh**, Ai K. Sigh*, O P. Sigh* *Depare of Applied Maheaics, Isiue of echology, Baaras Hidu

More information

Supplement for SADAGRAD: Strongly Adaptive Stochastic Gradient Methods"

Supplement for SADAGRAD: Strongly Adaptive Stochastic Gradient Methods Suppleme for SADAGRAD: Srogly Adapive Sochasic Gradie Mehods" Zaiyi Che * 1 Yi Xu * Ehog Che 1 iabao Yag 1. Proof of Proposiio 1 Proposiio 1. Le ɛ > 0 be fixed, H 0 γi, γ g, EF (w 1 ) F (w ) ɛ 0 ad ieraio

More information

Research Article A Generalized Nonlinear Sum-Difference Inequality of Product Form

Research Article A Generalized Nonlinear Sum-Difference Inequality of Product Form Joural of Applied Mahemaics Volume 03, Aricle ID 47585, 7 pages hp://dx.doi.org/0.55/03/47585 Research Aricle A Geeralized Noliear Sum-Differece Iequaliy of Produc Form YogZhou Qi ad Wu-Sheg Wag School

More information

The Solution of the One Species Lotka-Volterra Equation Using Variational Iteration Method ABSTRACT INTRODUCTION

The Solution of the One Species Lotka-Volterra Equation Using Variational Iteration Method ABSTRACT INTRODUCTION Malaysia Joural of Mahemaical Scieces 2(2): 55-6 (28) The Soluio of he Oe Species Loka-Volerra Equaio Usig Variaioal Ieraio Mehod B. Baiha, M.S.M. Noorai, I. Hashim School of Mahemaical Scieces, Uiversii

More information

EXAMPLE SHEET B (Partial Differential Equations) SPRING 2014

EXAMPLE SHEET B (Partial Differential Equations) SPRING 2014 Copuaioal Mechaics Eaples B - David Apsle EXAMPLE SHEET B Parial Differeial Equaios SPRING 0 B. Solve he parial differeial equaio 0 0 o, u u u u B. Classif he followig d -order PDEs as hperbolic, parabolic

More information

The Importance of Ordering the Number of Lattice Points Inside a Rational Polyhedron Using Generating Functions

The Importance of Ordering the Number of Lattice Points Inside a Rational Polyhedron Using Generating Functions Ieraioal Joural of Copuer Sciece ad Elecroics Egieerig (IJCSEE Volue Issue ( ISSN 48 (Olie he Iporace of Orderig he Nuber of Laice ois Iside a Raioal olyhedro Usig Geeraig Fucios Halil Sopce Absrac I pure

More information

Calculus BC 2015 Scoring Guidelines

Calculus BC 2015 Scoring Guidelines AP Calculus BC 5 Scorig Guidelies 5 The College Board. College Board, Advaced Placeme Program, AP, AP Ceral, ad he acor logo are regisered rademarks of he College Board. AP Ceral is he official olie home

More information

Approximately Quasi Inner Generalized Dynamics on Modules. { } t t R

Approximately Quasi Inner Generalized Dynamics on Modules. { } t t R Joural of Scieces, Islamic epublic of Ira 23(3): 245-25 (22) Uiversiy of Tehra, ISSN 6-4 hp://jscieces.u.ac.ir Approximaely Quasi Ier Geeralized Dyamics o Modules M. Mosadeq, M. Hassai, ad A. Nikam Deparme

More information

MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES. Boundary Value Problem for the Higher Order Equation with Fractional Derivative

MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES. Boundary Value Problem for the Higher Order Equation with Fractional Derivative Malaysia Joural of Maheaical Scieces 7(): 3-7 (3) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Joural hoepage: hp://eispe.up.edu.y/joural Boudary Value Proble for he Higher Order Equaio wih Fracioal Derivaive

More information

Review Exercises for Chapter 9

Review Exercises for Chapter 9 0_090R.qd //0 : PM Page 88 88 CHAPTER 9 Ifiie Series I Eercises ad, wrie a epressio for he h erm of he sequece..,., 5, 0,,,, 0,... 7,... I Eercises, mach he sequece wih is graph. [The graphs are labeled

More information

Sliding Mode Control for Robust Stabilization of Uncertain Input-Delay Systems

Sliding Mode Control for Robust Stabilization of Uncertain Input-Delay Systems 98 CSE: he siue of Corol, uoaio ad Syses Egieers, KORE Vol, No, Jue, Slidig Mode Corol for Robus Sabilizaio of Ucerai pu-delay Syses Youg-Hoo Roh ad Ju-Ho Oh bsrac: his paper is cocered wih a delay-depede

More information

Improvement Over General And Wider Class of Estimators Using Ranked Set Sampling

Improvement Over General And Wider Class of Estimators Using Ranked Set Sampling ITERATIOAL JOURAL OF SIETIFI & TEOLOG RESEAR VOLUME ISSUE 7 AUGUST ISS 77-866 Iprovee Over Geeral Ad ider lass of Esiaors Usi Raked Se Sapli V L Madoara iu Meha Raka Absrac: I his paper Iprovee over eeral

More information

Dynamic h-index: the Hirsch index in function of time

Dynamic h-index: the Hirsch index in function of time Dyamic h-idex: he Hirsch idex i fucio of ime by L. Egghe Uiversiei Hassel (UHassel), Campus Diepebeek, Agoralaa, B-3590 Diepebeek, Belgium ad Uiversiei Awerpe (UA), Campus Drie Eike, Uiversieisplei, B-260

More information

Lecture 9: Polynomial Approximations

Lecture 9: Polynomial Approximations CS 70: Complexiy Theory /6/009 Lecure 9: Polyomial Approximaios Isrucor: Dieer va Melkebeek Scribe: Phil Rydzewski & Piramaayagam Arumuga Naiar Las ime, we proved ha o cosa deph circui ca evaluae he pariy

More information

On Another Type of Transform Called Rangaig Transform

On Another Type of Transform Called Rangaig Transform Ieraioal Joural of Parial Differeial Equaios ad Applicaios, 7, Vol 5, No, 4-48 Available olie a hp://pubssciepubcom/ijpdea/5//6 Sciece ad Educaio Publishig DOI:69/ijpdea-5--6 O Aoher Type of Trasform Called

More information

Notes 03 largely plagiarized by %khc

Notes 03 largely plagiarized by %khc 1 1 Discree-Time Covoluio Noes 03 largely plagiarized by %khc Le s begi our discussio of covoluio i discree-ime, sice life is somewha easier i ha domai. We sar wih a sigal x[] ha will be he ipu io our

More information

An interesting result about subset sums. Nitu Kitchloo. Lior Pachter. November 27, Abstract

An interesting result about subset sums. Nitu Kitchloo. Lior Pachter. November 27, Abstract A ieresig resul abou subse sums Niu Kichloo Lior Pacher November 27, 1993 Absrac We cosider he problem of deermiig he umber of subses B f1; 2; : : :; g such ha P b2b b k mod, where k is a residue class

More information

APPROXIMATE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS WITH UNCERTAINTY

APPROXIMATE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS WITH UNCERTAINTY APPROXIMATE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS WITH UNCERTAINTY ZHEN-GUO DENG ad GUO-CHENG WU 2, 3 * School of Mahemaics ad Iformaio Sciece, Guagi Uiversiy, Naig 534, PR Chia 2 Key Laboraory

More information

Section 8 Convolution and Deconvolution

Section 8 Convolution and Deconvolution APPLICATIONS IN SIGNAL PROCESSING Secio 8 Covoluio ad Decovoluio This docume illusraes several echiques for carryig ou covoluio ad decovoluio i Mahcad. There are several operaors available for hese fucios:

More information

12 Getting Started With Fourier Analysis

12 Getting Started With Fourier Analysis Commuicaios Egieerig MSc - Prelimiary Readig Geig Sared Wih Fourier Aalysis Fourier aalysis is cocered wih he represeaio of sigals i erms of he sums of sie, cosie or complex oscillaio waveforms. We ll

More information

INTEGER INTERVAL VALUE OF NEWTON DIVIDED DIFFERENCE AND FORWARD AND BACKWARD INTERPOLATION FORMULA

INTEGER INTERVAL VALUE OF NEWTON DIVIDED DIFFERENCE AND FORWARD AND BACKWARD INTERPOLATION FORMULA Volume 8 No. 8, 45-54 ISSN: 34-3395 (o-lie versio) url: hp://www.ijpam.eu ijpam.eu INTEGER INTERVAL VALUE OF NEWTON DIVIDED DIFFERENCE AND FORWARD AND BACKWARD INTERPOLATION FORMULA A.Arul dass M.Dhaapal

More information

Fermat Numbers in Multinomial Coefficients

Fermat Numbers in Multinomial Coefficients 1 3 47 6 3 11 Joural of Ieger Sequeces, Vol. 17 (014, Aricle 14.3. Ferma Numbers i Muliomial Coefficies Shae Cher Deparme of Mahemaics Zhejiag Uiversiy Hagzhou, 31007 Chia chexiaohag9@gmail.com Absrac

More information

A note on Generalized Hermite polynomials

A note on Generalized Hermite polynomials INTERNATIONAL JOURNAL OF APPLIED MATHEMATICS AND INFORMATICS Volue 8 14 A oe o Geeralized Herie poloials Cleee Cesarao Absrac B sarig fro he geeral heor of he oevariable Herie poloials we will iroduce

More information

A NEW NUMERICAL TECHNIQUE FOR SOLVING FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS

A NEW NUMERICAL TECHNIQUE FOR SOLVING FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS A NEW NUMERICAL TECHNIQUE FOR SOLVING FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS Oer Aca *, Duiru Baleau,3 Depare of Maheaics, Faculy of Ar ad Sciece, Siir Uiversiy, 56, Siir, Turey Depare of Maheaics ad

More information

Completeness of Random Exponential System in Half-strip

Completeness of Random Exponential System in Half-strip 23-24 Prepri for School of Mahemaical Scieces, Beijig Normal Uiversiy Compleeess of Radom Expoeial Sysem i Half-srip Gao ZhiQiag, Deg GuaTie ad Ke SiYu School of Mahemaical Scieces, Laboraory of Mahemaics

More information

Homotopy Analysis Method for Solving Fractional Sturm-Liouville Problems

Homotopy Analysis Method for Solving Fractional Sturm-Liouville Problems Ausralia Joural of Basic ad Applied Scieces, 4(1): 518-57, 1 ISSN 1991-8178 Homoopy Aalysis Mehod for Solvig Fracioal Surm-Liouville Problems 1 A Neamay, R Darzi, A Dabbaghia 1 Deparme of Mahemaics, Uiversiy

More information

LIMITS OF FUNCTIONS (I)

LIMITS OF FUNCTIONS (I) LIMITS OF FUNCTIO (I ELEMENTARY FUNCTIO: (Elemeary fucios are NOT piecewise fucios Cosa Fucios: f(x k, where k R Polyomials: f(x a + a x + a x + a x + + a x, where a, a,..., a R Raioal Fucios: f(x P (x,

More information

Some Properties of Semi-E-Convex Function and Semi-E-Convex Programming*

Some Properties of Semi-E-Convex Function and Semi-E-Convex Programming* The Eighh Ieraioal Symposium o Operaios esearch ad Is Applicaios (ISOA 9) Zhagjiajie Chia Sepember 2 22 29 Copyrigh 29 OSC & APOC pp 33 39 Some Properies of Semi-E-Covex Fucio ad Semi-E-Covex Programmig*

More information

ECE 570 Session 7 IC 752-E Computer Aided Engineering for Integrated Circuits. Transient analysis. Discuss time marching methods used in SPICE

ECE 570 Session 7 IC 752-E Computer Aided Engineering for Integrated Circuits. Transient analysis. Discuss time marching methods used in SPICE ECE 570 Sessio 7 IC 75-E Compuer Aided Egieerig for Iegraed Circuis Trasie aalysis Discuss ime marcig meods used i SPICE. Time marcig meods. Explici ad implici iegraio meods 3. Implici meods used i circui

More information

K3 p K2 p Kp 0 p 2 p 3 p

K3 p K2 p Kp 0 p 2 p 3 p Mah 80-00 Mo Ar 0 Chaer 9 Fourier Series ad alicaios o differeial equaios (ad arial differeial equaios) 9.-9. Fourier series defiiio ad covergece. The idea of Fourier series is relaed o he liear algebra

More information

A Note on Random k-sat for Moderately Growing k

A Note on Random k-sat for Moderately Growing k A Noe o Radom k-sat for Moderaely Growig k Ju Liu LMIB ad School of Mahemaics ad Sysems Sciece, Beihag Uiversiy, Beijig, 100191, P.R. Chia juliu@smss.buaa.edu.c Zogsheg Gao LMIB ad School of Mahemaics

More information

A note on deviation inequalities on {0, 1} n. by Julio Bernués*

A note on deviation inequalities on {0, 1} n. by Julio Bernués* A oe o deviaio iequaliies o {0, 1}. by Julio Berués* Deparameo de Maemáicas. Faculad de Ciecias Uiversidad de Zaragoza 50009-Zaragoza (Spai) I. Iroducio. Le f: (Ω, Σ, ) IR be a radom variable. Roughly

More information

STK4080/9080 Survival and event history analysis

STK4080/9080 Survival and event history analysis STK48/98 Survival ad eve hisory aalysis Marigales i discree ime Cosider a sochasic process The process M is a marigale if Lecure 3: Marigales ad oher sochasic processes i discree ime (recap) where (formally

More information

Extended Laguerre Polynomials

Extended Laguerre Polynomials I J Coemp Mah Scieces, Vol 7, 1, o, 189 194 Exeded Laguerre Polyomials Ada Kha Naioal College of Busiess Admiisraio ad Ecoomics Gulberg-III, Lahore, Pakisa adakhaariq@gmailcom G M Habibullah Naioal College

More information

Moment Generating Function

Moment Generating Function 1 Mome Geeraig Fucio m h mome m m m E[ ] x f ( x) dx m h ceral mome m m m E[( ) ] ( ) ( x ) f ( x) dx Mome Geeraig Fucio For a real, M () E[ e ] e k x k e p ( x ) discree x k e f ( x) dx coiuous Example

More information

BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS

BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS Opimal ear Forecasig Alhough we have o meioed hem explicily so far i he course, here are geeral saisical priciples for derivig he bes liear forecas, ad

More information

Boundary-to-Displacement Asymptotic Gains for Wave Systems With Kelvin-Voigt Damping

Boundary-to-Displacement Asymptotic Gains for Wave Systems With Kelvin-Voigt Damping Boudary-o-Displaceme Asympoic Gais for Wave Sysems Wih Kelvi-Voig Dampig Iasso Karafyllis *, Maria Kooriaki ** ad Miroslav Krsic *** * Dep. of Mahemaics, Naioal Techical Uiversiy of Ahes, Zografou Campus,

More information

Actuarial Society of India

Actuarial Society of India Acuarial Sociey of Idia EXAMINAIONS Jue 5 C4 (3) Models oal Marks - 5 Idicaive Soluio Q. (i) a) Le U deoe he process described by 3 ad V deoe he process described by 4. he 5 e 5 PU [ ] PV [ ] ( e ).538!

More information

If boundary values are necessary, they are called mixed initial-boundary value problems. Again, the simplest prototypes of these IV problems are:

If boundary values are necessary, they are called mixed initial-boundary value problems. Again, the simplest prototypes of these IV problems are: 3. Iiial value problems: umerical soluio Fiie differeces - Trucaio errors, cosisecy, sabiliy ad covergece Crieria for compuaioal sabiliy Explici ad implici ime schemes Table of ime schemes Hyperbolic ad

More information

Lecture 15: Three-tank Mixing and Lead Poisoning

Lecture 15: Three-tank Mixing and Lead Poisoning Lecure 15: Three-ak Miig ad Lead Poisoig Eigevalues ad eigevecors will be used o fid he soluio of a sysem for ukow fucios ha saisfy differeial equaios The ukow fucios will be wrie as a 1 colum vecor [

More information

Common Fixed Point Theorem in Intuitionistic Fuzzy Metric Space via Compatible Mappings of Type (K)

Common Fixed Point Theorem in Intuitionistic Fuzzy Metric Space via Compatible Mappings of Type (K) Ieraioal Joural of ahemaics Treds ad Techology (IJTT) Volume 35 umber 4- July 016 Commo Fixed Poi Theorem i Iuiioisic Fuzzy eric Sace via Comaible aigs of Tye (K) Dr. Ramaa Reddy Assisa Professor De. of

More information

λiv Av = 0 or ( λi Av ) = 0. In order for a vector v to be an eigenvector, it must be in the kernel of λi

λiv Av = 0 or ( λi Av ) = 0. In order for a vector v to be an eigenvector, it must be in the kernel of λi Liear lgebra Lecure #9 Noes This week s lecure focuses o wha migh be called he srucural aalysis of liear rasformaios Wha are he irisic properies of a liear rasformaio? re here ay fixed direcios? The discussio

More information

arxiv:math/ v1 [math.fa] 1 Feb 1994

arxiv:math/ v1 [math.fa] 1 Feb 1994 arxiv:mah/944v [mah.fa] Feb 994 ON THE EMBEDDING OF -CONCAVE ORLICZ SPACES INTO L Care Schü Abrac. I [K S ] i wa how ha Ave ( i a π(i) ) π i equivale o a Orlicz orm whoe Orlicz fucio i -cocave. Here we

More information

David Randall. ( )e ikx. k = u x,t. u( x,t)e ikx dx L. x L /2. Recall that the proof of (1) and (2) involves use of the orthogonality condition.

David Randall. ( )e ikx. k = u x,t. u( x,t)e ikx dx L. x L /2. Recall that the proof of (1) and (2) involves use of the orthogonality condition. ! Revised April 21, 2010 1:27 P! 1 Fourier Series David Radall Assume ha u( x,) is real ad iegrable If he domai is periodic, wih period L, we ca express u( x,) exacly by a Fourier series expasio: ( ) =

More information

DIFFERENTIAL EQUATIONS

DIFFERENTIAL EQUATIONS DIFFERENTIAL EQUATIONS M.A. (Previous) Direcorae of Disace Educaio Maharshi Dayaad Uiversiy ROHTAK 4 Copyrigh 3, Maharshi Dayaad Uiversiy, ROHTAK All Righs Reserved. No par of his publicaio may be reproduced

More information

The Moment Approximation of the First Passage Time for the Birth Death Diffusion Process with Immigraton to a Moving Linear Barrier

The Moment Approximation of the First Passage Time for the Birth Death Diffusion Process with Immigraton to a Moving Linear Barrier America Joural of Applied Mahemaics ad Saisics, 015, Vol. 3, No. 5, 184-189 Available olie a hp://pubs.sciepub.com/ajams/3/5/ Sciece ad Educaio Publishig DOI:10.1691/ajams-3-5- The Mome Approximaio of

More information

10.3 Autocorrelation Function of Ergodic RP 10.4 Power Spectral Density of Ergodic RP 10.5 Normal RP (Gaussian RP)

10.3 Autocorrelation Function of Ergodic RP 10.4 Power Spectral Density of Ergodic RP 10.5 Normal RP (Gaussian RP) ENGG450 Probabiliy ad Saisics for Egieers Iroducio 3 Probabiliy 4 Probabiliy disribuios 5 Probabiliy Desiies Orgaizaio ad descripio of daa 6 Samplig disribuios 7 Ifereces cocerig a mea 8 Comparig wo reames

More information

Chapter 2: Time-Domain Representations of Linear Time-Invariant Systems. Chih-Wei Liu

Chapter 2: Time-Domain Representations of Linear Time-Invariant Systems. Chih-Wei Liu Caper : Time-Domai Represeaios of Liear Time-Ivaria Sysems Ci-Wei Liu Oulie Iroucio Te Covoluio Sum Covoluio Sum Evaluaio Proceure Te Covoluio Iegral Covoluio Iegral Evaluaio Proceure Iercoecios of LTI

More information

MATH 507a ASSIGNMENT 4 SOLUTIONS FALL 2018 Prof. Alexander. g (x) dx = g(b) g(0) = g(b),

MATH 507a ASSIGNMENT 4 SOLUTIONS FALL 2018 Prof. Alexander. g (x) dx = g(b) g(0) = g(b), MATH 57a ASSIGNMENT 4 SOLUTIONS FALL 28 Prof. Alexader (2.3.8)(a) Le g(x) = x/( + x) for x. The g (x) = /( + x) 2 is decreasig, so for a, b, g(a + b) g(a) = a+b a g (x) dx b so g(a + b) g(a) + g(b). Sice

More information

State-Space Model. In general, the dynamic equations of a lumped-parameter continuous system may be represented by

State-Space Model. In general, the dynamic equations of a lumped-parameter continuous system may be represented by Sae-Space Model I geeral, he dyaic equaio of a luped-paraeer coiuou ye ay be repreeed by x & f x, u, y g x, u, ae equaio oupu equaio where f ad g are oliear vecor-valued fucio Uig a liearized echique,

More information

Chapter 3 Moments of a Distribution

Chapter 3 Moments of a Distribution Chaper 3 Moes of a Disribuio Epecaio We develop he epecaio operaor i ers of he Lebesgue iegral. Recall ha he Lebesgue easure λ(a) for soe se A gives he legh/area/volue of he se A. If A = (3; 7), he λ(a)

More information

ON THE INITIAL-BOUNDARY VALUE PROBLEM FOR EQUATIONS OF ANOMALOUS DIFFUSION IN POLYMERS. D. A. Vorotnikov * Voronezh State University

ON THE INITIAL-BOUNDARY VALUE PROBLEM FOR EQUATIONS OF ANOMALOUS DIFFUSION IN POLYMERS. D. A. Vorotnikov * Voronezh State University UDC 579 ON THE INITIAL-BOUNDAY VALUE POBLEM FO EQUATIONS OF ANOMALOUS DIFFUSION IN POLYMES D A Voroikov * Voroezh Sae Uiversiy We sudy he syse of parial differeial equaios which describes he diffusio of

More information

Research Article Finite Difference and Sinc-Collocation Approximations to a Class of Fractional Diffusion-Wave Equations

Research Article Finite Difference and Sinc-Collocation Approximations to a Class of Fractional Diffusion-Wave Equations Joural of Applied Maheaics Volue 4, Aricle ID 5363, pages hp://dx.doi.org/.55/4/5363 Research Aricle Fiie Differece ad Sic-Collocaio Approxiaios o a Class of Fracioal Diffusio-Wave Equaios Zhi Mao,, Aiguo

More information

Approximating Solutions for Ginzburg Landau Equation by HPM and ADM

Approximating Solutions for Ginzburg Landau Equation by HPM and ADM Available a hp://pvamu.edu/aam Appl. Appl. Mah. ISSN: 193-9466 Vol. 5, No. Issue (December 1), pp. 575 584 (Previously, Vol. 5, Issue 1, pp. 167 1681) Applicaios ad Applied Mahemaics: A Ieraioal Joural

More information

Samuel Sindayigaya 1, Nyongesa L. Kennedy 2, Adu A.M. Wasike 3

Samuel Sindayigaya 1, Nyongesa L. Kennedy 2, Adu A.M. Wasike 3 Ieraioal Joural of Saisics ad Aalysis. ISSN 48-9959 Volume 6, Number (6, pp. -8 Research Idia Publicaios hp://www.ripublicaio.com The Populaio Mea ad is Variace i he Presece of Geocide for a Simple Birh-Deah-

More information

Existence Of Solutions For Nonlinear Fractional Differential Equation With Integral Boundary Conditions

Existence Of Solutions For Nonlinear Fractional Differential Equation With Integral Boundary Conditions Reserch Ivey: Ieriol Jourl Of Egieerig Ad Sciece Vol., Issue (April 3), Pp 8- Iss(e): 78-47, Iss(p):39-6483, Www.Reserchivey.Com Exisece Of Soluios For Nolier Frciol Differeil Equio Wih Iegrl Boudry Codiios,

More information

NONLINEAR SCHRÖDINGER EQUATION

NONLINEAR SCHRÖDINGER EQUATION TMA Nuerical Soluio of Differeial Equaios by Differece Mehods NONLINEAR SCHRÖDINGER EQUATION Group Cadidaes 8, 37 ad 55 Depare of Maheaical Scieces Norwegia Uiversiy of Sciece ad Techology April, 3 Absrac

More information

SOLVING OF THE FRACTIONAL NON-LINEAR AND LINEAR SCHRÖDINGER EQUATIONS BY HOMOTOPY PERTURBATION METHOD

SOLVING OF THE FRACTIONAL NON-LINEAR AND LINEAR SCHRÖDINGER EQUATIONS BY HOMOTOPY PERTURBATION METHOD SOLVING OF THE FRACTIONAL NON-LINEAR AND LINEAR SCHRÖDINGER EQUATIONS BY HOMOTOPY PERTURBATION METHOD DUMITRU BALEANU, ALIREZA K. GOLMANKHANEH,3, ALI K. GOLMANKHANEH 3 Deparme of Mahemaics ad Compuer Sciece,

More information

th m m m m central moment : E[( X X) ] ( X X) ( x X) f ( x)

th m m m m central moment : E[( X X) ] ( X X) ( x X) f ( x) 1 Trasform Techiques h m m m m mome : E[ ] x f ( x) dx h m m m m ceral mome : E[( ) ] ( ) ( x) f ( x) dx A coveie wa of fidig he momes of a radom variable is he mome geeraig fucio (MGF). Oher rasform echiques

More information

Pure Math 30: Explained!

Pure Math 30: Explained! ure Mah : Explaied! www.puremah.com 6 Logarihms Lesso ar Basic Expoeial Applicaios Expoeial Growh & Decay: Siuaios followig his ype of chage ca be modeled usig he formula: (b) A = Fuure Amou A o = iial

More information

Solutions to selected problems from the midterm exam Math 222 Winter 2015

Solutions to selected problems from the midterm exam Math 222 Winter 2015 Soluios o seleced problems from he miderm eam Mah Wier 5. Derive he Maclauri series for he followig fucios. (cf. Pracice Problem 4 log( + (a L( d. Soluio: We have he Maclauri series log( + + 3 3 4 4 +...,

More information

Math-303 Chapter 7 Linear systems of ODE November 16, Chapter 7. Systems of 1 st Order Linear Differential Equations.

Math-303 Chapter 7 Linear systems of ODE November 16, Chapter 7. Systems of 1 st Order Linear Differential Equations. Mah-33 Chaper 7 Liear sysems of ODE November 6, 7 Chaper 7 Sysems of s Order Liear Differeial Equaios saddle poi λ >, λ < Mah-33 Chaper 7 Liear sysems of ODE November 6, 7 Mah-33 Chaper 7 Liear sysems

More information

Compact Finite Difference Schemes for Solving a Class of Weakly- Singular Partial Integro-differential Equations

Compact Finite Difference Schemes for Solving a Class of Weakly- Singular Partial Integro-differential Equations Ma. Sci. Le. Vol. No. 53-0 (0 Maemaical Scieces Leers A Ieraioal Joural @ 0 NSP Naural Scieces Publisig Cor. Compac Fiie Differece Scemes for Solvig a Class of Weakly- Sigular Parial Iegro-differeial Equaios

More information

METHOD OF THE EQUIVALENT BOUNDARY CONDITIONS IN THE UNSTEADY PROBLEM FOR ELASTIC DIFFUSION LAYER

METHOD OF THE EQUIVALENT BOUNDARY CONDITIONS IN THE UNSTEADY PROBLEM FOR ELASTIC DIFFUSION LAYER Maerials Physics ad Mechaics 3 (5) 36-4 Received: March 7 5 METHOD OF THE EQUIVAENT BOUNDARY CONDITIONS IN THE UNSTEADY PROBEM FOR EASTIC DIFFUSION AYER A.V. Zemsov * D.V. Tarlaovsiy Moscow Aviaio Isiue

More information

Basic Results in Functional Analysis

Basic Results in Functional Analysis Preared by: F.. ewis Udaed: Suday, Augus 7, 4 Basic Resuls i Fucioal Aalysis f ( ): X Y is coiuous o X if X, (, ) z f( z) f( ) f ( ): X Y is uiformly coiuous o X if i is coiuous ad ( ) does o deed o. f

More information

Robust Dynamic Output Feedback Second-Order Sliding Mode Controller for Uncertain Systems

Robust Dynamic Output Feedback Second-Order Sliding Mode Controller for Uncertain Systems Ieraioal Joural of Corol, Auoaio, ad Syses (3 (5:878-884 DOI.7/s555--7-9 ISSN:598-6446 eissn:5-49 hp://www.spriger.co/555 Robus Dyaic Oupu Feedbac Secod-Order Slidig Mode Coroller for cerai Syses Jeag-Li

More information

Solution. 1 Solutions of Homework 6. Sangchul Lee. April 28, Problem 1.1 [Dur10, Exercise ]

Solution. 1 Solutions of Homework 6. Sangchul Lee. April 28, Problem 1.1 [Dur10, Exercise ] Soluio Sagchul Lee April 28, 28 Soluios of Homework 6 Problem. [Dur, Exercise 2.3.2] Le A be a sequece of idepede eves wih PA < for all. Show ha P A = implies PA i.o. =. Proof. Noice ha = P A c = P A c

More information

The Dynamics of a One-Dimensional Parabolic Problem versus the Dynamics of Its Discretization

The Dynamics of a One-Dimensional Parabolic Problem versus the Dynamics of Its Discretization Joural of Differeial Equaios 168, 6792 (2000) doi:10.1006jdeq.2000.3878, available olie a hp:www.idealibrary.com o The Dyamics of a Oe-Dimesioal Parabolic Problem versus he Dyamics of Is Discreizaio Simoe

More information

Mathematical Statistics. 1 Introduction to the materials to be covered in this course

Mathematical Statistics. 1 Introduction to the materials to be covered in this course Mahemaical Saisics Iroducio o he maerials o be covered i his course. Uivariae & Mulivariae r.v s 2. Borl-Caelli Lemma Large Deviaios. e.g. X,, X are iid r.v s, P ( X + + X where I(A) is a umber depedig

More information

Multifarious Implicit Summation Formulae of Hermite-Based Poly-Daehee Ploynomials

Multifarious Implicit Summation Formulae of Hermite-Based Poly-Daehee Ploynomials Appl. Mah. If. Sci. 12, No. 2, 305-310 (2018 305 Applied Maheaics & Iforaio Scieces A Ieraioal Joural hp://dx.doi.org/10.18576/ais/120204 Mulifarious Iplici Suaio Forulae of Herie-Based Poly-Daehee Ployoials

More information

Problems and Solutions for Section 3.2 (3.15 through 3.25)

Problems and Solutions for Section 3.2 (3.15 through 3.25) 3-7 Problems ad Soluios for Secio 3 35 hrough 35 35 Calculae he respose of a overdamped sigle-degree-of-freedom sysem o a arbirary o-periodic exciaio Soluio: From Equaio 3: x = # F! h "! d! For a overdamped

More information

Oscillator with Distributed Nonlinear Structure on a Segment of Lossy Transmission Line

Oscillator with Distributed Nonlinear Structure on a Segment of Lossy Transmission Line Ope Access Library oural 26, Volue 3, e36 ISSN Olie: 2333-972 ISSN Pri: 2333-975 Oscillaor wih Disribued Noliear Srucure o a Sege of Lossy rasissio Lie Vasil G Agelo Depare of Maheaics, Uiersiy of Miig

More information