Estimation Approach to Ratio of Two Inventory Population Means in Stratified Random Sampling

Size: px
Start display at page:

Download "Estimation Approach to Ratio of Two Inventory Population Means in Stratified Random Sampling"

Transcription

1 American Journal of Operational Researc 05, 5(4): 96-0 DOI: 0.593/j.ajor Estimation Approac to Ratio of Two Inventor Population Means in Stratified Random Sampling Subas Kumar Yadav, S. S. Misra,*, Alok Kumar Sukla Department of Matematics and Statistics (A Centre of Excellence), Dr. RM Avad Universit, Faizabad, U.P., India Department of Statistics, D.A-V, College, Kanpur, U.P., India Abstract Te present paper proposes a new approac of estimation to te inventor sstem wic undergoes te process of stratification b using te tecnique of stratified random sampling. Tis seeks to develop te estimation of ratio of two population means using auxiliar variable under stratified random sampling troug generalized ratio tpe estimator. Te statistical analses of existing estimators and proposed estimator ave been discussed and condition of efficient estimator as been attained. Wit te elp of computing algoritm, a numerical illustration of te problem as been also presented to meet te scientific footings of teor of estimation. Kewords Stratified Inventor, Ratio Estimator, Auxiliar Variable and Computing Algoritm. Introduction Usable and idle resources referring to man, material, macine and mone are called inventor items; vide Ackoff and Sasieni (993). Heterogeneous inventor items are stratified for optimal control of inventor sstem. Stratification of inventor sstem attempts to provide te mecanism to optimize te inventor control for eac stratum wic is believed to be ver efficient and easil minimizes te costs. Tis also makes andling of inventor eas for eac stratum and accordingl te stratum domain experts are supposed to supervise te field works of te inventor andling and management, vide for example Alex and Benn (009). In tis situation, stratified random sampling is used to develop te efficient estimator to estimate te population parameters of stratified inventor items for caracteristics under stud is not omogeneous, vide for example Cocran (940, 977). A few researces are available in tis field and terefore demands deep exploration in tis area of interdisciplinar researc. Under tis sampling sceme, te wole population of inventor items is divided into relativel omogeneous groups. Ten te required sample is drawn b taking appropriate subsamples from tese stratums using simple random sampling tecnique. Te fres objective of tis paper is to develop estimation approac to te ratio of two population means in stratified inventor sstem b evolving te efficient estimator using * Corresponding autor: sant_x003@aoo.co.in (S. S. Misra) Publised online at ttp://journal.sapub.org/ajor Coprigt 05 Scientific & Academic Publising. All Rigts Reserved stratified random sampling tecnique for inventor populations under consideration. Tis proposes estimation of ratio of two population means using auxiliar variable under stratified random sampling troug generalized ratio tpe estimator. Te expressions for te bias and mean square error (MSE) ave been obtained up to te first order of approximation. Te minimum MSE of proposed estimator is also obtained. An empirical stud is also carried out to meet te teoretical findings.. Statistical Analsis Te estimation of te ratio of two population means as been studied b man autors in te literature including Murt and Sing (965), Rao and Pareira (968), Sa and Sa (978), Ra and Sing (985), Upadaa and Sing (985), Upadaa et.al (985), Sing and Rani (005, 006) and Sindu et. al (009) etc. et us consider te population of inventor items under stud consisting of units and tis population of size is divided into stratums eac of size (,,..., ) and te required sample of size n is drawn b taking n (,,..., ) units from corresponding strata. Tus = and n= n. et 0 and be te two main variables under stud and x be te auxiliar variable. 0i, i and x i (,,...,, i =,,..., ) are te observations on te i t unit of te t stratum for te variables 0, and x respectivel.

2 American Journal of Operational Researc 05, 5(4): ow we ave te following notations wic ave been used trougout tis paper. Y0 = 0i, te t strata mean for te main i = variable under stud 0. Y = i, te t strata mean for te main i = variable under stud. X = xi, te t strata mean for te auxiliar i = variable under stud x. Y0 = Y 0 = WY 0, te inventor = = population mean for te main variable under stud 0. Y = Y = WY, te inventor = = population mean for te main variable under stud. X= X = WX, te inventor = = population mean for te auxiliar variable under stud x. n 0 = 0i n, te sample mean of 0 for t i = strata. n = i n, te sample mean of for t strata. i = n x = xi n, te sample mean of x for t strata. i = W =, te weigt of t strata. Y0 R =, te ratio of two population means of stud Y variables. It is well known tat te appropriate estimators for te estimation of inventor population parameters are te corresponding statistics. Tus appropriate estimators for population means Y 0, Y and X are te usual unbiased estimators, te sample means in stratified random sampling = W, st = W and 0st 0 x = Wx respectivel. st Te conventional estimator for te ratio of two population means R in stratified random sampling is defined as, ˆ 0st Rst = (.) st As it is well known tat te use of auxiliar information supplied b te auxiliar variable enances te efficienc of te estimator of an parameter, so Sing (965) proposed te traditional ratio tpe estimator of ratio of two population means using auxiliar information in stratified random sampling as, 0st X G st = (.) st xst Te bias and mean square error of G st, up to te first order of approximations respectivel are, S Sx S0 + Y X YY 0 st = γ (.3) S0x S x BG ( ) R W + + YX 0 YX S0 S Sx + + Y0 Y X st = γ S0 S0x Sx MSE( G ) R W + YY Y X YX 0 0 (.4) Jan et.al (04) proposed te following modified estimator using median (M d ) of te auxiliar information as, G st X + M = 0st d st xst (.5) Te bias and mean square error of G st, up to te first order of approximation is, S S x + Y ( X + Md ) S0 S 0x BG ( st ) = R W γ + (.6) YY 0 Y0( X+ Md ) S x + Y ( X + Md )

