Jacobi Similarity Transformation for SVD and Tikhonov Regularization for Least Squares Problem: The Theoretical foundation

Similar documents
Parameter Estimation and Hypothesis Testing of Two Negative Binomial Distribution Population with Missing Data

Generalisation on the Zeros of a Family of Complex Polynomials

Lecture 3 summary. C4 Lecture 3 - Jim Libby 1

Maximum likelihood estimate of phylogeny. BIOL 495S/ CS 490B/ MATH 490B/ STAT 490B Introduction to Bioinformatics April 24, 2002

Multiquadrics method for Couette flow of a yield-stress fluid under imposed torques

The Stability of High Order Max-Type Difference Equation

Technical Appendix for Inventory Management for an Assembly System with Product or Component Returns, DeCroix and Zipkin, Management Science 2005.

Chapter Linear Regression

African Journal of Science and Technology (AJST) Science and Engineering Series Vol. 4, No. 2, pp GENERALISED DELETION DESIGNS

FRACTIONAL MELLIN INTEGRAL TRANSFORM IN (0, 1/a)

TWO INTERFACIAL COLLINEAR GRIFFITH CRACKS IN THERMO- ELASTIC COMPOSITE MEDIA

EGN 3321 Final Exam Review Spring 2017

Suppose we have observed values t 1, t 2, t n of a random variable T.

- 1 - Processing An Opinion Poll Using Fuzzy Techniques

Nonlocal Boundary Value Problem for Nonlinear Impulsive q k Symmetric Integrodifference Equation

ANSWERS TO ODD NUMBERED EXERCISES IN CHAPTER 2

International Mathematical Forum, Vol. 9, 2014, no. 13, HIKARI Ltd,

-HYBRID LAPLACE TRANSFORM AND APPLICATIONS TO MULTIDIMENSIONAL HYBRID SYSTEMS. PART II: DETERMINING THE ORIGINAL

Convergence tests for the cluster DFT calculations

Exact Moments of Record Values from Burr Distribution. with Applications

Applying p-balanced Energy Technique to Solve Liouville-Type Problems in Calculus

Calculus 241, section 12.2 Limits/Continuity & 12.3 Derivatives/Integrals notes by Tim Pilachowski r r r =, with a domain of real ( )

ME 501A Seminar in Engineering Analysis Page 1

BINOMIAL THEOREM OBJECTIVE PROBLEMS in the expansion of ( 3 +kx ) are equal. Then k =

Integral Equations and their Relationship to Differential Equations with Initial Conditions

LIPSCHITZ ESTIMATES FOR MULTILINEAR COMMUTATOR OF MARCINKIEWICZ OPERATOR

4. Runge-Kutta Formula For Differential Equations. A. Euler Formula B. Runge-Kutta Formula C. An Example for Fourth-Order Runge-Kutta Formula

_ J.. C C A 551NED. - n R ' ' t i :. t ; . b c c : : I I .., I AS IEC. r '2 5? 9

One of the common descriptions of curvilinear motion uses path variables, which are measurements made along the tangent t and normal n to the path of

On the hydrogen wave function in Momentum-space, Clifford algebra and the Generating function of Gegenbauer polynomial

On the Quasi-Hyperbolic Kac-Moody Algebra QHA7 (2)

Modeling and Predicting Sequences: HMM and (may be) CRF. Amr Ahmed Feb 25

The sphere of radius a has the geographical form. r (,)=(acoscos,acossin,asin) T =(p(u)cos v, p(u)sin v,q(u) ) T.

4. Runge-Kutta Formula For Differential Equations

Exponential Synchronization of the Hopfield Neural Networks with New Chaotic Strange Attractor

DUALITY IN MULTIPLE CRITERIA AND MULTIPLE CONSTRAINT LEVELS LINEAR PROGRAMMING WITH FUZZY PARAMETERS

Studying the Problems of Multiple Integrals with Maple Chii-Huei Yu

