A new Approach for Solving Linear Ordinary Differential Equations

Size: px
Start display at page:

Download "A new Approach for Solving Linear Ordinary Differential Equations"

Transcription

1 , ISSN X (Onlne), ISSN (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of Mathematcs, Faculty of Scence, Unversty of Benghaz, Benghaz, Lbya; Emal: fawzabd@hotmal.com ABSTRACT Recently a new numercal scheme, based on seres soluton method, has been ntroduced for solvng Volterra ntegral equatons. In ths wor, we mae use of the natural relaton between ordnary dfferental equatons and Volterra ntegral equatons and apply ths approach to ntal value problems of lnear ordnary dfferental equatons. Ths study led to appromate solutons, whch can be used as algorthms to compute numercal solutons for the ntal value problem. Ths study also ntroduced specal case related to the order of the lnear ordnary dfferental equaton and the ernel of the converted Volterra ntegral equaton. At the end, we nvestgate the sgnfcant features of ths approach, by solvng selected numercal eamples and comparng our results wth the results obtaned by other methods. Under certan condtons, as n the case f ntegral equatons, we found, that ths approach leads to the Taylor epanson of the eact soluton. Keywords: Ordnary dfferental equatons, Integral equatons, Numercal methods. Mathematcs Subject Classfcaton: 65N1, 3C3 1. Introducton Numercal ordnary dfferental equatons studes the numercal soluton of ordnary dfferental equatons [1-]. A varety of methods, analytcal, sem-analytcal, appromate and purely numercal are avalable for the soluton of ordnary dfferental equatons such as the dfferental transformaton method [3], Adoman decomposton method [4] and Runge Kutta methods [5]. Recently a recursve method based on seres soluton method has been ntroduced for solvng Volterra ntegral equatons [6,7]. The man am of ths wor s to apply ths approach to fnd appromate solutons for the ntal value problems of lnear ordnary dfferental equatons. To do that, we use the natural relaton between ordnary dfferental equatons and ntegral equatons [8]. The general form of ntal value problems of lnear ordnary dfferental equatons can be wrtten n the form

2 n y ( n) a ( n ) ( ) ( ) y g, 1 () y () c,,,3,..., n 1 (1.1) where the coeffcents a ( ) are contnuous functons of and c are constants. To convert the ntal value problem (1.1) to Volterra ntegral equatons, we set ( n) y ( ) u ( ), (1.) Net, f we ntegrate the both sdes of (1.) j tmes from to, we can fnd the followng formula j 1 n 1 n j ( n j) ( t) y ( ) u( t) dt c j n j n j, j 1,,..., n (1.3) ( 1)! ( )! Now, by usng (1.) and (1.3), we can convert the ntal value problem (1.1) to the a lnear Volterra ntegral equaton of second nd of the form u ( ) f( ) tutdt (, ) ( ), u() f() (1.4) wth n ( t) t (, ) aj( ) ( j 1)! j 1 j 1 n n 1 j n f( ) g( ) aj( ) c ( j n )! (1.5) j 1 n j. A Numercal Algorthm In ths secton, we follow ref. [6] and apply the recursve method to the ntal value problem (1.1). But for smplcty, let us frst assume that t (, ) (, t),,1,...,. (.1) 1

3 Net, we dfferentate (1.4) tmes wth respect of. Hence for 1,, 3, 4, we have u ( ) f ( ) ( u, ) ( ) ( tutdt, ) ( ) (.) 1 u ( ) f ( ) ( u, ) ( ) ( u, ) ( ) ( u, ) ( ) ( tutdt, ) ( ) (.3) 1 u ( ) f ( ) (, ) u ( ) (, ) u ( ) (, ) u( ) (, ) u( ) 1 (, ) u ( ) (, ) u( ) (, t) u( t) dt 1 3 (.4) u ( ) f ( ) (, ) u ( ) 3 (, ) u ( ) 3 (, ) u ( ) (, ) u( ) (, ) u( ) (, ) u ( ) (, ) u ( ) ( u, ) ( ) ( u, ) ( ) ( u, ) ( ) ( tutdt, ) ( ) 3 4 (.5) To fnd the 1 st appromate soluton of (1.4), we substtute u 1 ( ) E E 1, wth E f() nto (.) and then settng. Ths leads to a formula for the unnown parameter E 1 of the form (1!) E f () (, ) E (.6) 1 Smlarly, to fnd the nd appromate soluton of (1.4), we substtute u ( ) E E E 1 nto (.4) and then settng. Ths leads to a formula for the unnown parameter E of the form (!) E f () (,) E (,) E (,) E (.7) 1 1 Hence, by repeatng the same procedure above, we can fnd the followng formulas for the unnown parameters E 3 and E4 (3!) E f () (,) E (,) E (,) E 3 1 (,) E (,) E (,) E (.8) 13

4 (4!) E f () 6 (,) E 6 (,) E 3 (,) E (,) E (,) E (,) E (,) E (,) E (,) E (,) E 1 3 (.9) Ths leads to 4 th appromate soluton of the ntegral equaton (1.5) of the form u ( ) E E E E E (.11) Note that, the appromate soluton (.1) can be used as an algorthm to fnd numercal soluton for the ntegral equaton (1.5) and hence the ntal value problem (1.1). At the end, f we suppose that coeffcents of the ntal value problem (1.1) and g ( ) are analytcal functons, then t s clear that f ( ) and t (, ) wll be analytc functons wth respect to all ther arguments. In ths case, f u ( ) h ( ) s an eact soluton of the ntegral equaton (1.5), whch we assume t s nfntely dfferentable functon n the neghborhood of, then ths procedure obtans the Taylor epanson of the eact soluton u h( ) around, and hence the Taylor epanson of the eact soluton of ntal value problem (1.1). Ths result can be proven by usng the same procedure of refs. [6,7]. 3. Specal Case In ths secton, we nvestgate ths procedure for frst order lnear ordnary dfferental equaton, or ndependently of the ntal value problem (1.1), for the case when the ernel of the Volterra ntegral equaton (1.4) s a functon of only, say t (, ) H ( ). In ths case, we can rewrte the formulas (.6-9) n the forms (1!) E f () H() E (3.1) 1 (!) E f () H () E H() E (3.) 1 (3!) E f () 3 H () E 3 H () E H() E (3.3)

5 (4!) E f () 4 H () E 6 H () E 8 H () E 6 H() E (3.4) In ths case we can ntroduce a close form formula for the parameters E ; 1,,3,...,. To do that, we compare the coeffcents of the equatons (3.1-4) wth the bnomal coeffcents. Ths comprsng study ntroduces the followng recursve formula for the parameters E f () 1!, 1,,3,... (3.5) ( ) ( ) E H E 1!! 1 ( )( )! Ths leads to the th appromate soluton u ( ) E E (3.6) 1 whch can be used as an algorthm to compute the numercal solutons for the converted ntegral equaton. Furthermore, n ths case, f the coeffcents of the ntal value problem (1.1) and g ( ) are analytcal functons, then (3.6), wth coeffcents gven by the formula (3.5), leads to the Taylor epanson of the eact soluton of converted ntegral equaton. 4. Numercal Applcatons In ths secton, we demonstrate the sgnfcant features of ths method. To do that, we apply ths technque on selected lnear ordnary dfferental equatons and compare our results wth results obtaned by other methods. Eample (1) consder the ntal value problem y y ; y() 3 (4.1) The eact soluton of the ntal value problem (3.1) s y ( ) 1/ (7/) e, whch has the Taylor seres epanson y ( ) 3... (4.)

6 Now (4.1) can be converted to an ntegral equaton of the form (1.4) wth t (, ) and f ( ) 7. Hence by usng the recursve formula (3.5) wth E and H( ), we were able to fnd E1 7, E, E3 7, E4, E5 7/, E6, E7 7/6, E8 and E9 7/4. Ths leads to the followng 9 th appromate soluton u 9( ) 7 7, (4.3) Now by usng (1.) and the ntal condtons, we can wrte the 1 th the ntal value problem (4.1) n the form appromate soluton of ( ) y 1 (4.4) At the end, wth the help of the Mathematca Pacages, we can wrte the soluton of the ntal value problem (4.1) n the form th appromate ( 1) y ( ) (1... ) (4.5) 6 4 1! whch can be wrtten n the form 1 7 ( 1) y ( ) (4.6)! Hence, for ths eample, we can deduce the eact soluton as 1 7 ( 1) 1 7 y ( ) Lm y ( ) e! Eample (): Consder the ntal values problem y 6y 3y 63y, y ( ) y () y () 1 (4.7) 16

7 Now a 1 ( ) 6, a ( ) 3, ( ) 63 3 a, c c1 c 1, g ( ) and n 3. Usng ths data and (1.4) and (1.5), we can convert the ntal value problem (4.7) to an ntegral equaton of the form u( ) ( t) 315( t) u( t) dt (4.8) Hence, f () 433, f () 47, f () 63, () f () ; (,) 6, 1 (,) 33, (,) 63 and (,) ; ( 3,4,...). To fnd th appromate soluton of the ntegral equaton (4.8), we used ths data wth E 433. Now, the formulas (.6-9) leads to the 4 th appromate soluton u ( ) (4.9) In addton, f mae use of (1.) and the ntal condtons, we can wrte the 7 th soluton of the ntal value problem (4.7) n the form appromate y ( ) (4.1) 7 54 At the end, we can use (4.1) as an algorthm to fnd a numercal soluton of the ntal value problem (4.7). Furthermore, the graphs (4.1) and (4.), show the plottng of the eact soluton (red) ( y( ) 1e 7e 16e ), aganst the appromate solutons y 3 ( ) and y ( ) 4 respectvely. 17

8 Graph (4.1) Plottng the eact soluton aganst y 3 ( ) Graph (4.) Plottng the eact soluton aganst y 4 ( ) Eample (3): Consder the ntal values problem y 3y y 3, y () = 1, y () = (4.11) In ref.[9], ths ntal value problem has been solved by usng the dfferental transformaton method. Ths soluton was gven n term of the power seres y ( ) 1... (4.1)! 3! 4! 5! Now by usng a 1 ( ) 3, a ( ), c 1 c, g ( ) 3 convert the ntal value problem (4.11) to an ntegral equaton of the form and n, we can 18

9 u ( ) 1 3 ( t) utdt ( ) (4.13) Usng the same prevous proceedng we can fnd f () 1, f (), () f (), (,) 3, (,) 1, and (,) ; (,3,...). Hence, the formulas (.6-9) wth E 1, gves E1 1, E 1/, E3 1/6 and E4 1/4. Ths leads to the 4 th appromate soluton u ( ) 1 (4.14) At the end, f mae use of (1.) and the ntal condtons, we can wrte the 6 th soluton of the ntal value problem (4.11) n the form appromate y ( ) 1 (4.15) Wth the help of the Mathematca Pacages, ths proceedng enables us to wrte the appromate soluton of the ntal value problem (4.11) n the form th y ( ) 1...! 3! 4!!! (4.16) Hence, t s clear that the eact soluton s gven by 1 y( ) Lm y ( ) e (4.17)! Eample (4): Consder the ntal value problem y y, (4.18) y() d =, 1 y () = d (4.19) 19

10 The equaton (4.18), whch s nown as the Ary equaton, has no soluton n term of elementary functons, as many other mportant ordnary dfferental equatons. The most mportant soluton of ths equaton s the Ary functon of the frst nd A(, ) whch ts propertes have been etensvely studed [1,11]. Adoman [4] appled the decomposton method technque to fnd the seres soluton of the Ary equaton (4.18), however, after hard wor, the seres soluton method gves the seres soluton n the form æ 3 ö ç å = 1 ( 3)(5 6) ((3-1) (3 )) y ( ) = d 1+ + è ø æ 3+ 1 ö 1 + ç å = 1 (3 4)(6 7) ((3 ) (3+ 1)) d è ø (4.) Now a 1 ( ), a ( ), c d c d1, g ( ) and n. Usng ths data, we can convert the ntal value problem ( ) to an ntegral equaton of the form u ( ) d d ( tutdt ) ( ) 1 (4.1) As n the prevous eample, we can calculate f (), f () d, f () d1 () f (), (, ), (, ), 1 (, ) and (,) ; ( 3,4,...). Hence the formulas (.6-9) wth E, gves E 1 d, E d1, E3 and E d /6 4. Ths leads to the 4th appromate soluton d 4 u ( ) d d (4.) Now, f mae use of (1.) and the ntal condtons, we can wrte the 6 th soluton n the form appromate d 3 d 1 4 d 6 y ( ) d d (4.3) Ths can be wrtten n the form

11 y ( ) d 1 6 d ( 3) ( 3)(5 6) 1 (4.4) (3 4) whch agrees wth the seres soluton (4.). Note that, ths study showed the accuracy, the relablty and the effcency of ths method comparng wth dfferental transformaton method, the Adoman method and the seres soluton method. 5. Concluson In ths wor, we studed the soluton of the ntal value problems of ordnary dfferental equatons by applyng a recursve method technque. Ths study led to appromate solutons, whch can be used as algorthms to compute numercal solutons for the ntal value problem. In addton, ths study enables us to ntroduce a formula for specal case related to the order of the ordnary dfferental equatons and the ernel of the converted ntegral equaton. To demonstrate the sgnfcant features of ths approach, we appled ths technque on selected lnear ordnary dfferental equatons and compare our results wth the eact soluton and the results obtaned by dfferent technques such as the dfferental transformaton method the Adoman method and the seres soluton method. Ths comparson study shows the accuracy, the relablty and the effcency of ths method. Certanly much remans to be done on ntroducng the same procedure for other types of dfferental equatons, whch can be vewed as a fascnatng challenge for further research. 6. References [1] J. C. Butcher, Numercal Methods for Ordnary Dfferental Systems, John Wley & Sons, 3. [] J. C. Butcher, Numercal methods for ordnary dfferental equatons n the th century, Journal of Computatonal and Appled Mathematcs 15 () 1-9. [3] S. Chang and I. Chang, A new algorthm for calculatng one-dmensonal dfferental transform of nonlnear functons, Appled Mathematcs and Computaton 195 (8) [4] G. Adoman, A revew of the decomposng method and some recent results for non-lnear equatons, Math. Comput. Model. 13(7) (199) [5] D. Zngg, T. Chsholmm, Runge Kutta methods for lnear ordnary dfferental equatons, Appled Numercal Mathematcs 31 (1999)

12 [6] A. Tahmasb, O. S. Fard, Numercal Soluton of Nonlnear Volterra Integral Equaton of the Second Knd by Power Seres, Journal of Informaton and Computng Scence, Vol. 3, No. 1 (8) [7] A. Tahmasb, O. S. Fard, A Numercal Scheme to solve Nonlnear Volterra Integral Equatons System of the Second Knd, Journal of Informaton and Computng Scence, Vol. 3, No. (8) [8] A. M. Wazwaz, A Frst Course n Integral Equatons, WSPC, New Jersey, (1997). [9] K. Batha and B. Batha, A New Algorthm for Solvng Lnear Ordnary Dfferental Equatons, World Appled Scences Journal 15 (1): , 11. [1] M. Abramowtz and I. Stegun, Handboo of mathematcal functons, Dover boos n mathematcs, [11] Danel Zwllnger, Handboo of dfferental equatons, 3rd edton, Academc Press, 1997.

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

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

More information

A Hybrid Variational Iteration Method for Blasius Equation

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

More information

arxiv: v1 [math.co] 12 Sep 2014

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

More information

The binomial transforms of the generalized (s, t )-Jacobsthal matrix sequence

The binomial transforms of the generalized (s, t )-Jacobsthal matrix sequence Int. J. Adv. Appl. Math. and Mech. 6(3 (2019 14 20 (ISSN: 2347-2529 Journal homepage: www.jaamm.com IJAAMM Internatonal Journal of Advances n Appled Mathematcs and Mechancs The bnomal transforms of the

More information

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

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

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

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

More information

The Analytical Solution of a System of Nonlinear Differential Equations

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

More information

Convexity preserving interpolation by splines of arbitrary degree

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

More information

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

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

More information

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method Soluton of Lnear System of Equatons and Matr Inverson Gauss Sedel Iteraton Method It s another well-known teratve method for solvng a system of lnear equatons of the form a + a22 + + ann = b a2 + a222

More information

Beyond Zudilin s Conjectured q-analog of Schmidt s problem

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

More information

VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES

VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES BÂRZĂ, Slvu Faculty of Mathematcs-Informatcs Spru Haret Unversty barza_slvu@yahoo.com Abstract Ths paper wants to contnue

More information

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

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

More information

Erbakan University, Konya, Turkey. b Department of Mathematics, Akdeniz University, Antalya, Turkey. Published online: 28 Nov 2013.

Erbakan University, Konya, Turkey. b Department of Mathematics, Akdeniz University, Antalya, Turkey. Published online: 28 Nov 2013. Ths artcle was downloaded by: [Necmettn Erbaan Unversty] On: 24 March 2015, At: 05:44 Publsher: Taylor & Francs Informa Ltd Regstered n England and Wales Regstered Number: 1072954 Regstered offce: Mortmer

More information

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

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

More information

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

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

More information

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS

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

More information

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

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

More information

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

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

More information

Binomial transforms of the modified k-fibonacci-like sequence

Binomial transforms of the modified k-fibonacci-like sequence Internatonal Journal of Mathematcs and Computer Scence, 14(2019, no. 1, 47 59 M CS Bnomal transforms of the modfed k-fbonacc-lke sequence Youngwoo Kwon Department of mathematcs Korea Unversty Seoul, Republc

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

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

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

More information

Finding Dense Subgraphs in G(n, 1/2)

Finding Dense Subgraphs in G(n, 1/2) Fndng Dense Subgraphs n Gn, 1/ Atsh Das Sarma 1, Amt Deshpande, and Rav Kannan 1 Georga Insttute of Technology,atsh@cc.gatech.edu Mcrosoft Research-Bangalore,amtdesh,annan@mcrosoft.com Abstract. Fndng

More information

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

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

More information

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

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

More information

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

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

More information

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

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

More information

Numerical Solution of Ordinary Differential Equations

Numerical Solution of Ordinary Differential Equations Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples

More information

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

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

More information

NUMERICAL DIFFERENTIATION

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

More information

Solution for singularly perturbed problems via cubic spline in tension

Solution for singularly perturbed problems via cubic spline in tension ISSN 76-769 England UK Journal of Informaton and Computng Scence Vol. No. 06 pp.6-69 Soluton for sngularly perturbed problems va cubc splne n tenson K. Aruna A. S. V. Rav Kant Flud Dynamcs Dvson Scool

More information

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

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

More information

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

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

More information

Inexact Newton Methods for Inverse Eigenvalue Problems

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

More information

Numerical Heat and Mass Transfer

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

More information

Appendix B. The Finite Difference Scheme

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

More information

A New Recursive Method for Solving State Equations Using Taylor Series

A New Recursive Method for Solving State Equations Using Taylor Series I J E E E C Internatonal Journal of Electrcal, Electroncs ISSN No. (Onlne) : 77-66 and Computer Engneerng 1(): -7(01) Specal Edton for Best Papers of Mcael Faraday IET Inda Summt-01, MFIIS-1 A New Recursve

More information

Remarks on the Properties of a Quasi-Fibonacci-like Polynomial Sequence

Remarks on the Properties of a Quasi-Fibonacci-like Polynomial Sequence Remarks on the Propertes of a Quas-Fbonacc-lke Polynomal Sequence Brce Merwne LIU Brooklyn Ilan Wenschelbaum Wesleyan Unversty Abstract Consder the Quas-Fbonacc-lke Polynomal Sequence gven by F 0 = 1,

More information

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

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

More information

Difference Equations

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

More information

Bezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0

Bezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0 Bezer curves Mchael S. Floater August 25, 211 These notes provde an ntroducton to Bezer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of the

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

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

More information

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

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

More information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

Haar wavelet collocation method to solve problems arising in induction motor

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

More information

Lecture 12: Discrete Laplacian

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

More information

Computers and Mathematics with Applications

Computers and Mathematics with Applications Computers and Mathematcs wth Applcatons 56 (2008 3204 3220 Contents lsts avalable at ScenceDrect Computers and Mathematcs wth Applcatons journal homepage: www.elsever.com/locate/camwa An nnovatve egenvalue

More information

Chapter 3 Differentiation and Integration

Chapter 3 Differentiation and Integration MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton

More information

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

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

More information

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence

More information

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse

More information

Grid Generation around a Cylinder by Complex Potential Functions

Grid Generation around a Cylinder by Complex Potential Functions Research Journal of Appled Scences, Engneerng and Technolog 4(): 53-535, 0 ISSN: 040-7467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around

More information

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

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

More information

6.3.4 Modified Euler s method of integration

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

More information

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

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

More information

International Conference on Advanced Computer Science and Electronics Information (ICACSEI 2013) equation. E. M. E. Zayed and S. A.

International Conference on Advanced Computer Science and Electronics Information (ICACSEI 2013) equation. E. M. E. Zayed and S. A. Internatonal Conference on Advanced Computer Scence and Electroncs Informaton (ICACSEI ) The two varable (G'/G/G) -expanson method for fndng exact travelng wave solutons of the (+) dmensonal nonlnear potental

More information

Modelli Clamfim Equazioni differenziali 7 ottobre 2013

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

More information

Comparison Solutions Between Lie Group Method and Numerical Solution of (RK4) for Riccati Differential Equation

Comparison Solutions Between Lie Group Method and Numerical Solution of (RK4) for Riccati Differential Equation Appled and Computatonal Mathematcs 06; 5(: 64-7 http://www.scencepublshnggroup.com/j/acm do: 0.648/j.acm.06050.5 ISSN: 8-5605 (Prnt; ISSN: 8-56 (Onlne Comparson Solutons Between Le Group Method and Numercal

More information

New Method for Solving Poisson Equation. on Irregular Domains

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

More information

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0 Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of

More information

Significance of Dirichlet Series Solution for a Boundary Value Problem

Significance of Dirichlet Series Solution for a Boundary Value Problem IOSR Journal of Engneerng (IOSRJEN) ISSN (e): 5-3 ISSN (p): 78-879 Vol. 6 Issue 6(June. 6) V PP 8-6 www.osrjen.org Sgnfcance of Drchlet Seres Soluton for a Boundary Value Problem Achala L. Nargund* and

More information

Improvement in Estimating the Population Mean Using Exponential Estimator in Simple Random Sampling

Improvement in Estimating the Population Mean Using Exponential Estimator in Simple Random Sampling Bulletn of Statstcs & Economcs Autumn 009; Volume 3; Number A09; Bull. Stat. Econ. ISSN 0973-70; Copyrght 009 by BSE CESER Improvement n Estmatng the Populaton Mean Usng Eponental Estmator n Smple Random

More information

2-π STRUCTURES ASSOCIATED TO THE LAGRANGIAN MECHANICAL SYSTEMS UDC 531.3: (045)=111. Victor Blãnuţã, Manuela Gîrţu

2-π STRUCTURES ASSOCIATED TO THE LAGRANGIAN MECHANICAL SYSTEMS UDC 531.3: (045)=111. Victor Blãnuţã, Manuela Gîrţu FACTA UNIVERSITATIS Seres: Mechancs Automatc Control and Robotcs Vol. 6 N o 1 007 pp. 89-95 -π STRUCTURES ASSOCIATED TO THE LAGRANGIAN MECHANICAL SYSTEMS UDC 531.3:53.511(045)=111 Vctor Blãnuţã Manuela

More information

LECTURE 9 CANONICAL CORRELATION ANALYSIS

LECTURE 9 CANONICAL CORRELATION ANALYSIS LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of

More information

Randić Energy and Randić Estrada Index of a Graph

Randić Energy and Randić Estrada Index of a Graph EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 5, No., 202, 88-96 ISSN 307-5543 www.ejpam.com SPECIAL ISSUE FOR THE INTERNATIONAL CONFERENCE ON APPLIED ANALYSIS AND ALGEBRA 29 JUNE -02JULY 20, ISTANBUL

More information

Trees and Order Conditions

Trees and Order Conditions Trees and Order Condtons Constructon of Runge-Kutta order condtons usng Butcher trees and seres. Paul Tranqull 1 1 Computatonal Scence Laboratory CSL) Department of Computer Scence Vrgna Tech. Trees and

More information

ME 501A Seminar in Engineering Analysis Page 1

ME 501A Seminar in Engineering Analysis Page 1 umercal Solutons of oundary-value Problems n Os ovember 7, 7 umercal Solutons of oundary- Value Problems n Os Larry aretto Mechancal ngneerng 5 Semnar n ngneerng nalyss ovember 7, 7 Outlne Revew stff equaton

More information

2 Finite difference basics

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

More information

Formal solvers of the RT equation

Formal solvers of the RT equation Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan

More information

1 Introduction We consider a class of singularly perturbed two point singular boundary value problems of the form: k x with boundary conditions

1 Introduction We consider a class of singularly perturbed two point singular boundary value problems of the form: k x with boundary conditions Lakshm Sreesha Ch. Non Standard Fnte Dfference Method for Sngularly Perturbed Sngular wo Pont Boundary Value Problem usng Non Polynomal Splne LAKSHMI SIREESHA CH Department of Mathematcs Unversty College

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve

More information

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

More information

A FORMULA FOR COMPUTING INTEGER POWERS FOR ONE TYPE OF TRIDIAGONAL MATRIX

A FORMULA FOR COMPUTING INTEGER POWERS FOR ONE TYPE OF TRIDIAGONAL MATRIX Hacettepe Journal of Mathematcs and Statstcs Volume 393 0 35 33 FORMUL FOR COMPUTING INTEGER POWERS FOR ONE TYPE OF TRIDIGONL MTRIX H Kıyak I Gürses F Yılmaz and D Bozkurt Receved :08 :009 : ccepted 5

More information

Global Sensitivity. Tuesday 20 th February, 2018

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

More information

Section 3.6 Complex Zeros

Section 3.6 Complex Zeros 04 Chapter Secton 6 Comple Zeros When fndng the zeros of polynomals, at some pont you're faced wth the problem Whle there are clearly no real numbers that are solutons to ths equaton, leavng thngs there

More information

1 GSW Iterative Techniques for y = Ax

1 GSW Iterative Techniques for y = Ax 1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn

More information

A Computational Quadruple Laplace Transform for the Solution of Partial Differential Equations

A Computational Quadruple Laplace Transform for the Solution of Partial Differential Equations Appled Mathematcs,, 5, 337-338 Publshed Onlne December n ScRes http://wwwscrporg/journal/am http://dxdoorg/36/am53 A Computatonal Quadruple Laplace Transform for the Soluton of Partal Dfferental Equatons

More information

Hidden Markov Models

Hidden Markov Models Hdden Markov Models Namrata Vaswan, Iowa State Unversty Aprl 24, 204 Hdden Markov Model Defntons and Examples Defntons:. A hdden Markov model (HMM) refers to a set of hdden states X 0, X,..., X t,...,

More information

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function A Local Varatonal Problem of Second Order for a Class of Optmal Control Problems wth Nonsmooth Objectve Functon Alexander P. Afanasev Insttute for Informaton Transmsson Problems, Russan Academy of Scences,

More information

Integral Transforms and Dual Integral Equations to Solve Heat Equation with Mixed Conditions

Integral Transforms and Dual Integral Equations to Solve Heat Equation with Mixed Conditions Int J Open Probles Copt Math, Vol 7, No 4, Deceber 214 ISSN 1998-6262; Copyrght ICSS Publcaton, 214 www-csrsorg Integral Transfors and Dual Integral Equatons to Solve Heat Equaton wth Mxed Condtons Naser

More information

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

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

More information

Ballot Paths Avoiding Depth Zero Patterns

Ballot Paths Avoiding Depth Zero Patterns Ballot Paths Avodng Depth Zero Patterns Henrch Nederhausen and Shaun Sullvan Florda Atlantc Unversty, Boca Raton, Florda nederha@fauedu, ssull21@fauedu 1 Introducton In a paper by Sapounaks, Tasoulas,

More information

THE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD

THE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD Journal of Appled Mathematcs and Computatonal Mechancs 06, 5(), 7-36 www.amcm.pcz.pl p-iss 99-9965 DOI: 0.75/jamcm.06..4 e-iss 353-0588 THE STURM-LIOUVILLE EIGEVALUE PROBLEM - A UMERICAL SOLUTIO USIG THE

More information

Projective change between two Special (α, β)- Finsler Metrics

Projective change between two Special (α, β)- Finsler Metrics Internatonal Journal of Trend n Research and Development, Volume 2(6), ISSN 2394-9333 www.jtrd.com Projectve change between two Specal (, β)- Fnsler Metrcs Gayathr.K 1 and Narasmhamurthy.S.K 2 1 Assstant

More information

ORDINARY DIFFERENTIAL EQUATIONS EULER S METHOD

ORDINARY DIFFERENTIAL EQUATIONS EULER S METHOD Numercal Analss or Engneers German Jordanan Unverst ORDINARY DIFFERENTIAL EQUATIONS We wll eplore several metods o solvng rst order ordnar derental equatons (ODEs and we wll sow ow tese metods can be appled

More information

SHIFTED JACOBI COLLOCATION METHOD BASED ON OPERATIONAL MATRIX FOR SOLVING THE SYSTEMS OF FREDHOLM AND VOLTERRA INTEGRAL EQUATIONS

SHIFTED JACOBI COLLOCATION METHOD BASED ON OPERATIONAL MATRIX FOR SOLVING THE SYSTEMS OF FREDHOLM AND VOLTERRA INTEGRAL EQUATIONS Mathematcal and Computatonal Applcatons, Vol., o., pp. 76-93, 5 http://d.do.org/.9/mca-5-7 SHIFED JACOBI COLLOCAIO MEHOD BASED O OPERAIOAL MARIX FOR SOLVIG HE SYSEMS OF FREDHOLM AD VOLERRA IEGRAL EQUAIOS

More information

Fundamental loop-current method using virtual voltage sources technique for special cases

Fundamental loop-current method using virtual voltage sources technique for special cases Fundamental loop-current method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

Maejo International Journal of Science and Technology

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

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

On a direct solver for linear least squares problems

On a direct solver for linear least squares problems ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear

More information

Integrals and Invariants of Euler-Lagrange Equations

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

More information

Modelli Clamfim Equazioni differenziali 22 settembre 2016

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

More information

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

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

More information

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George

More information

The Study of Teaching-learning-based Optimization Algorithm

The Study of Teaching-learning-based Optimization Algorithm Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute

More information

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

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

More information

MATH 829: Introduction to Data Mining and Analysis The EM algorithm (part 2)

MATH 829: Introduction to Data Mining and Analysis The EM algorithm (part 2) 1/16 MATH 829: Introducton to Data Mnng and Analyss The EM algorthm (part 2) Domnque Gullot Departments of Mathematcal Scences Unversty of Delaware Aprl 20, 2016 Recall 2/16 We are gven ndependent observatons

More information

Solution of Singularly Perturbed Differential Difference Equations Using Higher Order Finite Differences

Solution of Singularly Perturbed Differential Difference Equations Using Higher Order Finite Differences Amercan Journal of Numercal Analyss, 05, Vol. 3, No., 8-7 Avalable onlne at http://pubs.scepub.com/ajna/3// Scence and Educaton Publshng DOI:0.69/ajna-3-- Soluton of Sngularly Perturbed Dfferental Dfference

More information