3 98 Subas Kumar Yadav et al.: Estimation Approac to Ratio of Two Inventor Population Means in Stratified Random Sampling S0 S Sx + + Y0 Y θ st = γ S0 S0x Sx MSE( G ) R W were, θ = ( X + M d ), 0 = 0i 0 = i x = ( i ) ( ) S, ( ) S, S x x, + YY Yθ Yθ 0 0i 0 i 0 0 (.7) S = ( )( ), S = ( )( x x ), 0x 0i 0 i S = x x ( )( ) x i i 3. Statistical Analsis of Proposed Estimator Motivated b Jan et.al (04), we propose te following generalized estimator for te ratio of two inventor population means as, α 0st X + Md τ st = (3.) st xst + Md Were α is a suitable constant be determined suc tat te mean square error of τ st is minimum. In order to stud te large sample properties of te proposed estimator we ave assumed tat, 0 = Y0( + e0), = Y( + e) and x = X( + e) suc tat Ee ( 0) = Ee ( ) = Ee ( ) And S0 Ee ( 0) Wγ Y0 S Ee ( ) Wγ Y =, =, Sx Ee ( ) = Wγ, X S0 Ee ( 0e) = Wγ, YY 0 S0x Ee ( 0e) = Wγ, YX 0 Sx Ee ( e) = Wγ. YX Using above expressions, te bias and te mean square error of te proposed estimator, up to te first order of approximation respectivel are, S α( + αθ ) Sx + Y X st = γ S0 αθs0 x αθs x B( τ ) R W + YY 0 Y0X YX (3.) S0 S αθ S x + + Y0 Y X S ( ) 0 S MSE τ 0x st = R W γ + αθ YY 0 Y0X Sx αθ YX (3.3) X were, θ =. X + Md Tis mean square error is minimum for, S0x Sx Wγ YX 0 YX α = Sx θ Wγ X And te minimum mean square error is, S0 S S0 Wγ + Y0 Y YY 0 MSEmin ( τ st ) = R S0x Sx Wγ YX 0 YX Sx Wγ X (3.4) (3.5)

4 American Journal of Operational Researc 05, 5(4): Efficienc Comparison From (.4) and (3.5) we ave tat te proposed estimator τ st is better tan te estimator G st, if MSE( G st ) MSEmin ( τ st ) > 0, if S0x Sx W γ YX 0 YX Sx R Wγ 0 > (4.) X Sx S0x Sx + Wγ X YX 0 YX Under te above laid down condition, te proposed estimator is believed to perform more efficientl tan te estimator G st considered in (.) of ratio of two population means for it will ave lesser mean squared error as compared to G st. From (.7) and (3.5) we ave tat te proposed estimator τ st is better tan te estimator G st, if MSE( Gst ) MSEmin ( τ st ) > 0, if S0x Sx W γ YX 0 YX Sx R Wγ 0 > (4.) X Sx S0x Sx + Wγ θ Y0θ Yθ Te proposed estimator τ st given in (3.) is supposed to provide to lesser mean squared error as compared to G st under te laid down condition in (4.). In brief, te psical significance of expressions given in (4.) and (4.) attempts to seek exceedingl better estimator aving least mean square error as compared to previousl existing ones. 5. Computing Algoritm and umerical Illustration Te following algoritm as been developed to compute te estimator and its efficienc. i. Begin ii. Data input iii. Compute sample mean of first inventor population iv. Compute sample mean of second inventor population v. Compute sample and population means of auxiliar inventor population vi. Compute estimator for ratio of two inventor populations vii. Compute biases for all estimators viii. Compute MSE ix. Compute efficienc (Percentage Relative Efficienc-PRE) x. If PRE is greater tan previous ones xi. Find efficient estimator xii. Data output xiii. End o Start Input Data Define Means Input Sample Data Compute MSE Compute Diff Is diff greater tan zero Yes Condition Attained Find Estimator Input for PRE Compute PRE Output PRE End Figure 5.. Tabular form of te algoritm (Computing flow cart)