( ) ( ) Weibull Distribution: k ti. u u. Suppose t 1, t 2, t n are times to failure of a group of n mechanisms. The likelihood function is

Universitatea Transilvania din Braşov HABILITATION THESIS

1. Consider an economy of identical individuals with preferences given by the utility function

Conquering kings their titles take ANTHEM FOR CONGREGATION AND CHOIR

Boyce/DiPrima 9 th ed, Ch 7.6: Complex Eigenvalues

() t. () t r () t or v. ( t) () () ( ) = ( ) or ( ) () () () t or dv () () Section 10.4 Motion in Space: Velocity and Acceleration

APPLICATION OF A Z-TRANSFORMS METHOD FOR INVESTIGATION OF MARKOV G-NETWORKS

Integral Solutions of Non-Homogeneous Biquadratic Equation With Four Unknowns

ON THE EXTENSION OF WEAK ARMENDARIZ RINGS RELATIVE TO A MONOID

Fakultas Matematika dan Ilmu Pengetahuan Alam, Institut Teknologi Bandung, Bandung, 40132, Indonesia *

SOME REMARKS ON HORIZONTAL, SLANT, PARABOLIC AND POLYNOMIAL ASYMPTOTE

Difference Sets of Null Density Subsets of

T h e C S E T I P r o j e c t

X-Ray Notes, Part III

ON TOTAL TIME ON TEST TRANSFORM ORDER ABSTRACT

The Solution of Heat Conduction Equation with Mixed Boundary Conditions

The Eigenvalue Problem of the Symmetric Toeplitz Matrix

Giant density fluctuations and other out-of-equilibrium properties of active matter. Francesco Ginelli

STOCHASTIC CALCULUS I STOCHASTIC DIFFERENTIAL EQUATION

Asymptotic Behavior of Solutions of Nonlinear Delay Differential Equations With Impulse

Interval Estimation. Consider a random variable X with a mean of X. Let X be distributed as X X

INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY

A TREATY OF SYMMETRIC FUNCTION PART 1

Reza Poultangari a a, b, * Mansour Nikkhah-Bahrami

On Certain Classes of Analytic and Univalent Functions Based on Al-Oboudi Operator

Duration Notes 1. To motivate this measure, observe that the duration may also be expressed as. a a T a

Brownian Motion and Stochastic Calculus. Brownian Motion and Stochastic Calculus

Physics 15 Second Hour Exam

Handout on. Crystal Symmetries and Energy Bands

Calculation of Effective Resonance Integrals

P-Convexity Property in Musielak-Orlicz Function Space of Bohner Type

Rotations.

Chapter 17. Least Square Regression

Decompression diagram sampler_src (source files and makefiles) bin (binary files) --- sh (sample shells) --- input (sample input files)

BEST PATTERN OF MULTIPLE LINEAR REGRESSION

Forward Kinematics of a Stewart Platform Mechanism

8. Queueing systems lect08.ppt S Introduction to Teletraffic Theory - Fall

Chapter Simpson s 1/3 Rule of Integration. ( x)

Patterns of Continued Fractions with a Positive Integer as a Gap

4.1 Schrödinger Equation in Spherical Coordinates

Laplace Transform. Definition of Laplace Transform: f(t) that satisfies The Laplace transform of f(t) is defined as.

On Fractional Operational Calculus pertaining to the product of H- functions

The z-transform. LTI System description. Prof. Siripong Potisuk

5-1. We apply Newton s second law (specifically, Eq. 5-2). F = ma = ma sin 20.0 = 1.0 kg 2.00 m/s sin 20.0 = 0.684N. ( ) ( )

Generalized Fibonacci-Type Sequence and its Properties

Eurasian International Center of Theoretical Physics, Eurasian National University, Astana , Kazakhstan

= y and Normed Linear Spaces

3. REVIEW OF PROPERTIES OF EIGENVALUES AND EIGENVECTORS

