THE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD
|
|
- Earl Parks
- 6 years ago
- Views:
Transcription
1 Journal of Appled Mathematcs and Computatonal Mechancs 06, 5(), p-iss DOI: 0.75/jamcm e-iss THE STURM-LIOUVILLE EIGEVALUE PROBLEM - A UMERICAL SOLUTIO USIG THE COTROL VOLUME METHOD Jarosław Sedleck, Marusz Ceselsk, Tomasz Błaszczyk Insttute of Mathematcs, Czestochowa Unversty of Technology Częstochowa, Poland Insttute of Computer and Informaton Scences Czestochowa Unversty of Technology, Częstochowa, Poland jaroslaw.sedleck@m.pcz.pl, marusz.ceselsk@cs.pcz.pl, tomasz.blaszczyk@m.pcz.pl Abstract. The soluton of the D Sturm-Louvlle problem usng the Control Volume Method s dscussed. The second order lnear dfferental equaton wth homogeneous boundary condtons s dscretzed and converted to the system of lnear algebrac equatons. The matrx assocated wth ths system s trdagonal and egenvalues of ths system are an approxmaton of the real egenvalues of the boundary value problem. The numercal results of the egenvalues for varous cases and the expermental rate of convergence are presented. Keywords: Sturm-Louvlle problem, egenvalues, numercal methods, Control Volume Method. Introducton Ths paper s concerned wth the computaton of egenvalues of regular egenvalue problems occurrng n ordnary dfferental equatons. The Sturm-Louvlle problem arses wthn n many areas of scence, engneerng and appled mathematcs. It has been studed for more than two decades. Many physcal, bologcal and chemcal processes are descrbed usng models based on the Sturm-Louvlle equatons. The Sturm-Louvlle problem appears drectly as the egenvalue problem n a one-dmensonal space. It also arses when lnear partal dfferental equatons are separable n a certan coordnate system. For a more detaled study of the nteger order Sturm-Louvlle theory, we refer the reader to [-5]. The Sturm-Louvlle problem can be solved by usng ether analytcal or numercal methods. One of the most common approaches to a numercal soluton of the consdered problem s the fnte dfference method [3, 6-8] where each dervatve s dscretzed at each grd pont wth an adequate dfference scheme. Apart from the
2 8 J. Sedleck, M. Ceselsk, T. Błaszczyk fnte dfference method, several analytcal ones, such as the varatonal or decomposton methods, are proposed to fnd an approxmate soluton.. Statement of the problem We consder the followng problem defned on the bounded nterval x [a, b] d dy x p( x) + q x y x = w x y x dx dx α y a + α y ' a = 0 β y b + β y ' b = 0 () () where p(x) > 0, dp/dx and q(x) are contnuous, w(x) > 0 on [a, b], α + α 0 and β + β 0. The above problem s called the regular Sturm-Louvlle Problem (SLP). The soluton to SLP conssts of a par and y, where s a constant - called an egenvalue, whle y s a nontrval (nonzero) functon - called an egenfuncton, and together they satsfy the gven SLP. For each SLP, the egenvalues form an nfnte ncreasng sequence: < < 3 < and lm k k =. For arbtrary choces of the functons p(x), q(x) and w(x) n Eq. (), the computaton of the exact values of the egenvalues for whch SLP ()-() has a nontrval egensoluton y(x) whch s very complcated or t s practcally mpossble to determne. umercal methods should, therefore, be used for computng the approxmate values of. In many cases, the mportance of a numercal approxmaton to the SLP descrbed by a dfferental egenvalue problem s to reduce the problem to that of solvng the egenvalue problem of a matrx equaton (an algebrac problem). In ths paper, we apply the Control Volume Method (also known as the Fnte Volume Method) to compute the egenvalues of the Sturm-Louvlle problem numercally. 3. Control Volume Method In ths numercal method [9] the consdered doman of SLP: x [a, b] s dvded nto control volumes Ω for =,, wth the central nodes ξ. The mesh s presented n Fgure. The auxlary nodes x = a + x, for = 0,, and x = = (b a) / have also been ntroduced. Then, the postons of central nodes are the followng: ξ = a + ( 0.5) x, =,,.
3 The Sturm-Louvlle egenvalue problem - a numercal soluton usng the Control Volume Method 9 Fg.. Mesh of control volumes The ntegraton of Eq. () wth respect to volume Ω leads to d dy x p( x) dv + q( x) y( x) dv = w( x) y( x) dv dx dx (3) Ω Ω Ω or wrtten n the form (assumng that Ω : [x -, x ]) x x x d dy x p( x) dx + q( x) y( x) dx = w( x) y( x) dx dx dx (4) x x x All the components n Eq. (4) can be approxmated as follows: x x q( x) y( x) dx q( ξ) y( ξ) x (5) x x w( x) y( x) dx w( ξ) y( ξ) x (6) x d dy x dy x dy x p( x) dx = p( x) + p( x) dx dx dx dx x x= x x= x y ξ y ξ y ξ+ y ξ for > for < (7) x x p( x ) + p( x) y( ξ) ya yb y( ξ ) for = for = 0.5 x 0.5 x The values of y a = y(a) and y b = y(b) are determned on the bass of approxmatons of the boundary condtons () ( ξ ) y y( ξ ) y ya b α ya + α = 0, β yb + β = x 0.5 x and hence, these values are equal to (8) α β y = y ξ y = y ξ, a α α b x β + β x (9)
4 30 J. Sedleck, M. Ceselsk, T. Błaszczyk Substtutng (9) nto (8), the followng approxmaton of (7) has been obtaned where x x p( x ) ( ) ( ) dx y d dy x y ξ y ξ, for > p( x) dx dx x γ a ξ, for = 4α ( ) ( + ) ( ) x p x y ξ y ξ, for < + γ b y ξ, for = γ a =, b α α x (0) 4β γ = () β + β x For example, n the case of the Drchlet boundary condtons at boundares x = a (.e. α =, α = 0) and/or x = b (.e. β =, β = 0), the coeffcents γ a and γ b are equal to, and n the case of the eumann boundary condtons at x = a (α = 0, α = ) and/or x = b (β = 0, β = ), these coeffcents take the values of 0, respectvely. After substtuton of (5), (6) and (0) nto (4), the followng system of the dscrete equatons for every control volume: for Ω : ( x) ( x) γ p x a + p x + q( ξ ) y( ξ ) p x y ξ = w ξ y ξ 0 () for Ω, for =,, : ( ) = w( ξ ) y( ξ ) ( ) + p x p x p x p x y ξ + + q( ξ ) y( ξ ) y ξ x x x + (3) for Ω : ( ) ( x) ( ) ( x) p x γ p x + p x y( ξ ) + + q ξ y ξ = w ξ y ξ b (4) s obtaned. Ths system can be also wrtten n the matrx form as ( P+ Q) y=wy (5)
5 where P The Sturm-Louvlle egenvalue problem - a numercal soluton usng the Control Volume Method 3 ( x) ( 0) γa p x p( x) p x p( x) p( x) p( x) p( x) p( x ) 0 p( x ) p( x) 0 p( x) (6) p( x ) 0 0 p( x ) p( x ) + p( x ) p( x ) p( x ) +γb p( x) = + and y y, y,..., y ( q( ) q q) Q = dag ξ, ξ,..., ξ (7) ( w( ) w w) W = dag ξ, ξ,..., ξ (8) T = ξ ξ ξ. Then, SLP ()-() s equvalent to the matrx egenvalue problem (5). The ordered egenvalues of (5) are denoted by, =,,. In order to evaluate egenvalues of large matrx egenvalue problem (5), the numercal methods mplemented n mathematcal software can be used. 3.. Partcular case Assumng the functons p(x) =, q(x) = 0 and w(x) = n the SLP problem, then matrces P, Q and W are reduced to the followng forms: Q = O (Zero matrx), W = I (Identty matrx) and γ a P = ( x) γ b + (9)
6 3 J. Sedleck, M. Ceselsk, T. Błaszczyk In the case of the mxed Drchlet and/or eumann boundary condtons at both boundares, one can fnd n lterature (e.g. [4]) the explct formulas for the -th egenvalue of the SLP ()-() and they are presented n Table. For the dscrete case, the egenvalues of the matrx P can be determned n an analytcal way. The egenvalue problems of trdagonal matrces (n a smlar form as the matrx P) are consdered n [0, ]. On the bass of these results, the egenvalues of the matrx P:, =,, are adopted to the analysed problem for four mxed boundary condtons and they are also gven n Table. The Drchlet B.C. at x = a γ a = (for α =, α = 0) The eumann B.C. at x = a γ a = 0 (for α = 0, α = ) Egenvalues of Eq. () for functons p(x) =, q(x) = 0 and w(x) = and four mxed cases of boundary condtons (explct formulas for, =,, and dscrete one for, =,, ) The Drchlet B.C. at x = b γ b = (for β =, β = 0) π = b a ( ) 4 π = sn b a = ( b a) ( 0.5) π b a ( 0.5) 4 π ( ) = sn The eumann B.C. at x = b γ b = 0 (for β = 0, β = ) = ( b a) ( 0.5) π b a ( 0.5) 4 π ( ) = sn ( ) π = b a ( ) 4 ( ) π = sn ( b a) Table For more complcated cases such as non-constant functons p(x), q(x), w(x) occurrng n Eq. (), the egenvalues of system (5) should be determned n a numercal way (.e. usng mathematcal software - here the Maple s used). 3.. Error of numercal approxmaton of the egenvalues Theorem. For the Sturm-Louvlle problem ()-() there exsts a constant C such that the -th exact egenvalue and the approxmate egenvalue satsfy the followng relaton 4 4 C x = O x, x = b a /, =,..., (0) Proof: For the presented formulas n Table (n these cases, the exact and dscrete egenvalues are known), one can estmate the error of approxmaton of the dscrete egenvalues. Let us start from the Taylor seres expanson for the functon sn (x) at about pont x = 0
7 The Sturm-Louvlle egenvalue problem - a numercal soluton usng the Control Volume Method 33 sn ( x) n= ( ) ( n) n+ n n x = ()! ext, f we take nto account the frst two terms n the Taylor seres (), then we estmate the followng expresson as ( n) n+ n n x 4 4 x sn ( x) = x < x x x = x n=! 3 3 () For the case of the Drchlet-Drchlet boundary condtons, the error s evaluated by usng the estmaton () n the followng way: ( ) ( ) ( ) 4 sn sn ( ) π π π π = = b a b a b a ( ) π π π ( b a) ( b a) ( b a) 4 4 < = = x = C x 4 3 (3) In the other cases (from Table ), the results are smlar. Another method of error estmaton s based on the nvestgaton of the Expermental Rate of Convergence (ERC). The total error n the estmate of the egenvalues s composed of both the error resultng from the dscretzaton of the equaton and the error of the numercal algorthm for fndng egenvalues of system (5). Here, we assume that the error s r s ( ) = = = (4) ( ) O x, x b a /,,..., where the parameters r and s are to be determned expermentally. If the -th exact egenvalue to SLP s known, then the parameters r and s can be determned usng the followng formulas for the ERC for varable values of : r= ERCr(, ) = log ( /) (5) ( ) s(, ) log ( ) s= ERC = Whereas, f the exact egenvalue s unknown then, we determne the parameter r from the followng formula (6) ( ) ( /) r(, ) log ( ) ( ) r= ERC = (7)
8 34 J. Sedleck, M. Ceselsk, T. Błaszczyk The parameter s can be estmated usng (6) and assumng that s numercally determned for suffcently hgh value of. 4. Example of numercal smulatons In tests of verfcaton of the numercal solutons, three cases are taken nto account: Example : p(x) =, q(x) = 0, w(x) = ( + x), α =, α = 0, β =, β = 0. Example : p(x) = (x + ), q(x) = x, w(x) = exp(x), α =, α = 0, β = 0, β =. Example 3: p(x) = + sn(πx), q(x) = 0, w(x) = + sqrt(x), α =, α = 0, β = 5, β =. In all examples, the values of a = 0 and b = have been assumed. In the case of Example, the exact egenvalues are gven by ( ) = / 4 + π / ln [], whle for the remanng cases the exact egenvalues are unknown. In Tables, 4 and 5, the numercal values of the frst 8 egenvalues for dfferent values of and the calculated ((5) or (7)) for all examples are presented, respectvely. In addton, n Table 3, the calculated values of ERC s (6) for Example are shown. Egenvalues and for Example Table anal anal
9 The Sturm-Louvlle egenvalue problem - a numercal soluton usng the Control Volume Method 35 ERC s for Example Table 3 ERC s (,) ERC s (,) ERC s (,3) ERC s (,4) ERC s (,5) ERC s (,6) ERC s (,7) ERC s (,8) Egenvalues and for Example Table Egenvalues and for Example 3 4 Table The analyss of the results presented n Tables -5 ndcates that the rate r ( ) s close to, whle the rate s (ERC s ) s close to 4. Thus we can confrm that the relatonshp (0) s satsfed. The errors n the approxmaton of egenvalues ncrease rapdly as the ndex of the egenvalue grows. 3 4
10 36 J. Sedleck, M. Ceselsk, T. Błaszczyk 5. Conclusons In ths paper, the new approach based on the control volume method for fndng the egenvalues of the Sturm-Louvlle problem was dscussed. The contnuous problem descrbed by the dfferental equaton wth the adequate boundary condtons was converted to the correspondng dscrete one. The rate of convergence of the proposed numercal scheme s order. The presented results of the approxmaton of egenvalues are n close agreement wth the results obtaned n an analytcal way or n the mathematcal software feld. In the future, the presented approach can be extended to apply hgh order of accuracy dfference schemes for approxmatng egenvalues of the Sturm-Louvlle problem and can be appled to the fractonal Sturm-Louvlle problem whch s related to the correspondng fractonal Euler- -Lagrangan equaton [3]. References [] Agarwal R.P., O Regan D., An Introducton to Ordnary Dfferental Equatons, Sprnger, ew York 008. [] Atknson F.V., Dscrete and Contnuous Boundary Value Problems, Academc Press, ew York, London 964. [3] Pryce J.D., umercal Soluton of Sturm-Louvlle Problems, Oxford Unv. Press, London 993. [4] Zatsev V.F., Polyann A.D., Handbook of Exact Solutons for Ordnary Dfferental Equatons, CRC Press, ew York 995. [5] Pruess S., Estmatng the egenvalues of Sturm-Louvlle problems by approxmatng the dfferental equaton, SIAM J. umer. Anal. 973, 0, [6] Aceto L., Ghelardon P., Maghern C., Boundary value methods as an extenson of umerov s method for Sturm-Louvlle egenvalue estmates, Appled umercal Mathematcs 009, 59 (7), [7] Amodo P., Settann G., A matrx method for the soluton of Sturm-Louvlle problems, Journal of umercal Analyss, Industral and Appled Mathematcs (JAIAM) 0, 6 (-), -3. [8] Ascher U.M., umercal Methods for Evolutonary Dfferental Equatons, SIAM, 008. [9] Ascher U.M., Matthej R.M.M., Russell R.D., umercal Soluton of Boundary Value Problems for ODEs, Classcs n Appled Mathematcs 3, SIAM, Phladelpha 995. [0] Ellott J.F., The characterstc roots of certan real symmetrc matrces, Master s Thess, Unversty of Tennessee, 953. [] Gregory R.T., Karney D., A Collecton of Matrces for Testng Computatonal Algorthm, Wley- Interscence, 969. [] Akulenko L.D., esterov S.V., Hgh-Precson Methods n Egenvalue Problems and Ther Applcatons (seres: Dfferental and Integral Equatons and Ther Applcatons), Chapman and Hall/CRC, 004. [3] Ceselsk M., Blaszczyk T., umercal soluton of non-homogenous fractonal oscllator equaton n ntegral form, Journal of Theoretcal and Appled Mechancs 05, 53 (4), [4] Press W.H., Teukolsky S.A., Vetterlng W.T., Flannery B.P., umercal Recpes: The Art of Scentfc Computng (3rd ed.), Cambrdge Unversty Press, ew York 007.
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 informationNON-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 informationNUMERICAL 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 informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationConvexity 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 informationA 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 informationApplication 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 informationDETERMINATION 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 informationThe 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 informationProcedia Computer Science
Avalable onlne at www.scencedrect.com Proceda Proceda Computer Computer Scence Scence 1 (01) 00 (009) 589 597 000 000 Proceda Computer Scence www.elsever.com/locate/proceda Internatonal Conference on Computatonal
More informationThe Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites
7 Asa-Pacfc Engneerng Technology Conference (APETC 7) ISBN: 978--6595-443- The Two-scale Fnte Element Errors Analyss for One Class of Thermoelastc Problem n Perodc Compostes Xaoun Deng Mngxang Deng ABSTRACT
More informationLecture 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 informationNew 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 informationSTUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS
Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationGroup 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 informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More informationAn efficient algorithm for multivariate Maclaurin Newton transformation
Annales UMCS Informatca AI VIII, 2 2008) 5 14 DOI: 10.2478/v10065-008-0020-6 An effcent algorthm for multvarate Maclaurn Newton transformaton Joanna Kapusta Insttute of Mathematcs and Computer Scence,
More informationTHE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD
Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationMaejo 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 informationAppendix 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(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 informationAsymptotics 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 informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More informationNumerical 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 informationIrregular 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 informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More informationn α 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 informationON 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 informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationLecture 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 informationPART 8. Partial Differential Equations PDEs
he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef
More informationLecture 2: Numerical Methods for Differentiations and Integrations
Numercal Smulaton of Space Plasmas (I [AP-4036] Lecture 2 by Lng-Hsao Lyu March, 2018 Lecture 2: Numercal Methods for Dfferentatons and Integratons As we have dscussed n Lecture 1 that numercal smulaton
More informationResearch Article Cubic B-Spline Collocation Method for One-Dimensional Heat and Advection-Diffusion Equations
Appled Mathematcs Volume 22, Artcle ID 4587, 8 pages do:.55/22/4587 Research Artcle Cubc B-Splne Collocaton Method for One-Dmensonal Heat and Advecton-Dffuson Equatons Joan Goh, Ahmad Abd. Majd, and Ahmad
More informationFixed 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 informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationSolution 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 informationCHAPTER 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 informationHaar 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 informationSolving 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 informationNumerical Simulation of One-Dimensional Wave Equation by Non-Polynomial Quintic Spline
IOSR Journal of Matematcs (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 14, Issue 6 Ver. I (Nov - Dec 018), PP 6-30 www.osrournals.org Numercal Smulaton of One-Dmensonal Wave Equaton by Non-Polynomal
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
More informationNumerical 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 informationAn identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites
IOP Conference Seres: Materals Scence and Engneerng PAPER OPE ACCESS An dentfcaton algorthm of model knetc parameters of the nterfacal layer growth n fber compostes o cte ths artcle: V Zubov et al 216
More informationChapter 4 The Wave Equation
Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Finite Element Method - Jacques-Hervé SAIAC
COMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Fnte Element Method - Jacques-Hervé SAIAC FINITE ELEMENT METHOD Jacques-Hervé SAIAC Départment de Mathématques, Conservatore Natonal des Arts et Méters, Pars,
More informationALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION
ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LIX, 013, f.1 DOI: 10.478/v10157-01-00-y ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION BY ION
More informationCopyright 2014 Tech Science Press CMC, vol.43, no.2, pp.87-95, 2014
Copyrght 2014 Tech Scence Press CMC, vol.43, no.2, pp.87-95, 2014 Analytcal Treatment of the Isotropc and Tetragonal Lattce Green Functons for the Face-centered Cubc, Body-centered Cubc and Smple Cubc
More informationInternational 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 informationA 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 informationGlobal 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 informationLectures - 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 informationThe Finite Element Method: A Short Introduction
Te Fnte Element Metod: A Sort ntroducton Wat s FEM? Te Fnte Element Metod (FEM) ntroduced by engneers n late 50 s and 60 s s a numercal tecnque for solvng problems wc are descrbed by Ordnary Dfferental
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased
More informationFormal 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 informationDifference 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 informationA 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 informationModule 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 informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationRandić 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 informationHandout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system. and write EOM (1) as two first-order Eqs.
Handout # 6 (MEEN 67) Numercal Integraton to Fnd Tme Response of SDOF mechancal system State Space Method The EOM for a lnear system s M X + DX + K X = F() t () t = X = X X = X = V wth ntal condtons, at
More informationJournal of Universal Computer Science, vol. 1, no. 7 (1995), submitted: 15/12/94, accepted: 26/6/95, appeared: 28/7/95 Springer Pub. Co.
Journal of Unversal Computer Scence, vol. 1, no. 7 (1995), 469-483 submtted: 15/12/94, accepted: 26/6/95, appeared: 28/7/95 Sprnger Pub. Co. Round-o error propagaton n the soluton of the heat equaton by
More informationTransfer 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 information2 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 informationRELIABILITY ASSESSMENT
CHAPTER Rsk Analyss n Engneerng and Economcs RELIABILITY ASSESSMENT A. J. Clark School of Engneerng Department of Cvl and Envronmental Engneerng 4a CHAPMAN HALL/CRC Rsk Analyss for Engneerng Department
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12
REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationThe 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 informationYong 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 informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationPerron Vectors of an Irreducible Nonnegative Interval Matrix
Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of
More information1 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 informationA Robust Method for Calculating the Correlation Coefficient
A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal
More informationA Cartesian-grid integrated-rbf method for viscoelastic flows
Home Search Collectons Journals About Contact us My IOPscence A Cartesan-grd ntegrated-rbf method for vscoelastc flows Ths artcle has been downloaded from IOPscence. Please scroll down to see the full
More informationChapter 4: Root Finding
Chapter 4: Root Fndng Startng values Closed nterval methods (roots are search wthn an nterval o Bsecton Open methods (no nterval o Fxed Pont o Newton-Raphson o Secant Method Repeated roots Zeros of Hgher-Dmensonal
More informationTensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q
For orthogonal curvlnear coordnates, eˆ grad a a= ( aˆ ˆ e). h q (98) Expandng the dervatve, we have, eˆ aˆ ˆ e a= ˆ ˆ a h e + q q 1 aˆ ˆ ˆ a e = ee ˆˆ ˆ + e. h q h q Now expandng eˆ / q (some of the detals
More informationComputers 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 informationUsing 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 informationThe Discretization Process
FMIA F Moukalled L Mangan M Darwsh An Advanced Introducton wth OpenFOAM and Matlab Ths textbook explores both the theoretcal foundaton of the Fnte Volume Method (FVM) and ts applcatons n Computatonal Flud
More informationInexact 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 informationWorkshop: Approximating energies and wave functions Quantum aspects of physical chemistry
Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department
More informationThe internal structure of natural numbers and one method for the definition of large prime numbers
The nternal structure of natural numbers and one method for the defnton of large prme numbers Emmanul Manousos APM Insttute for the Advancement of Physcs and Mathematcs 3 Poulou str. 53 Athens Greece Abstract
More informationHongyi 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 informationA 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 informationSolving Singularly Perturbed Differential Difference Equations via Fitted Method
Avalable at ttp://pvamu.edu/aam Appl. Appl. Mat. ISSN: 193-9466 Vol. 8, Issue 1 (June 013), pp. 318-33 Applcatons and Appled Matematcs: An Internatonal Journal (AAM) Solvng Sngularly Perturbed Dfferental
More informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More informationNumerical Solution of two dimensional coupled viscous Burgers Equation using the Modified Cubic B-Spline Differential Quadrature Method
umercal Soluton of two dmensonal coupled vscous Burgers Equaton usng the odfed Cubc B-Splne Dfferental Quadrature ethod H. S. Shukla 1, ohammad Tamsr 1*, Vneet K. Srvastava, Ja Kumar 3 1 Department of
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationOn an Extension of Stochastic Approximation EM Algorithm for Incomplete Data Problems. Vahid Tadayon 1
On an Extenson of Stochastc Approxmaton EM Algorthm for Incomplete Data Problems Vahd Tadayon Abstract: The Stochastc Approxmaton EM (SAEM algorthm, a varant stochastc approxmaton of EM, s a versatle tool
More informationLOW 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 informationSupplementary Notes for Chapter 9 Mixture Thermodynamics
Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects
More informationThe Second Anti-Mathima on Game Theory
The Second Ant-Mathma on Game Theory Ath. Kehagas December 1 2006 1 Introducton In ths note we wll examne the noton of game equlbrum for three types of games 1. 2-player 2-acton zero-sum games 2. 2-player
More informationA 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 informationOn the Multicriteria Integer Network Flow Problem
BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 5, No 2 Sofa 2005 On the Multcrtera Integer Network Flow Problem Vassl Vasslev, Marana Nkolova, Maryana Vassleva Insttute of
More information