5 00 Subas Kumar Yadav et al.: Estimation Approac to Ratio of Two Inventor Population Means in Stratified Random Sampling Table 5.. Computed Data Statistics Y 0 = 9.40 Y 0 = 09.0 Y = Y = X = 3757 X = S 0 ( ) = S 0 ( ) = S ( ) = S ( ) = S x = S x = 6.79 S 0 = S 0 = 45.0 S 0x = S 0x = S x = S x = Md ( x ) = Table 5.. Computed Bias, MSE and PRE of different estimators w.r.t. G st Estimators BIAS M S E PRE G st.4e e G st 8.4E E τ st -9.94E E To meet out te teoretical findings, we ave considered te data in Murt (967) were te main variables 0 and are te number of workers and te fixed capital respectivel along wit te output as te auxiliar variable. Te size of te population is 0 and is divided into two stratums eac of size 5. Te sample size is 4 b taking from eac stratum. Following are te parameters computed. 6. Observations and Conclusions Wen inventor items are eterogeneous and uge, te estimation approac is onl panacea for estimating te caracteristics of inventor populations under consideration wic oterwise seems difficult to control and manage for an organization. So for te estimation of ratio of two inventor population means ave been discussed in paper, one is man-inventor population and anoter is mone-inventor population as defined b Ackoff and Sasieni (993). Under te section numerical illustration, two eterogeneous inventor populations are given wose estimators are given tereafter final estimator as ratio as been developed wose bias and mean square errors are lesser tan all previous estimators and percentage relative efficienc is given wic is iger tan previous ones, vide table 5., row 3. Tus proposed estimator is most efficient estimator of ratio of two inventor population means using auxiliar information among suc estimators as it as lesser mean square error. Tus te proposed estimator sould be preferred for te estimation of ratio of two inventor population means. REFERECES [] Ackoff R and Sasieni M W (993), Fundamentals of Operations Researc, Wile Eastern td. [] Cocran W G (940), te estimation of ields of cereal experiments b sampling for te ratio of grain to total produce, Jour. Agri. Sci., 59, 5-6. [3] Cocran W G (977), Sampling Tecniques, Wile Eastern td, tird edition. [4] Gerskov Alex and Moldovanu Benn (009), Dnamic Revenue Maximization wit Heterogeneous Objects: A Mecanism Design Approac, American Economic Journal: Microeconomics, Volume, umber, pp (3). [5] Jan, R Maqbool, S Amad, A and azir, A (04), Modified Ratio Tpe Estimator of Two Population Means in Stratified Sampling, Indian Streams Researc Journal, 4,5, -6. [6] Murti M and Sing M P (965), On te estimation of ratio and product of te population parameters, Sanka, B, 7, [7] Rao J..K. & Pereira.P. (968), On double ratio estimators, Sanka, A, 30, [8] Ra S.K. & Sing R.K. (985), Some estimators for te ratio and product of population parameters. Journal of te Indian Societ of Agricultural Statistics, 37 (), -0. [9] Sa S.M. & Sa D.. (978), Ratio-cum-product estimator for estimating ratio (product) of two population parameters, Sanka, C, 40 (), [0] Sindu S.S., Tailor R. & Sing S. (009), On te estimation of population proportion, Applied Matematical Science, 3(35),

6 American Journal of Operational Researc 05, 5(4): [] Sing G.. and Rani R. (005, 006), Some linear transformations on auxiliar variable for estimating te ratio of two population means in sample surves, Model Assisted Statistics and Applications, (), IOS Press, -5. [] Upadaa.. & Sing H.P. (985), A class of estimators using auxiliar information for estimating ratio of two finite means. Gujarat Statistical Review, (), 7-6. [3] Upadaa.. Sing H.P. & Vos, J.W.E. (985), On te estimation of population means and ratios using supplementar, Statistica. eerlandica, 39(3),

Estimation of Population Mean on Recent Occasion under Non-Response in h-occasion Successive Sampling

Estimation of Population Mean on Recent Occasion under Non-Response in h-occasion Successive Sampling Journal of Modern Applied Statistical Metods Volume 5 Issue Article --06 Estimation of Population Mean on Recent Occasion under Non-Response in -Occasion Successive Sampling Anup Kumar Sarma Indian Scool

More information

VARIANCE ESTIMATION FOR COMBINED RATIO ESTIMATOR

VARIANCE ESTIMATION FOR COMBINED RATIO ESTIMATOR Sankyā : Te Indian Journal of Statistics 1995, Volume 57, Series B, Pt. 1, pp. 85-92 VARIANCE ESTIMATION FOR COMBINED RATIO ESTIMATOR By SANJAY KUMAR SAXENA Central Soil and Water Conservation Researc

More information

Handling Missing Data on Asymmetric Distribution

Handling Missing Data on Asymmetric Distribution International Matematical Forum, Vol. 8, 03, no. 4, 53-65 Handling Missing Data on Asymmetric Distribution Amad M. H. Al-Kazale Department of Matematics, Faculty of Science Al-albayt University, Al-Mafraq-Jordan

More information

Variance Estimation in Stratified Random Sampling in the Presence of Two Auxiliary Random Variables

Variance Estimation in Stratified Random Sampling in the Presence of Two Auxiliary Random Variables International Journal of Science and Researc (IJSR) ISSN (Online): 39-7064 Impact Factor (0): 3.358 Variance Estimation in Stratified Random Sampling in te Presence of Two Auxiliary Random Variables Esubalew

More information

ESTIMATION OF A POPULATION MEAN OF A SENSITIVE VARIABLE IN STRATIFIED TWO-PHASE SAMPLING