Reflection from a surface depends on the quality of the surface and how much light is absorbed during the process. Rays

). So the estimators mainly considered here are linear

Bayesian Estimation of the parameters of the Weibull-Weibull Length-Biased mixture distributions using time censored data

Copyright Birkin Cars (Pty) Ltd

P a g e 5 1 of R e p o r t P B 4 / 0 9

Certain Expansion Formulae Involving a Basic Analogue of Fox s H-Function

Observations on the transcendental Equation

The Complete Graph: Eigenvalues, Trigonometrical Unit-Equations with associated t-complete-eigen Sequences, Ratios, Sums and Diagrams

Key words: Fractional difference equation, oscillatory solutions,

On Some Integral Inequalities of Hardy-Type Operators

On Probability Density Function of the Quotient of Generalized Order Statistics from the Weibull Distribution

Available online through

Upon Parametric Sensitivity Used in Damping Active Control for Human Body Protection Against Vibrations

Solution of Impulsive Differential Equations with Boundary Conditions in Terms of Integral Equations

Addition & Subtraction of Polynomials

Transcription:

IOSR Joul of Mhemcs IOSR-JM e-issn: 78-578 -ISSN: 39-765X Volume Issue 5 Ve VII Se - Oc6 PP 76-85 wwwosoulsog Jcob Smly sfomo fo SVD d hoov Regulzo fo Les Sues Poblem: he heoecl foudo Sehe Ehdmhe Uwmus Deme of Mhemcl Sceces Kog Se Uvesy Aygb Kog Se Absc: he e eses soluo o Les sues euo s dve by Sgul Vlue Decomoso SVD obed fom usg Jcob Smly sfomo whou ecouse o hoov egulzo mee We use he LU decomoso d QR lgohm s bses of comso fo he obed esuls A olyoml f of ode fve fo he Les sues oblem ws doed Smle umecl d oblem ws obed s my souce od Kog Se Uvesy Aygb showg he fucol elosh bewee emeue d Relve Humdy I s suggesed h he SVD s umeclly bcwd sble s esed o by he LU decomoso mehod elve o sml esuls comued by QR d ohe ow mehods Keywods: les sues euo sgul vlue decomoso Jcob smly sfomo lgohm vese oblems AMS subec Clssfco 65G4 65F5 I Ioduco I he e we cosde les sues euo s lcble o modellg my Scefc oblems usg SVD obed fom Jcob Smly sfomo []Les sues euo belogs o he clss of mos oweful vese oblems whch dels wh he ocess of clculg fom gve se of obsevos he cusl fcos whch oduce hem Such suos wee ledy ofe ecoueed ocs d couscs commuco heoy d lguge ocessg clculo of he desy of he eh mesuemes of s gvy feld medcl mgg comue vso oceoghy soomy emoe sesg mche leg d o-desucve esg [] Bsc ools fo dscussos e vese fuco heoems fo se vlued ms eguly os fo Lschz fuco chcezo of gel d mec eguly d Lschz behvou of mulfuco coeco wh closed cove gh heoem[3 4] he dsce bewee eguly d sog eguly of se vlued m F : X Y whch he Jcob m my be vesged fo oeess oologcl sce foms he bss of ou dscusso hs c be vefed usg heoem due o [5] whch ses h evl m [A] s sogly egul f md A] ] d [ ] [ A whch my be we u o fco whee s he sgul vlue of he m A decesg ode he elve dsce fo whch m A s sgul o he se of sgul A s s gul A mces some om s defed s ds A m A A he ecocl of codo umbe c be eeed s mesue of eess o sguly of he m A d heefoe ds A K A As esul [6] oved h fo such sgul m A A hee holds he esme: A h s fo b A A A K A b A A b A K A we hve Ioduced hus o he dscusso s he well ow Holde s om defed he fom DOI: 979578-577685 wwwosoulsog 76 Pge

Jcob Smly sfomo fo SVD d hoov Regulzo fo Les Sues Poblem: Whee s udesood [7] h C 3 4 C m e d e e e e s he cocl bss By fuhe eoso o Holde euly would yeld h: e R hs mles lm ; fo e hs bgs us o he coce of Lschz couy fo f me logous o Rem egble fuco [8] whch eles h fo evey om R hee s b g d g d b such h : f [ ] f d f d Md M I follows h fo cove fucol defed o D R we would eec h g f s couously dffeeble o ] [ fo whch f f g g g d f f [ ][ ] Wheefom esos due o [9] mlced h f f f [ ] M em vld Moe moly whe he m A s sue m he loghm of A s defed by he euo log A A I A I I d d s codo umbe K A K A log A log he emg he e s cegozed s follows: I seco he secl decomoso v Jcob smly sfomo s gve heoecl foudo he Gves ohogol m le oo fo he cosuco of oul QR decomoso s useful ool fo umecl comuo s hghlghed Seco 3 he e dscusses he sgul vlue decomoso SVD fo he les sues oblem s well s hoov egulzo mees he cse of sguly I seco4 we demose wh umecl emle wh descbed mehods hee mgh be eed o comue zeos of olyomls; heefoe seco 5 he e eses evl Newo-Le mehods fo smuleously efg ll zeos of olyoml euo wh cul efeece o he wo of [] Numecl lluso s demosed wh he mehods seco 6 Flly cocluso s dw bsed o he segh of ou fdgs II he Secl Decomoso v Jcob Smly sfomo he moce of comug egevlues d egevecos of el o comle m co be oveemhszed ou esech whch s ofe me hyscs d egeeg cces good emle s vbo oblems d lyss of vces sscs he Jcob smly sfomo ws deved by Jcob 946 d s lcble o el symmec m A [ ] I he cse he m s o symmec oe mus fs sfom he m A o symmec cse he Jcob smly sfomo m s bsed d DOI: 979578-577685 wwwosoulsog 77 Pge

Jcob Smly sfomo fo SVD d hoov Regulzo fo Les Sues Poblem: o he heoem of cl es whch eles h ech symmec m A c be sfomed by mes of ohogol m Q o sml dgol m D such h Q A Q dg he m Q s ohogol d s defed he fom: he lgohm fo Jcob smly sfomo s evously used [] s fuhe descbed s follows u o he dscusso: Fo m of ode : Defe ode of ccucy If he we egevlues e DOI: 979578-577685 wwwosoulsog 78 Pge so Else efom he followg oeos begg wh se Choose v If Se m wh < 4 Else ob he vlue fo by seg Whee d sg of sg of 4 v Ob he ew m A fom he Old m A by usg he followg udg fomul: Fo ech whee d Se cos s By symmey cos s By symmey cos s cos s By symmey

Jcob Smly sfomo fo SVD d hoov Regulzo fo Les Sues Poblem: s s cos cos By symmey e cos s cos s by v d symmey Elsewhee v Ree oeos sg fom se Fsh Fsh Ed I mus be oed ou h he sue oo of dgol m mus be flly comued ssumg o symmec m ws oglly sfomed o symmec m s he lgohm bove I well develoed sese he Gves Roo Ohogol m [] whch s efed Jcob Smly sfomo obed by Gves 954 fo develog he QR decomoso of symmec m whee cos d s esecvely s comued he fom: cos s Moe defve fuhe use of Gves oo s educo of symmec m A o -dgol fom whee Sum seuece becomes hd ool fo he culzo of olyoml euo o whch y of Newo s mehods s lcble he gves oo s fe bu Jcob s s fe fo ehusve dscusso o hs ese edes e efeed o []We lso oed ou h Householde Refleco c be used s leve ug comuol cos o cosdeo All he sme Gves le oo s wohwhle he sudes of oboc ms mulos III he Les Sues Poblem he dscezo of ll-osed oblem ofe leds o le les sues oblem he fom Fd m m A b such h m A b A R s he les sues soluo of he oblem I mos eemel wo whch s ofe he cse cce he m A my hve lge umbe of sgul vlues vey close o zeo hs ushes he ose whch s ese b he gh hd sde o be mlfed he seudo-vese soluo A b 3 wh esul effec of huge codo umbe d heeby edeg he ome soluo useless [3]We hus oduce he Sgul vlue decomoso SVD o he comuo s dydc decomoso of A : A U V dg whee U U U ] U s m [ m m wh ohogol colum vecos whee V s uy m As esul he soluo ocess by SVD s he fom: c v whee c u b 3 I he bsece of e sguly o sguly SVD s umeclly bcwd sble O he ohe hd f hee s sguly he soluo sysem hee s hus eed fo develome of good umbe of egulzo echues [4] fo whch comes o focus he ychoov egulzo h elces soluo of ogl sysem m A m b L A R LR by he eby sysem A A I A 34 b 33 DOI: 979578-577685 wwwosoulsog 79 Pge

Jcob Smly sfomo fo SVD d hoov Regulzo fo Les Sues Poblem: he mee s egulo whch cools he wegh gve o he mmzo of he mmzo of he esdul A L L elve o d b A I s suosed h L m he m L s ce m whch my be e s dey m o dscee omo o some devve oeo [6]Noe h he mmzo of oblem 33 he ssumo Null A Null L O s mo fo useful uose beg ohogol d R ue dgol he ue gul mes fo hee dffee cses e ouled below: I he me beg we defe he QR decomoso wh m Q m R m< R mm m 3 m> R o O m I oe ssg he eed hoov mehod wh esdul [6d 7] s he euo A A I A b A = 35 I he mlemeo of mehod 35 s dvsble s sg o by g A A I A b d Ideed mhemclly euo 35 hs fed cocve guge [5] ecessy codo fo Lschz couy fo whch he Bch Fed o heoem s deely mlced We shll o delve o dels hee hs dscusso DOI: 979578-577685 wwwosoulsog 8 Pge

Jcob Smly sfomo fo SVD d hoov Regulzo fo Les Sues Poblem: IV Numecl Eemes Emle he eemel d wee g s my souce fom RODAN Aygb Kog Se Uvesy s showed ble SN DAEIME -4-4 3-4 -8-8 -8 - - - - - 3- - - 3- -7-7 -7-6 -6-3 -3 3-3 -5-5 - - - - - - -9-9 -9 ble EMPERAURE X 396 333634 9887573 76463 643887 6986 33966 8355457 863878 367793 3435 363575 96744 886789 9588 648773 7984 69459 8987 787376 3395 334677 333333 363983 936996 947345 8766 889984 8458 7455 74785 75395 6664 659674 O RELAIVE HUMIDIY % Y 6796866 684359 749686 87566 83368 8345398 475699 373 59693 666683 6444 55738473 39557 4495453 4797735 835739 7994463 844995 789 7953775 636848 5478859 653497 73495656 753975 7653 66996 7698347 87376 86998 894369 83485 899554 835695 I s suggesed h emeue s deede vble X d h Relve Humdy s he deede vble Y We use wdows 7 veso of MALAB fo he clculos d se of esuls flog o hmec s dslyed ble ble showg Resuls comued fo he cosdeed mehods decomoso QR decomoso Sgul vlue hoov decomoso mehod LU mehod fo 3 ^ e+6* 467 449 64 5 mehod fo 3 ^ e+9* 34 6977 9778 3 ^ e+6* 49 38 59 5 Regulzo mehod 34 wh [] ^ 3 6 3 6 e * DOI: 979578-577685 wwwosoulsog 8 Pge

Jcob Smly sfomo fo SVD d hoov Regulzo fo Les Sues Poblem: I ble comued esuls fo LU decomoso d Sgul Vlue decomoso SVD e umeclly bcwd sble wh hgh degee of closeess whe comed wh esuls fom QR Fcozo O he ohe hd he hoov egulzo mehod seemed o hve efomed vey ooly hs cse whch ws ue ueeced We oo he mee he evl of [ ] d foud ou h hee ws o sgfc dffeece he esuls s ll comued esuls by hoov eguled mee mehod he evl sed wee he sme he gh llusg he bove esuls c be demosed usg MALAB lo fuco s dslyed Fgues -3 Ech of he mehods he fgues ws efeed o Sees 3 d 4 esecvely h s o sy he LU mehod Sees QR mehod sees SVD mehod sees 3 d hoov egulzo mehod sees 4Comuolly he ll-osed oblem mlced h A A he secum of A A h fo some hee s m hece he m A A cl vecos deve he Jod fom of dco h egevlues of A A [] A A e ubouded hs mes s degeeed hus he egevecos d Fgue Fgue DOI: 979578-577685 wwwosoulsog 8 Pge

Jcob Smly sfomo fo SVD d hoov Regulzo fo Les Sues Poblem: Fgue 3 V Comuo of olyoml zeos by he Ievl Newo Oeo Suosg he coeffces eg he olyoml f fo he les sues oblems e eessed he fom of ucey such h: P whee e el o comle 5 If we omlze he bove olyoml euo moc olyoml wh ledg coeffce uy s obed he Newo oeo fo sgle vble euo s lcble o he wo s gve he fom 5 Le be o emy cove subse of Bch sce d le be coco of o self If hee s N such h N : N comc se emy closed cove se K such h: K K d N d N hs les wo os he hee ess o K K s he comleme of I s suosed h of os s seed h s hey e K sce d Ldelof codo fo s : R eeds o couous fuco o If we ele closely comubly om holds By eze eeso heoem would follow h s oml d closed hus evey couous fuco wh Hh-Bch eeso heoem he we would l eguly sces wh omly v Uysoh s mezo heoem h s ffoded by ychooff so h he closed gh heoem d ufom boudedess cle he lls of Be s cegoy e foce fo whch he bsc cle fo fdg zeos of olyoml holds I s eeced he u bll B * D s lso comc he we s oology he heoecl foudo s ele meoed he begg dedced o he wo [] s s follows Fo moc olyoml of degee whose zeos of esecve mullces e ow we defe h P 5 Ioduce he followg oos he sese of [67] h: v v v ; DOI: 979578-577685 wwwosoulsog 83 Pge

Jcob Smly sfomo fo SVD d hoov Regulzo fo Les Sues Poblem: DOI: 979578-577685 wwwosoulsog 84 Pge ; ; By fuhe usg des due o [96d 7] would be obed h 53 he fed o elo of Lgueel s ye mehod c be develoed [6] wheefom: 54 I he cse of smle zeos fomul eessed euo 54 s he fml Lgueel s mehod of hd ode of covegece he choce of he wo oos o be e fom he wo dss euo 54 s guded by he dscusso deled [6] Imlemeo of he bove Ccul evl hmec by smuleously efg he zeos of he olyoml euo P whee he mdo d dus of he efg ds should esec d g o cosdeo h ^ s ew omo o he zeos of P s we he fom: f 55 md I I Euo 54 S S S S f We oed h veso of he eve mehod of Hlley ye he cse of smle zeos euo 5 s gve by he euo: = 56 he umbe of mullces of he oos gve ds R ID he sese of [7] s he Lgouell lmg fomul s deved [8] he fom: lm lm lm d d D d d 56 Whee deoes he ege closes o he el umbe

Jcob Smly sfomo fo SVD d hoov Regulzo fo Les Sues Poblem: VI Cocluso he e suded use of sgul vlue decomoso obed fom Jcob smly sfomo fo he soluo of les sues oblem We comed oes wh sml esuls comued by he LU Fcozo d QR lgohm Smle umecl emle ws e s my souce d fom RODAN Cee fo Lowe Amoshec Sudes Kog Se Uvesy Aygb Kog Se Nge fo he emeue d Relve Humdy fo eod of hy fou mohs whch ws he me ge he Cee fo Lowe moshec Sudes ws esblshed A olyoml f of ode fve ws used Numecl esuls comued by he bove med mehods e gve ble We lso llused he esuls comued fo he uose of foecsg s show Fgues -3 I s suggesed h hee s sog coelo he d s see fom Fgue 3 s esul of coley bewee LU decomoso d SVD A effce evl bsed mehod fo smuleously efg zeos of olyoml whe coeffces e eessed he fom of ucey ws dscussed he sese Peovc [ 6] fo he uose of esy ccessbly o he edes I comug wh Jcob mehod oe fs sfoms o symmec m o symmec m I s ecommeded h whe m A s ely dgol Jcob smly sfomo mehod wll lwys hve ue hd ove QR mehod ovded h oe sog ceo s secfed fo Jcob s mehod [6] Refeeces [] Wlso JH 965: he Algebc Egevlue Poblem Cldo Pess Ofod [] Wled : Ivese Poblem he fee ecycloed ; Auhoy Cool:NDL:57748 Jue 5 [3] Uwmus SE 5 Gh Comleo cluso sooe fo Ievl Les Sues Euo Amec Joul of Mhemcs d Sscs 5 4-3 [4] Neume A 99: Ievl mehods fo sysems of euos Cmbdge Uvesy Pess Cmbdge [5] Rum SM 996: he dsce bewee eguly d sog eguly I G Alefeld d B Lg edos Scefc Comug d Vlded Numecs Vol 9 of Mhemcl Resech P 5-7 [6] Boc A 9: Numecl Mehods Scefc Comug Vol SIAM Phldelh [7] L W 8: Lecue Noes of M Comuos Deme of Mhemcs Nol sg Hu Uvesy Hschu w 343ROC [8] Oeg JM Rhebold WC : Ieve Soluo of Nole euos sevel vbles Clsscs Aled Mhemcs SIAM Phldelh USA [9] Neume A: Ioduco o Numecl lyss Cmbdge Uvesy PessCmbdge [] Roh J 5: A Hd boo of esuls o evl le oblems Czech Acdemy of Sceces Pgue Czech Reublc Euoe Uo;5wwwcscsczoh [] Peovc MRcc L d Mlosevc D 6:Lgueel-Le mehods fo he smuleous omo of olyoml mulle zeos; Yugoslv Joul of Oeos Resech 6 3-44 [] Soe J d Bulsch R 98: Ioduco o Numecl Alyss Sge Velg New Yo [3] Uwmus SE d Ou FO : Comuo of Egevlues of Hem M v Gves Ple Roos Nge Joul of Aled Scece Vol P -6 [4] Golub G d V-Lo F983:M Comuoshe Johs Hos Uvesy PessBlmoe Myld [5] Sef W 8: ol vo egulzo fo le ll-osed vese oblems: Eesos d lco PhD hess Gdue College Azo Se Uvesy USA [6] Uwmus SE4:Regulzo of No-smooh fuco whose Jcob s ely sgul o sgul Amec Ieol Joul of Coemoy Scefc Resech-8 ISSN 349-445 [7] Mye G 998: Eslo-flo wh cocve evl fucos Alcos of Mhemcs 434 4-54 [8] Peovc MS 989: Ieve mehod fo smuleous cluso of olyoml zeos 387 Sge Velg [9] Peovc MS d covc S 995: Zeo fdg mehods of Fouh Ode Joul of Comuol d Aled Mhemcs No 64 P 9-94 [] Fme MR d Lozu G 977: A Algohm fo he ol o Pl Fcozo of olyoml Mhs Poc Cmb PhlSoc No 8 P 47-437 DOI: 979578-577685 wwwosoulsog 85 Pge