ESTIMATION OF A POPULATION MEAN OF A SENSITIVE VARIABLE IN STRATIFIED TWO-PHASE SAMPLING Pak. J. Stati. 06 Vol. 3(5), 393-404 ESTIMATION OF A POPUATION MEAN OF A SENSITIVE VARIABE IN STRATIFIED TWO-PHASE SAMPING Nadia Muaq, Muammad Noor-ul-Amin and Muammad Hanif National College of Business

More information

Department of Statistics & Operations Research, Aligarh Muslim University, Aligarh, India

Department of Statistics & Operations Research, Aligarh Muslim University, Aligarh, India Open Journal of Optimization, 04, 3, 68-78 Publised Online December 04 in SciRes. ttp://www.scirp.org/ournal/oop ttp://dx.doi.org/0.436/oop.04.34007 Compromise Allocation for Combined Ratio Estimates of

More information

Improved Rotated Finite Difference Method for Solving Fractional Elliptic Partial Differential Equations

Improved Rotated Finite Difference Method for Solving Fractional Elliptic Partial Differential Equations American Scientific Researc Journal for Engineering, Tecnolog, and Sciences (ASRJETS) ISSN (Print) 33-44, ISSN (Online) 33-44 Global Societ of Scientific Researc and Researcers ttp://asrjetsjournal.org/

More information

Chapter-2: A Generalized Ratio and Product Type Estimator for the Population Mean in Stratified Random Sampling CHAPTER-2

Chapter-2: A Generalized Ratio and Product Type Estimator for the Population Mean in Stratified Random Sampling CHAPTER-2 Capter-: A Generalized Ratio and Product Tpe Eimator for te Population Mean in tratified Random ampling CHAPTER- A GEERALIZED RATIO AD PRODUCT TYPE ETIMATOR FOR THE POPULATIO MEA I TRATIFIED RADOM AMPLIG.

More information

EFFICIENT REPLICATION VARIANCE ESTIMATION FOR TWO-PHASE SAMPLING

EFFICIENT REPLICATION VARIANCE ESTIMATION FOR TWO-PHASE SAMPLING Statistica Sinica 13(2003), 641-653 EFFICIENT REPLICATION VARIANCE ESTIMATION FOR TWO-PHASE SAMPLING J. K. Kim and R. R. Sitter Hankuk University of Foreign Studies and Simon Fraser University Abstract:

More information

How to Find the Derivative of a Function: Calculus 1

How to Find the Derivative of a Function: Calculus 1 Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te

More information

Exponential Ratio Type Estimators In Stratified Random Sampling

Exponential Ratio Type Estimators In Stratified Random Sampling Eponential atio Tpe Eimators In tratified andom ampling ajes ing, Mukes Kumar,. D. ing, M. K. Caudar Department of tatiics, B.H.U., Varanasi (U.P.-India Corresponding autor Abract Kadilar and Cingi (003

More information

On the Concept of Returns to Scale: Revisited

On the Concept of Returns to Scale: Revisited 3 J. Asian Dev. Stud, Vol. 5, Issue, (Marc 206) ISSN 2304-375X On te Concept of Returns to Scale: Revisited Parvez Azim Abstract Tis paper sows w it is tat in Economics text books and literature we invariabl

More information

De-Coupler Design for an Interacting Tanks System

De-Coupler Design for an Interacting Tanks System IOSR Journal of Electrical and Electronics Engineering (IOSR-JEEE) e-issn: 2278-1676,p-ISSN: 2320-3331, Volume 7, Issue 3 (Sep. - Oct. 2013), PP 77-81 De-Coupler Design for an Interacting Tanks System

More information

Polynomial Interpolation

Polynomial Interpolation Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc

More information

Ratio estimation using stratified ranked set sample

Ratio estimation using stratified ranked set sample METRON - International Journal of Statistics 003, vol. LXI, n. 1, pp. 75-90 HANI M. SAMAWI MAHMOUD I. SIAM Ratio estimation using stratified ranked set sample Summary - Ratio estimation metod is used to

More information

EFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS

EFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS Statistica Sinica 24 2014, 395-414 doi:ttp://dx.doi.org/10.5705/ss.2012.064 EFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS Jun Sao 1,2 and Seng Wang 3 1 East Cina Normal University,

More information

Model Selection in Functional Networks via Genetic Algorithms

Model Selection in Functional Networks via Genetic Algorithms Model Selection in Functional Networs via Genetic Algoritms R E Pruneda and B Lacruz Abstract Several statistical tools and most recentl Functional Networs (FN) ave been used to solve nonlinear regression

More information

June : 2016 (CBCS) Body. Load

June : 2016 (CBCS) Body. Load Engineering Mecanics st Semester : Common to all rances Note : Max. marks : 6 (i) ttempt an five questions (ii) ll questions carr equal marks. (iii) nswer sould be precise and to te point onl (iv) ssume

More information

Artificial Neural Network Model Based Estimation of Finite Population Total

Artificial Neural Network Model Based Estimation of Finite Population Total International Journal of Science and Researc (IJSR), India Online ISSN: 2319-7064 Artificial Neural Network Model Based Estimation of Finite Population Total Robert Kasisi 1, Romanus O. Odiambo 2, Antony

More information

arxiv: v1 [math.na] 11 Apr 2016

arxiv: v1 [math.na] 11 Apr 2016 Hig order approximation to non-smoot multivariate functions Anat Amir David Levin arxiv:164.281v1 [mat.na] 11 Apr 216 April 12, 216 Abstract Approximations of non-smoot multivariate functions return low-order

More information

ESCUELA LATINOAMERICANA DE COOPERACIÓN Y DESARROLLO Especialización en Cooperación Internacional para el Desarrollo

ESCUELA LATINOAMERICANA DE COOPERACIÓN Y DESARROLLO Especialización en Cooperación Internacional para el Desarrollo ESCUELA LATINOAMERICANA DE COOPERACIÓN Y DESARROLLO Especialización en Cooperación Internacional para el Desarrollo A SURVIVAL KIT IN CASE O.MATHEMATICS b Marco Missaglia * Universit of Pavia September

More information

Numerical Differentiation

Numerical Differentiation Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function

More information

An Accurate Self-Starting Initial Value Solvers for. Second Order Ordinary Differential Equations

An Accurate Self-Starting Initial Value Solvers for. Second Order Ordinary Differential Equations International Journal of Contemporar Matematical Sciences Vol. 9, 04, no. 5, 77-76 HIKARI Ltd, www.m-iari.com ttp://dx.doi.org/0.988/icms.04.4554 An Accurate Self-Starting Initial Value Solvers for Second

More information

lecture 26: Richardson extrapolation

lecture 26: Richardson extrapolation 43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)

More information

REVIEW LAB ANSWER KEY

REVIEW LAB ANSWER KEY REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g

More information

Investigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001

Investigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001 Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x

More information

Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1

Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1 Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.

More information

Optimal Search of Developed Class of Modified Ratio Estimators for Estimation of Population Variance

Optimal Search of Developed Class of Modified Ratio Estimators for Estimation of Population Variance American Journal of Operational Research 05 5(4): 8-95 DOI: 0.593/j.ajor.050504.0 Optimal earch of Developed Class of Modified Ratio Estimators for Estimation of Population Variance Alok Kumar hukla heela

More information

Overdispersed Variational Autoencoders

Overdispersed Variational Autoencoders Overdispersed Variational Autoencoders Harsil Sa, David Barber and Aleksandar Botev Department of Computer Science, University College London Alan Turing Institute arsil.sa.15@ucl.ac.uk, david.barber@ucl.ac.uk,

More information

ch (for some fixed positive number c) reaching c

ch (for some fixed positive number c) reaching c GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan

More information

Lecture 21. Numerical differentiation. f ( x+h) f ( x) h h

Lecture 21. Numerical differentiation. f ( x+h) f ( x) h h Lecture Numerical differentiation Introduction We can analytically calculate te derivative of any elementary function, so tere migt seem to be no motivation for calculating derivatives numerically. However

More information

Copyright c 2008 Kevin Long

Copyright c 2008 Kevin Long Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula

More information

232 Calculus and Structures

232 Calculus and Structures 3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE

More information

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of

More information

Polynomial Interpolation

Polynomial Interpolation Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x

More information

2.11 That s So Derivative

2.11 That s So Derivative 2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point

More information

Desalination by vacuum membrane distillation: sensitivity analysis

Desalination by vacuum membrane distillation: sensitivity analysis Separation and Purification Tecnology 33 (2003) 75/87 www.elsevier.com/locate/seppur Desalination by vacuum membrane distillation: sensitivity analysis Fawzi Banat *, Fami Abu Al-Rub, Kalid Bani-Melem

More information

Fractional Derivatives as Binomial Limits

Fractional Derivatives as Binomial Limits Fractional Derivatives as Binomial Limits Researc Question: Can te limit form of te iger-order derivative be extended to fractional orders? (atematics) Word Count: 669 words Contents - IRODUCIO... Error!

More information

A = h w (1) Error Analysis Physics 141

A = h w (1) Error Analysis Physics 141 Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.

More information

lecture 35: Linear Multistep Mehods: Truncation Error

lecture 35: Linear Multistep Mehods: Truncation Error 88 lecture 5: Linear Multistep Meods: Truncation Error 5.5 Linear ultistep etods One-step etods construct an approxiate solution x k+ x(t k+ ) using only one previous approxiation, x k. Tis approac enoys

More information

Definition of the Derivative

Definition of the Derivative Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of

More information

Estimating Peak Bone Mineral Density in Osteoporosis Diagnosis by Maximum Distribution

Estimating Peak Bone Mineral Density in Osteoporosis Diagnosis by Maximum Distribution International Journal of Clinical Medicine Researc 2016; 3(5): 76-80 ttp://www.aascit.org/journal/ijcmr ISSN: 2375-3838 Estimating Peak Bone Mineral Density in Osteoporosis Diagnosis by Maximum Distribution

More information

Financial Econometrics Prof. Massimo Guidolin

Financial Econometrics Prof. Massimo Guidolin CLEFIN A.A. 2010/2011 Financial Econometrics Prof. Massimo Guidolin A Quick Review of Basic Estimation Metods 1. Were te OLS World Ends... Consider two time series 1: = { 1 2 } and 1: = { 1 2 }. At tis

More information

The derivative function

The derivative function Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative

More information

Fuzzy Geometric Programming in Multivariate Stratified Sample Surveys in Presence of Non-Response with Quadratic Cost Function

Fuzzy Geometric Programming in Multivariate Stratified Sample Surveys in Presence of Non-Response with Quadratic Cost Function American Journal of Operations Researc, 04, 4, 73-88 Publised Online May 04 in SciRes. ttp://.scirp.org/ournal/aor ttp://dx.doi.org/0.436/aor.04.4307 Fuzzy Geometric Programming in Multivariate Stratified

More information

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note

More information

These errors are made from replacing an infinite process by finite one.

These errors are made from replacing an infinite process by finite one. Introduction :- Tis course examines problems tat can be solved by metods of approximation, tecniques we call numerical metods. We begin by considering some of te matematical and computational topics tat

More information

On the modified Reynolds equation for journal bearings in a case of non-newtonian Rabinowitsch fluid model

On the modified Reynolds equation for journal bearings in a case of non-newtonian Rabinowitsch fluid model MATEC Web of Conferences 45, 7 (8) NCTAM 7 ttps://doi.org/.5/matecconf/8457 On te modified Renolds equation for journal bearings in a case of non-newtonian Rabinowitsc fluid model Juliana Javorova,*, and

More information

Test 2 Review. 1. Find the determinant of the matrix below using (a) cofactor expansion and (b) row reduction. A = 3 2 =

Test 2 Review. 1. Find the determinant of the matrix below using (a) cofactor expansion and (b) row reduction. A = 3 2 = Test Review Find te determinant of te matrix below using (a cofactor expansion and (b row reduction Answer: (a det + = (b Observe R R R R R R R R R Ten det B = (((det Hence det Use Cramer s rule to solve:

More information

Strati cation by Size Revisited

Strati cation by Size Revisited Journal of Of cial Statistics, Vol. 16, No. 2, 2000, pp. 139±154 Strati cation by Size Revisited Alan H. Dorfman 1 and Ricard Valliant 2 Strati cation by size is used in nite population sampling as a means

More information

WYSE Academic Challenge 2004 State Finals Mathematics Solution Set

WYSE Academic Challenge 2004 State Finals Mathematics Solution Set WYSE Academic Callenge 00 State Finals Matematics Solution Set. Answer: c. We ave a sstem of tree equations and tree unknowns. We ave te equations: x + + z 0, x + 6 + 7z 9600, and 7x + + z 90. Wen we solve,

More information

Differentiation. Area of study Unit 2 Calculus

Differentiation. Area of study Unit 2 Calculus Differentiation 8VCE VCEco Area of stud Unit Calculus coverage In tis ca 8A 8B 8C 8D 8E 8F capter Introduction to limits Limits of discontinuous, rational and brid functions Differentiation using first

More information

The total error in numerical differentiation

The total error in numerical differentiation AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and

More information

The Priestley-Chao Estimator

The Priestley-Chao Estimator Te Priestley-Cao Estimator In tis section we will consider te Pristley-Cao estimator of te unknown regression function. It is assumed tat we ave a sample of observations (Y i, x i ), i = 1,..., n wic are

More information

Fundamentals of Concept Learning

Fundamentals of Concept Learning Aims 09s: COMP947 Macine Learning and Data Mining Fundamentals of Concept Learning Marc, 009 Acknowledgement: Material derived from slides for te book Macine Learning, Tom Mitcell, McGraw-Hill, 997 ttp://www-.cs.cmu.edu/~tom/mlbook.tml

More information

Chapter. Differentiation: Basic Concepts. 1. The Derivative: Slope and Rates. 2. Techniques of Differentiation. 3. The Product and Quotient Rules

Chapter. Differentiation: Basic Concepts. 1. The Derivative: Slope and Rates. 2. Techniques of Differentiation. 3. The Product and Quotient Rules Differentiation: Basic Concepts Capter 1. Te Derivative: Slope and Rates 2. Tecniques of Differentiation 3. Te Product and Quotient Rules 4. Marginal Analsis: Approimation b Increments 5. Te Cain Rule

More information

A Multiaxial Variable Amplitude Fatigue Life Prediction Method Based on a Plane Per Plane Damage Assessment

A Multiaxial Variable Amplitude Fatigue Life Prediction Method Based on a Plane Per Plane Damage Assessment American Journal of Mecanical and Industrial Engineering 28; 3(4): 47-54 ttp://www.sciencepublisinggroup.com/j/ajmie doi:.648/j.ajmie.2834.2 ISSN: 2575-679 (Print); ISSN: 2575-66 (Online) A Multiaxial

More information

The effects of shear stress on the lubrication performances of oil film of large-scale mill bearing

The effects of shear stress on the lubrication performances of oil film of large-scale mill bearing Universit of Wollongong Researc Online Facult of Engineering - Papers (Arcive) Facult of Engineering and Information Sciences 9 Te effects of sear stress on te lubrication performances of oil film of large-scale

More information

School of Geomatics and Urban Information, Beijing University of Civil Engineering and Architecture, Beijing, China 2

School of Geomatics and Urban Information, Beijing University of Civil Engineering and Architecture, Beijing, China 2 Examination Metod and Implementation for Field Survey Data of Crop Types Based on Multi-resolution Satellite Images Yang Liu, Mingyi Du, Wenquan Zu, Scool of Geomatics and Urban Information, Beijing University

More information

2.8 The Derivative as a Function

2.8 The Derivative as a Function .8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open

More information

Symmetry Labeling of Molecular Energies

Symmetry Labeling of Molecular Energies Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry

More information

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225 THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:

More information

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department

More information

Parametric Spline Method for Solving Bratu s Problem

Parametric Spline Method for Solving Bratu s Problem ISSN 749-3889 print, 749-3897 online International Journal of Nonlinear Science Vol4202 No,pp3-0 Parametric Spline Metod for Solving Bratu s Problem M Zarebnia, Z Sarvari 2,2 Department of Matematics,

More information

(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a?

(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a? Solutions to Test 1 Fall 016 1pt 1. Te grap of a function f(x) is sown at rigt below. Part I. State te value of eac limit. If a limit is infinite, state weter it is or. If a limit does not exist (but is

More information

Lines, Conics, Tangents, Limits and the Derivative

Lines, Conics, Tangents, Limits and the Derivative Lines, Conics, Tangents, Limits and te Derivative Te Straigt Line An two points on te (,) plane wen joined form a line segment. If te line segment is etended beond te two points ten it is called a straigt

More information

Function Composition and Chain Rules

Function Composition and Chain Rules Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function

More information

Solutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014

Solutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014 Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.

More information

1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow

1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow 1.7, Groundwater Hydrology Prof. Carles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow Simulation: Te prediction of quantities of interest (dependent variables) based upon an equation

More information

Two Step Hybrid Block Method with Two Generalized Off-step Points for Solving Second Ordinary Order Differential Equations Directly

Two Step Hybrid Block Method with Two Generalized Off-step Points for Solving Second Ordinary Order Differential Equations Directly Global Journal of Pure and Applied Matematics. ISSN 0973-768 Volume 2, Number 2 (206), pp. 59-535 Researc India Publications ttp://www.ripublication.com/gjpam.tm Two Step Hybrid Block Metod wit Two Generalized

More information

University Mathematics 2

University Mathematics 2 University Matematics 2 1 Differentiability In tis section, we discuss te differentiability of functions. Definition 1.1 Differentiable function). Let f) be a function. We say tat f is differentiable at

More information

Diffraction. S.M.Lea. Fall 1998

Diffraction. S.M.Lea. Fall 1998 Diffraction.M.Lea Fall 1998 Diffraction occurs wen EM waves approac an aperture (or an obstacle) wit dimension d > λ. We sall refer to te region containing te source of te waves as region I and te region

More information

3.1 Extreme Values of a Function

3.1 Extreme Values of a Function .1 Etreme Values of a Function Section.1 Notes Page 1 One application of te derivative is finding minimum and maimum values off a grap. In precalculus we were only able to do tis wit quadratics by find

More information

Wind Turbine Micrositing: Comparison of Finite Difference Method and Computational Fluid Dynamics

Wind Turbine Micrositing: Comparison of Finite Difference Method and Computational Fluid Dynamics IJCSI International Journal of Computer Science Issues, Vol. 9, Issue 1, No 1, January 01 ISSN (Online): 169-081 www.ijcsi.org 7 Wind Turbine Micrositing: Comparison of Finite Difference Metod and Computational

More information

Combining functions: algebraic methods

Combining functions: algebraic methods Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)

More information

Lecture 10: Carnot theorem

Lecture 10: Carnot theorem ecture 0: Carnot teorem Feb 7, 005 Equivalence of Kelvin and Clausius formulations ast time we learned tat te Second aw can be formulated in two ways. e Kelvin formulation: No process is possible wose

More information

Tail Conditional Expectations for Extended Exponential Dispersion Models

Tail Conditional Expectations for Extended Exponential Dispersion Models American Researc Journal of Matematics Original Article ISSN 378-704 Volume 1 Issue 4 015 Tail Conditional Expectations for Extended Exponential Dispersion Models Ye (Zoe) Ye Qiang Wu and Don Hong 1 Program

More information

Continuity and Differentiability Worksheet

Continuity and Differentiability Worksheet Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;

More information

3.4 Worksheet: Proof of the Chain Rule NAME

3.4 Worksheet: Proof of the Chain Rule NAME Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are

More information

Name: Sept 21, 2017 Page 1 of 1

Name: Sept 21, 2017 Page 1 of 1 MATH 111 07 (Kunkle), Eam 1 100 pts, 75 minutes No notes, books, electronic devices, or outside materials of an kind. Read eac problem carefull and simplif our answers. Name: Sept 21, 2017 Page 1 of 1

More information

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator. Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative

More information

Effect of the Dependent Paths in Linear Hull

Effect of the Dependent Paths in Linear Hull 1 Effect of te Dependent Pats in Linear Hull Zenli Dai, Meiqin Wang, Yue Sun Scool of Matematics, Sandong University, Jinan, 250100, Cina Key Laboratory of Cryptologic Tecnology and Information Security,

More information

NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE

NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE Global and Stocastic Analysis Vol. 4 No. 1, January 2017, 1-10 NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE K. PHANEENDRA AND E. SIVA PRASAD Abstract.

More information

Click here to see an animation of the derivative

Click here to see an animation of the derivative Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,

More information

Application of Quintic B-splines Collocation Method on Some Rosenau Type Nonlinear Higher Order Evolution Equations.

Application of Quintic B-splines Collocation Method on Some Rosenau Type Nonlinear Higher Order Evolution Equations. ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.13(2012) No.2,pp.142-152 Application of Quintic B-splines Collocation Metod on Some Rosenau Type Nonlinear Higer

More information

A STATIC PDE APPROACH FOR MULTI-DIMENSIONAL EXTRAPOLATION USING FAST SWEEPING METHODS

A STATIC PDE APPROACH FOR MULTI-DIMENSIONAL EXTRAPOLATION USING FAST SWEEPING METHODS A STATIC PDE APPROACH FOR MULTI-DIMENSIONAL EXTRAPOLATION USING FAST SWEEPING METHODS TARIQ ASLAM, SONGTING LUO, AND HONGKAI ZHAO Abstract. A static Partial Differential Equation (PDE) approac is presented

More information

Optimal parameters for a hierarchical grid data structure for contact detection in arbitrarily polydisperse particle systems

Optimal parameters for a hierarchical grid data structure for contact detection in arbitrarily polydisperse particle systems Comp. Part. Mec. 04) :357 37 DOI 0.007/s4057-04-000-9 Optimal parameters for a ierarcical grid data structure for contact detection in arbitrarily polydisperse particle systems Dinant Krijgsman Vitaliy

More information

A finite element approximation for the quasi-static Maxwell Landau Lifshitz Gilbert equations

A finite element approximation for the quasi-static Maxwell Landau Lifshitz Gilbert equations ANZIAM J. 54 (CTAC2012) pp.c681 C698, 2013 C681 A finite element approximation for te quasi-static Maxwell Landau Lifsitz Gilbert equations Kim-Ngan Le 1 T. Tran 2 (Received 31 October 2012; revised 29

More information

Inf sup testing of upwind methods

Inf sup testing of upwind methods INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Met. Engng 000; 48:745 760 Inf sup testing of upwind metods Klaus-Jurgen Bate 1; ;, Dena Hendriana 1, Franco Brezzi and Giancarlo

More information

Grid-independent large-eddy simulation of compressible turbulent flows using explicit filtering

Grid-independent large-eddy simulation of compressible turbulent flows using explicit filtering Center for Turbulence Researc Proceedings of te Summer Program 2010 203 Grid-independent large-edd simulation of compressible turbulent flows using explicit filtering B D. You, S. T. Bose AND P. Moin Te

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat

More information

The cluster problem in constrained global optimization

The cluster problem in constrained global optimization Te cluster problem in constrained global optimization Te MIT Faculty as made tis article openly available. Please sare ow tis access benefits you. Your story matters. Citation As Publised Publiser Kannan,

More information

Kernel Density Estimation

Kernel Density Estimation Kernel Density Estimation Univariate Density Estimation Suppose tat we ave a random sample of data X 1,..., X n from an unknown continuous distribution wit probability density function (pdf) f(x) and cumulative

More information

Section 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is

Section 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit

More information

Simulation and verification of a plate heat exchanger with a built-in tap water accumulator

Simulation and verification of a plate heat exchanger with a built-in tap water accumulator Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation

More information

Department of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801

Department of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801 RESEARCH SUMMARY AND PERSPECTIVES KOFFI B. FADIMBA Department of Matematical Sciences University of Sout Carolina Aiken Aiken, SC 29801 Email: KoffiF@usca.edu 1. Introduction My researc program as focused

More information

1 ode.mcd. Find solution to ODE dy/dx=f(x,y). Instructor: Nam Sun Wang

1 ode.mcd. Find solution to ODE dy/dx=f(x,y). Instructor: Nam Sun Wang Fin solution to ODE /=f(). Instructor: Nam Sun Wang oe.mc Backgroun. Wen a sstem canges wit time or wit location, a set of ifferential equations tat contains erivative terms "/" escribe suc a namic sstem.

More information

Higher Derivatives. Differentiable Functions

Higher Derivatives. Differentiable Functions Calculus 1 Lia Vas Higer Derivatives. Differentiable Functions Te second derivative. Te derivative itself can be considered as a function. Te instantaneous rate of cange of tis function is te second derivative.

More information

THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA EXAMINATION MODULE 5

THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA EXAMINATION MODULE 5 THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA EXAMINATION NEW MODULAR SCHEME introduced from te examinations in 009 MODULE 5 SOLUTIONS FOR SPECIMEN PAPER B THE QUESTIONS ARE CONTAINED IN A SEPARATE FILE

More information

Dedicated to the 70th birthday of Professor Lin Qun

Dedicated to the 70th birthday of Professor Lin Qun Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,

More information