18 Numerical Integration of Functions
|
|
- Clement Strickland
- 5 years ago
- Views:
Transcription
1 Slightly moifie //9, /8/6 Firstly written at Marh 5 8 Numerial Integration of Funtions Introution Romberg Integration Gauss Quarature Aaptive Quarature Case Stuy: Root-Mean-Square Current DM869/Computational Numerial Analysis/7_na.o 7
2 Slightly moifie //9, /8/6 Firstly written at Marh 5 Introution two forms of numerial integration table of values funtion Riharson extrapolation - ombining two numerial integral estimates to obtain a thir, more aurate value - an effiient omputational algorithm for Riharson's extrapolation Romberg integration Gauss quarature aaptive quarature Romberg Integration one tehnique esigne to attain effiient numerial integrals of funtions base on suessive appliation of the trapezoial rule Riharson's extrapolation using two estimates of an integral to ompute a thir, more aurate approximation Let us onsier: I I ( h E( h where I = the exat value of the integral I(h = the approximation from an n-segment appliation of the trapezoial rule h = (b a/n, step size E(h = the trunation error Making two separate estimates an having exat values for the error: I( h E( h I( h E( h Reall that b a E h f an assume that f is onstant regarless of step size. Then we an have: E( h h or E( h h whih yiels: E ( h h I( h E( h I( h E( h h I( h I( h or E( h ( h / h E( h 7 DM869/Computational Numerial Analysis/7_na.o h h
3 Slightly moifie //9, /8/6 Firstly written at Marh 5 This estimate of the truation error an be substitute into I I( h E( h to yiel an improve estimate of the integral: I I( h [ I ( h I( h ] error ~ O(h ( h / h for the speial ase where h = h /: I I( h I( h Example Use Riharson extrapolation to evaluate the integral of f ( x. 5x x 675x 9x x from a = to b =.8. Note that the exat value is Sol. Single an omposite trapezoial rule: Segments h Integral t %.9% 9.5% for & segments: I (.688 ( E t = =.767 ( t = 6.6% for & segments: I (.88 ( E t = =.767 ( t =.% Riharson extrapolation again for two improve integrals 6 I I m I l O(h an again Riharson extrapolation I I m I l O(h where I m = the more aurate estimate I l = the less aurate estimate DM869/Computational Numerial Analysis/7_na.o 75
4 Slightly moifie //9, /8/6 Firstly written at Marh 5 Example Repeat the above example to ompute an integral with O(h 6. Sol. 6 I (.667 ( the exat value of the integral The Romberg integration algorithm general form well suite for omputer implementation I j, k k I j, k k I j, k where I j+,k- = the more aurate integral I j,k- = the less aurate integral I j,k = the improve integral k signifying the level of the integration e.g., k = the original trapezoial rule estimates k = the O(h estimates k = the O(h 6 estimates I, I, ex for k = an j = : I, I( h I( h Graphial epition of the sequene of integral estimates generate using Romberg integration. (a First iteration. (b Seon iteration. ( Thir iteration. 76 DM869/Computational Numerial Analysis/7_na.o
5 Slightly moifie //9, /8/6 Firstly written at Marh 5 termination (or stopping riterion - an estimate of the perent relative error ε a I, k I I, k, k % Romberg integration is more effiient than the trapezoial rule an Simpson's rules - Simpson's rule requires a 8-segment appliation for the above example - but, Romberg integration only requires 5 funtion evaluation <An M-file to implement Romberg integration> DM869/Computational Numerial Analysis/7_na.o 77
6 Slightly moifie //9, /8/6 Firstly written at Marh 5 Gauss Quarature (a Graphial epition of the trapezoial rule as the area uner the straight line joining fixe en points. (b An improve integral estimate obtaine by taking the area uner the straight line passing through two intermeiate points. By positioning these points wisely, the positive an negative errors are better balane, an an improve integral estimate results. the trapezoial rule preetermine or fixe points for the integration interval sometimes resulting in a large error if taking the enpoints freely an wisely an improve estimate by balaning the positive an negative errors Metho of unetermine oeffiients Reall the trapezoial rule: Two integrals that shoul be evaluate exatly by the trapezoial rule: (a a onstant an (b a straight line. f ( a f ( b I ( b a or I f a f ( ( b The trapezoial rule yiels exat results when the funtion being integrate is a onstant or a straight line. Thus, an ( b a / ( ba / x b a b a b a ( b a / b a b a xx ( ba / 78 DM869/Computational Numerial Analysis/7_na.o
7 Slightly moifie //9, /8/6 Firstly written at Marh 5 Solving for the two equations with two unknowns: b a Finally, we have b a b a I f ( a f ( b equivalent to the trapezoial rule Derivation of the two-point Gauss-Legenre formula the objet of Gauss quarature - to etermine the oeffiients of an equation of the form I f ( x f ( x where the 's = the unknown oeffiients - note that not fixe x an x four unknowns requiring four onitions y = onstant, y = x, y = x, an y = x x x x xx x x x x x x x x x x = x sine x x x x.5775 Graphial epition of the unknown variables x an x for integration by Gauss quarature. x.5775 Therefore, the two-point Gauss-Legenre formula is I f ( f ( DM869/Computational Numerial Analysis/7_na.o 79
8 Slightly moifie //9, /8/6 Firstly written at Marh 5 general formulation in the integral limits assuming that x a x a - the lower limit: x = a then x = a a ( a - the upper limit: x = b then x = b a ( a Solving: a a b a b a Then we have: ( b a ( b a x x b a x x Example Use the two-point Gauss-Legenre formula to evaluate the integral of f ( x. 5x x 675x 9x x between the limits x = to.8. Note that the exat value is Sol. Changing the limits to an +: x x.. an x. x.8 (. 5x x - 675x [. 5(..x (..x 9(..x 9x (..x 5 x x 675(..x 5 ].x.567 at x / an.587 at x / =.8578 a perent relative error of.% DM869/Computational Numerial Analysis/7_na.o
9 Slightly moifie //9, /8/6 Firstly written at Marh 5 Higher-point formulas I f ( x f ( x n f ( xn where n = the number of points <Weighting fators an funtion arguments use in Gauss-Legenre formulas.> Points Weighting fators Funtion arguments Trunation error 5 6 =. =. = = = =.7858 =.655 =.655 =.7858 =.6969 = = = =.6969 =.75 =.6766 =.6799 =.6799 = =.75 x = x = x = x =. x = x =.866 x =.998 x =.998 x =.866 x = x =.5869 x =. x =.5869 x = x =.9695 x = x = x = x = x 5 =.9695 f ( ( f (6 ( f (8 ( f ( ( f ( ( Example Use the three-point Gauss-Legenre formula to evaluate the integral of f ( x. 5x x 675x 9x x between the limits x = to.8. Note that the exat value is Sol. I f ( f ( f ( I not appropriate when the funtion is unknown not suite for engineering problem with tabulate ata DM869/Computational Numerial Analysis/7_na.o 8
10 Slightly moifie //9, /8/6 Firstly written at Marh 5 Aaptive Quarature omposite Simpson's / rule is applie for equally spae points not take into aount relatively abruptly hange regions where more refine spaing might be require aaptive quarature methos - automatially ajusting the step size - small steps in regions of sharp variations - larger steps in regions of graual variations - applying the / rule at two levels of refinement - estimating the trunation errors for the two levels - refining the step size an repeating until the error falls to aeptable levels, if the error is too large MATLAB funtions: qua an qual qua - aaptive Simpson quarature - more effiient for low auaries or nonsmooth funtions qual - Lobatto quarature - more effiient for high auraies an smooth funtions q = qua(fun, a, b, tol, trae, p, p,... where fun the funtion to be integrate a an b the integration bouns tol the esire absolute error tolerane (efault = -6 trae the optional variable p, p,... parameters passe to fun 8 DM869/Computational Numerial Analysis/7_na.o
11 Slightly moifie //9, /8/6 Firstly written at Marh 5 Example Use qua to integrate the following funtion: f ( x s ( x q. ( x r. between the limits x = to. Note that for q =., r =.9, an s = 6, this is the built-in humps funtion that MATLAB uses to emonstrate some of its numerial apabilities. The humps funtion exhibits both flat an steep regions over a relatively short x range. Hene, it is useful for emonstrating an testing routines like qua an qual. Note that the humps funtion an be integrate analytially between the given limits to yiel an exat integral of Sol. Using built-in version of humps along with the efault tolerane: >> format long >> qua(@humps,, ans = Using a looser tolerane: funtion y = myhumps(x,q,r,s y =./((x-q.^ +. +./((x-r.^ +. s; >> qua(@myhumps,,,e-,[],.,.9,6 ans = DM869/Computational Numerial Analysis/7_na.o 8
12 8 DM869/Computational Numerial Analysis/7_na.o HK Kim Slightly moifie //9, /8/6 Firstly written at Marh 5
Numerical Methods School of Mechanical Engineering Chung-Ang University
Part 5 Chapter 8 Numerical Integration of Functions Prof. Hae-Jin Choi hjchoi@cau.ac.kr Numerical Methos 00- Chapter Objectives l Unerstaning how Richarson extrapolation provies a means to create a more
More informationNumerical Integration of Functions
Numerical Integration of Functions Berlin Chen Department of Computer Science & Information Engineering National Taiwan Normal University Reference: 1. Applied Numerical Methods with MATLAB for Engineers,
More informationChapter 2: One-dimensional Steady State Conduction
1 Chapter : One-imensional Steay State Conution.1 Eamples of One-imensional Conution Eample.1: Plate with Energy Generation an Variable Conutivity Sine k is variable it must remain insie the ifferentiation
More informationMath 225B: Differential Geometry, Homework 6
ath 225B: Differential Geometry, Homework 6 Ian Coley February 13, 214 Problem 8.7. Let ω be a 1-form on a manifol. Suppose that ω = for every lose urve in. Show that ω is exat. We laim that this onition
More informationThe Islamic University of Gaza Faculty of Engineering Civil Engineering Department. Numerical Analysis ECIV Chapter 5. Bracketing Methods
The Islami University of Gaza Faulty of Engineering Civil Engineering Department Numerial Analysis ECIV 3306 Chapter 5 Braketing Methods Assoiate Prof. Mazen Abualtayef Civil Engineering Department, The
More informationChapter 5: Numerical Integration and Differentiation
Chapter 5: Numerical Integration and Differentiation PART I: Numerical Integration Newton-Cotes Integration Formulas The idea of Newton-Cotes formulas is to replace a complicated function or tabulated
More informationAn Integer Solution of Fractional Programming Problem
Gen. Math. Notes, Vol. 4, No., June 0, pp. -9 ISSN 9-784; Copyright ICSRS Publiation, 0 www.i-srs.org Available free online at http://www.geman.in An Integer Solution of Frational Programming Problem S.C.
More informationProject 10.5C Circular Membrane Vibrations
Projet 1.5C Cirular Membrane Vibrations In problems involving regions that enjoy irular symmetry about the origin in the plane (or the vertial z-axis in spae), the use of polar (or ylindrial) oordinates
More informationReliability Optimization With Mixed Continuous-Discrete Random Variables and Parameters
Subroto Gunawan Researh Fellow Panos Y. Papalambros Professor e-mail: pyp@umih.eu Department of Mehanial Engineering, University of Mihigan, Ann Arbor, MI 4809 Reliability Optimization With Mixe Continuous-Disrete
More informationOptimal torque control of permanent magnet synchronous machines using magnetic equivalent circuits
This oument ontains a post-print version of the paper Optimal torque ontrol of permanent magnet synhronous mahines using magneti equivalent iruits authore by W. Kemmetmüller, D. Faustner, an A. Kugi an
More informationChapter 8 Thermodynamic Relations
Chapter 8 Thermodynami Relations 8.1 Types of Thermodynami roperties The thermodynami state of a system an be haraterized by its properties that an be lassified as measured, fundamental, or deried properties.
More informationProcess engineers are often faced with the task of
Fluids and Solids Handling Eliminate Iteration from Flow Problems John D. Barry Middough, In. This artile introdues a novel approah to solving flow and pipe-sizing problems based on two new dimensionless
More informationFast Evaluation of Canonical Oscillatory Integrals
Appl. Math. Inf. Si. 6, No., 45-51 (01) 45 Applie Mathematis & Information Sienes An International Journal 01 NSP Natural Sienes Publishing Cor. Fast Evaluation of Canonial Osillatory Integrals Ying Liu
More informationMcCreight s Suffix Tree Construction Algorithm. Milko Izamski B.Sc. Informatics Instructor: Barbara König
1. Introution MCreight s Suffix Tree Constrution Algorithm Milko Izamski B.S. Informatis Instrutor: Barbara König The main goal of MCreight s algorithm is to buil a suffix tree in linear time. This is
More informationThe optimization of kinematical response of gear transmission
Proeeings of the 7 WSEAS Int. Conferene on Ciruits, Systems, Signal an Teleommuniations, Gol Coast, Australia, January 7-9, 7 The optimization of inematial response of gear transmission VINCENZO NIOLA
More informationTo derive the other Pythagorean Identities, divide the entire equation by + = θ = sin. sinθ cosθ tanθ = 1
Syllabus Objetives: 3.3 The student will simplify trigonometri expressions and prove trigonometri identities (fundamental identities). 3.4 The student will solve trigonometri equations with and without
More informationThe Chebyshev Wavelet Method for Numerical Solutions of A Fractional Oscillator
Shiraz University of Tehnology From the SeleteWorks of Habibolla Latifizaeh 01 The Chebyshev Wavelet Metho for Numerial Solutions of A Frational Osillator E. Hesameini, Shiraz University of Tehnology S.
More informationSome facts you should know that would be convenient when evaluating a limit:
Some fats you should know that would be onvenient when evaluating a it: When evaluating a it of fration of two funtions, f(x) x a g(x) If f and g are both ontinuous inside an open interval that ontains
More informationmax min z i i=1 x j k s.t. j=1 x j j:i T j
AM 221: Advaned Optimization Spring 2016 Prof. Yaron Singer Leture 22 April 18th 1 Overview In this leture, we will study the pipage rounding tehnique whih is a deterministi rounding proedure that an be
More informationEcon 455 Answers - Problem Set Consider a small country (Belgium) with the following demand and supply curves for corn:
Spring 004 Eon 455 Harvey Lapan Eon 455 Answers - Problem Set 4 1. Consier a small ountry (Belgium with the ollowing eman an supply urves or orn: Supply = 4P s ; Deman = 1000 Assume Belgium an import steel
More informationDetermination the Invert Level of a Stilling Basin to Control Hydraulic Jump
Global Avane Researh Journal of Agriultural Siene Vol. (4) pp. 074-079, June, 0 Available online http://garj.org/garjas/inex.htm Copyright 0 Global Avane Researh Journals Full Length Researh Paper Determination
More informationJournal of Mathematical Analysis and Applications
J. Math. Anal. Appl. 371 (010) 759 763 Contents lists available at SieneDiret Journal of Mathematial Analysis an Appliations www.elsevier.om/loate/jmaa Singular Sturm omparison theorems Dov Aharonov, Uri
More informationZero-Free Region for ζ(s) and PNT
Contents Zero-Free Region for ζs an PN att Rosenzweig Chebyshev heory ellin ransforms an Perron s Formula Zero-Free Region of Zeta Funtion 6. Jensen s Inequality..........................................
More informationFINITE WORD LENGTH EFFECTS IN DSP
FINITE WORD LENGTH EFFECTS IN DSP PREPARED BY GUIDED BY Snehal Gor Dr. Srianth T. ABSTRACT We now that omputers store numbers not with infinite preision but rather in some approximation that an be paed
More informationChapter 9. There are 7 out of 50 measurements that are greater than or equal to 5.1; therefore, the fraction of the
Pratie questions 6 1 a y i = 6 µ = = 1 i = 1 y i µ i = 1 ( ) = 95 = s n 95 555. x i f i 1 1+ + 5+ n + 5 5 + n µ = = = f 11+ n 11+ n i 7 + n = 5 + n = 6n n = a Time (minutes) 1.6.1.6.1.6.1.6 5.1 5.6 6.1
More informationForce Reconstruction for Nonlinear Structures in Time Domain
Fore Reonstrution for Nonlinear Strutures in ime Domain Jie Liu 1, Bing Li 2, Meng Li 3, an Huihui Miao 4 1,2,3,4 State Key Laboratory for Manufaturing Systems Engineering, Xi an Jiaotong niversity, Xi
More informationWavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013
Ultrafast Pulses and GVD John O Hara Created: De. 6, 3 Introdution This doument overs the basi onepts of group veloity dispersion (GVD) and ultrafast pulse propagation in an optial fiber. Neessarily, it
More informationSingular Event Detection
Singular Event Detetion Rafael S. Garía Eletrial Engineering University of Puerto Rio at Mayagüez Rafael.Garia@ee.uprm.edu Faulty Mentor: S. Shankar Sastry Researh Supervisor: Jonathan Sprinkle Graduate
More informationDifferential Equations 8/24/2010
Differential Equations A Differential i Equation (DE) is an equation ontaining one or more derivatives of an unknown dependant d variable with respet to (wrt) one or more independent variables. Solution
More informationREFINED UPPER BOUNDS FOR THE LINEAR DIOPHANTINE PROBLEM OF FROBENIUS. 1. Introduction
Version of 5/2/2003 To appear in Advanes in Applied Mathematis REFINED UPPER BOUNDS FOR THE LINEAR DIOPHANTINE PROBLEM OF FROBENIUS MATTHIAS BECK AND SHELEMYAHU ZACKS Abstrat We study the Frobenius problem:
More informationA Queueing Model for Call Blending in Call Centers
A Queueing Model for Call Blending in Call Centers Sandjai Bhulai and Ger Koole Vrije Universiteit Amsterdam Faulty of Sienes De Boelelaan 1081a 1081 HV Amsterdam The Netherlands E-mail: {sbhulai, koole}@s.vu.nl
More informationIs the Free Vacuum Energy Infinite?
Is the Free Vauum Energy Infite? H. Razmi () an S. M. Shirazi () Department of Physis, the University of Qom, Qom, I. R. Iran. () razmi@qom.a.ir & razmiha@hotmail.om () sms0@gmail.om Abstrat Consierg the
More informationSensitivity Analysis of Resonant Circuits
1 Sensitivity Analysis of Resonant Ciruits Olivier Buu Abstrat We use first-orer perturbation theory to provie a loal linear relation between the iruit parameters an the poles of an RLC network. The sensitivity
More informationHe s Semi-Inverse Method and Ansatz Approach to look for Topological and Non-Topological Solutions Generalized Nonlinear Schrödinger Equation
Quant. Phys. Lett. 3, No. 2, 23-27 2014) 23 Quantum Physis Letters An International Journal http://x.oi.org/10.12785/qpl/030202 He s Semi-Inverse Metho an Ansatz Approah to look for Topologial an Non-Topologial
More informationLabeling Workflow Views with Fine-Grained Dependencies
Labeling Workflow Views with Fine-Graine Depenenies Zhuowei Bao Department of omputer an Information iene University of Pennsylvania Philaelphia, P 1914, U zhuowei@is.upenn.eu usan B. Davison Department
More informationUPPER-TRUNCATED POWER LAW DISTRIBUTIONS
Fratals, Vol. 9, No. (00) 09 World Sientifi Publishing Company UPPER-TRUNCATED POWER LAW DISTRIBUTIONS STEPHEN M. BURROUGHS and SARAH F. TEBBENS College of Marine Siene, University of South Florida, St.
More informationSPLINE ESTIMATION OF SINGLE-INDEX MODELS
SPLINE ESIMAION OF SINGLE-INDEX MODELS Li Wang and Lijian Yang University of Georgia and Mihigan State University Supplementary Material his note ontains proofs for the main results he following two propositions
More informationSolutions to Problem Set 1
Eon602: Maro Theory Eonomis, HKU Instrutor: Dr. Yulei Luo September 208 Solutions to Problem Set. [0 points] Consider the following lifetime optimal onsumption-saving problem: v (a 0 ) max f;a t+ g t t
More informationComplexity of Regularization RBF Networks
Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw
More informationMethods of evaluating tests
Methods of evaluating tests Let X,, 1 Xn be i.i.d. Bernoulli( p ). Then 5 j= 1 j ( 5, ) T = X Binomial p. We test 1 H : p vs. 1 1 H : p>. We saw that a LRT is 1 if t k* φ ( x ) =. otherwise (t is the observed
More informationExtended Spectral Nonlinear Conjugate Gradient methods for solving unconstrained problems
International Journal of All Researh Euation an Sientifi Methos IJARESM ISSN: 55-6 Volume Issue 5 May-0 Extene Spetral Nonlinear Conjuate Graient methos for solvin unonstraine problems Dr Basim A Hassan
More informationProbabilistic analysis of the face stability of circular tunnels
Probabilisti analysis of the fae stability of irular tunnels Guilhem Mollon 1, Daniel Dias 2 and Abdul-Hamid Soubra 3, M. ASCE 1 PhD Student, INSA Lyon, LGCIE Site Coulomb 3, Géotehnique, Bât. J.C.A. Coulomb,
More informationThe First Integral Method for Solving a System of Nonlinear Partial Differential Equations
ISSN 179-889 (print), 179-897 (online) International Journal of Nonlinear Siene Vol.5(008) No.,pp.111-119 The First Integral Method for Solving a System of Nonlinear Partial Differential Equations Ahmed
More informationON THE MOVING BOUNDARY HITTING PROBABILITY FOR THE BROWNIAN MOTION. Dobromir P. Kralchev
Pliska Stud. Math. Bulgar. 8 2007, 83 94 STUDIA MATHEMATICA BULGARICA ON THE MOVING BOUNDARY HITTING PROBABILITY FOR THE BROWNIAN MOTION Dobromir P. Kralhev Consider the probability that the Brownian motion
More informationBrazilian Journal of Physics, vol. 29, no. 1, March,
Brazilian Journal of hysis, vol. 29, no., Marh, 999 79 Computational Methos Inspire by Tsallis Statistis: Monte Carlo an Moleular Dynamis Algorithms for the Simulation of Classial an Quantum Systems John
More informationMost results in this section are stated without proof.
Leture 8 Level 4 v2 he Expliit formula. Most results in this setion are stated without proof. Reall that we have shown that ζ (s has only one pole, a simple one at s =. It has trivial zeros at the negative
More informationAdvances in Radio Science
Advanes in adio Siene 2003) 1: 99 104 Copernius GmbH 2003 Advanes in adio Siene A hybrid method ombining the FDTD and a time domain boundary-integral equation marhing-on-in-time algorithm A Beker and V
More informationLecture 3 - Lorentz Transformations
Leture - Lorentz Transformations A Puzzle... Example A ruler is positioned perpendiular to a wall. A stik of length L flies by at speed v. It travels in front of the ruler, so that it obsures part of the
More information( ) ( ) Volumetric Properties of Pure Fluids, part 4. The generic cubic equation of state:
CE304, Spring 2004 Leture 6 Volumetri roperties of ure Fluids, part 4 The generi ubi equation of state: There are many possible equations of state (and many have been proposed) that have the same general
More informationEFFECT OF PITCH NUMBER IN OVERALL HEAT TRANSFER RATE IN DOUBLE PIPE HELICAL HEAT EXCHANGER
Volume 116 No. 5 2017, 1-6 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu EFFECT OF PITCH NUMBER IN OVERALL HEAT TRANSFER RATE IN DOUBLE PIPE HELICAL
More informationCMSC 451: Lecture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017
CMSC 451: Leture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017 Reading: Chapt 11 of KT and Set 54 of DPV Set Cover: An important lass of optimization problems involves overing a ertain domain,
More information11.4 Molecular Orbital Description of the Hydrogen Molecule Electron Configurations of Homonuclear Diatomic Molecules
Chap Moleular Eletroni Struture Table of Contents. The orn-oppenheimer pproximation -. The Hydrogen Moleule Ion.3 Calulation of the Energy of the Hydrogen Moleule Ion.4 Moleular Orbital Desription of the
More informationProceedings of Meetings on Acoustics
Proeedings of Meetings on Aoustis Volume 19, 2013 http://aoustialsoiety.org/ ICA 2013 Montreal Montreal, Canada 2-7 June 2013 Arhitetural Aoustis Session 4aAAa: Room Aoustis Computer Simulation I 4aAAa8.
More informationDeveloping Excel Macros for Solving Heat Diffusion Problems
Session 50 Developing Exel Maros for Solving Heat Diffusion Problems N. N. Sarker and M. A. Ketkar Department of Engineering Tehnology Prairie View A&M University Prairie View, TX 77446 Abstrat This paper
More informationGEOMETRIC AND STOCHASTIC ERROR MINIMISATION IN MOTION TRACKING. Karteek Alahari, Sujit Kuthirummal, C. V. Jawahar, P. J. Narayanan
GEOMETRIC AND STOCHASTIC ERROR MINIMISATION IN MOTION TRACKING Karteek Alahari, Sujit Kuthirummal, C. V. Jawahar, P. J. Narayanan Centre for Visual Information Tehnology International Institute of Information
More informationA NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS
Vietnam Journal of Mehanis, VAST, Vol. 4, No. (), pp. A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS Le Thanh Tung Hanoi University of Siene and Tehnology, Vietnam Abstrat. Conventional ship autopilots are
More informationA MATLAB Method of Lines Template for Evolution Equations
A MATLAB Metho of Lines Template for Evolution Equations H.S. Lee a, C.J. Matthews a, R.D. Braok a, G.C. Saner b an F. Ganola a a Faulty of Environmental Sienes, Griffith University, Nathan, QLD, 4111
More informationHankel Optimal Model Order Reduction 1
Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both
More informationHeat exchangers: Heat exchanger types:
Heat exhangers: he proess of heat exhange between two fluids that are at different temperatures and separated by a solid wall ours in many engineering appliations. he devie used to implement this exhange
More informationReliability-Based Approach for the Determination of the Required Compressive Strength of Concrete in Mix Design
Reliability-Based Approah for the Determination of the Required Compressive Strength of Conrete in Mix Design Nader M Okasha To ite this version: Nader M Okasha. Reliability-Based Approah for the Determination
More informationSOA/CAS MAY 2003 COURSE 1 EXAM SOLUTIONS
SOA/CAS MAY 2003 COURSE 1 EXAM SOLUTIONS Prepared by S. Broverman e-mail 2brove@rogers.om website http://members.rogers.om/2brove 1. We identify the following events:. - wathed gymnastis, ) - wathed baseball,
More informationHOW TO FACTOR. Next you reason that if it factors, then the factorization will look something like,
HOW TO FACTOR ax bx I now want to talk a bit about how to fator ax bx where all the oeffiients a, b, and are integers. The method that most people are taught these days in high shool (assuming you go to
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 4
Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas
More informationEECS 120 Signals & Systems University of California, Berkeley: Fall 2005 Gastpar November 16, Solutions to Exam 2
EECS 0 Signals & Systems University of California, Berkeley: Fall 005 Gastpar November 6, 005 Solutions to Exam Last name First name SID You have hour and 45 minutes to omplete this exam. he exam is losed-book
More informationCSE 5311 Notes 18: NP-Completeness
SE 53 Notes 8: NP-ompleteness (Last upate 7//3 8:3 PM) ELEMENTRY ONEPTS Satisfiability: ( p q) ( p q ) ( p q) ( p q ) Is there an assignment? (Deision Problem) Similar to ebugging a logi iruit - Is there
More informationDetermination of the reaction order
5/7/07 A quote of the wee (or amel of the wee): Apply yourself. Get all the eduation you an, but then... do something. Don't just stand there, mae it happen. Lee Iaoa Physial Chemistry GTM/5 reation order
More informationECE Microwave Engineering
ECE 5317-6351 Mirowave Engineering Aapte from notes by Prof. Jeffery T. Williams Fall 18 Prof. Davi R. Jakson Dept. of ECE Notes 7 Waveguiing Strutures Part : Attenuation ε, µσ, 1 Attenuation on Waveguiing
More information10.2 The Occurrence of Critical Flow; Controls
10. The Ourrene of Critial Flow; Controls In addition to the type of problem in whih both q and E are initially presribed; there is a problem whih is of pratial interest: Given a value of q, what fators
More informationA simple expression for radial distribution functions of pure fluids and mixtures
A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.
More informationF = F x x + F y. y + F z
ECTION 6: etor Calulus MATH20411 You met vetors in the first year. etor alulus is essentially alulus on vetors. We will need to differentiate vetors and perform integrals involving vetors. In partiular,
More information1. A dependent variable is also known as a(n). a. explanatory variable b. control variable c. predictor variable d. response variable ANSWER:
1. A epenent variale is also known as a(n). a. explanatory variale. ontrol variale. preitor variale. response variale FEEDBACK: A epenent variale is known as a response variale. Definition of the Simple
More informationMulti-scale Godunov-type method for cell-centered discrete Lagrangian hydrodynamics.
Multi-sale Gounov-type metho for ell-entere isrete Lagrangian hyroynamis. Pierre-Henri Maire, Bonifae Nkonga To ite this version: Pierre-Henri Maire, Bonifae Nkonga. Multi-sale Gounov-type metho for ell-entere
More informationThe tanh - coth Method for Soliton and Exact Solutions of the Sawada - Kotera Equation
Volume 117 No. 13 2017, 19-27 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu The tanh - oth Method for Soliton and Exat Solutions of the Sawada -
More informationProblem set 6 for the course Theoretical Optics Sample Solutions
Karlsruher Institut für Tehnologie KIT) Institut für theoretishe Festkörperphysik SS01 Prof. Dr. G. Shön, Dr. R. Frank 15.06.01 http://www.tfp.kit.eu/stuium-lehre.php Tutorial: Group 1, Name: Group, Group
More informationArray Design for Superresolution Direction-Finding Algorithms
Array Design for Superresolution Diretion-Finding Algorithms Naushad Hussein Dowlut BEng, ACGI, AMIEE Athanassios Manikas PhD, DIC, AMIEE, MIEEE Department of Eletrial Eletroni Engineering Imperial College
More informationEE 321 Project Spring 2018
EE 21 Projet Spring 2018 This ourse projet is intended to be an individual effort projet. The student is required to omplete the work individually, without help from anyone else. (The student may, however,
More informationSimultaneous and Sequential Auctions of Oligopoly Licenses
Simultaneous an Sequential Autions of Oligopoly Lienses Georgios Katsenos Institut für Mikroökonomik, Leibniz Universität Hannover September 1, 2007 Abstrat This paper ompares two proeures for alloating
More informationFREQUENCY DOMAIN FEEDFORWARD COMPENSATION. F.J. Pérez Castelo and R. Ferreiro Garcia
FREQUENCY DOMAIN FEEDFORWARD COMPENSATION F.J. Pérez Castelo and R. Ferreiro Garia Dept. Ingeniería Industrial. Universidad de La Coruña javierp@ud.es, Phone: 98 7.Fax: -98-7 ferreiro@ud.es, Phone: 98
More informationRecursive integral time extrapolation methods for scalar waves Paul J. Fowler*, Xiang Du, and Robin P. Fletcher, WesternGeco
Reursive integral time extrapolation methods for salar waves Paul J. Fowler*, Xiang Du, and Robin P. Flether, WesternGeo Summary We derive and ompare a variety of algorithms for reursive time extrapolation
More informationChapter 3 Lecture 7. Drag polar 2. Topics. Chapter-3
hapter 3 eture 7 Drag polar Topis 3..3 Summary of lift oeffiient, drag oeffiient, pithing moment oeffiient, entre of pressure and aerodynami entre of an airfoil 3..4 Examples of pressure oeffiient distributions
More informationTime Domain Method of Moments
Time Domain Method of Moments Massahusetts Institute of Tehnology 6.635 leture notes 1 Introdution The Method of Moments (MoM) introdued in the previous leture is widely used for solving integral equations
More informationLecture Notes 4 MORE DYNAMICS OF NEWTONIAN COSMOLOGY
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physis Department Physis 8.286: The Early Universe Otober 1, 218 Prof. Alan Guth Leture Notes 4 MORE DYNAMICS OF NEWTONIAN COSMOLOGY THE AGE OF A FLAT UNIVERSE: We
More informationThe Hanging Chain. John McCuan. January 19, 2006
The Hanging Chain John MCuan January 19, 2006 1 Introdution We onsider a hain of length L attahed to two points (a, u a and (b, u b in the plane. It is assumed that the hain hangs in the plane under a
More information7.1 Roots of a Polynomial
7.1 Roots of a Polynomial A. Purpose Given the oeffiients a i of a polynomial of degree n = NDEG > 0, a 1 z n + a 2 z n 1 +... + a n z + a n+1 with a 1 0, this subroutine omputes the NDEG roots of the
More informationGreen s Function for Potential Field Extrapolation
Green s Funtion for Potential Field Extrapolation. Soe Preliinaries on the Potential Magneti Field By definition, a potential agneti field is one for whih the eletri urrent density vanishes. That is, J
More informationInfluence of transverse cracks on the onset of delamination: application to L-angle specimens. F. Laurin*, A. Mavel, P. Nuñez, E.
Influene of transverse raks on the onset of delamination: appliation to L-angle speimens F. Laurin*, A. Mavel, P. Nuñez, E. Auguste Composite strutures subjeted to 3D loading Wings Strutures under 3D loadings
More informationSubject: Introduction to Component Matching and Off-Design Operation % % ( (1) R T % (
16.50 Leture 0 Subjet: Introdution to Component Mathing and Off-Design Operation At this point it is well to reflet on whih of the many parameters we have introdued (like M, τ, τ t, ϑ t, f, et.) are free
More informationOn Certain Singular Integral Equations Arising in the Analysis of Wellbore Recharge in Anisotropic Formations
On Certain Singular Integral Equations Arising in the Analysis of Wellbore Reharge in Anisotropi Formations C. Atkinson a, E. Sarris b, E. Gravanis b, P. Papanastasiou a Department of Mathematis, Imperial
More informationCOMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION
4 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION Jiri Nozika*, Josef Adame*, Daniel Hanus** *Department of Fluid Dynamis and
More informationarxiv:math/ v4 [math.ca] 29 Jul 2006
arxiv:math/0109v4 [math.ca] 9 Jul 006 Contiguous relations of hypergeometri series Raimundas Vidūnas University of Amsterdam Abstrat The 15 Gauss ontiguous relations for F 1 hypergeometri series imply
More informationSensitivity analysis for linear optimization problem with fuzzy data in the objective function
Sensitivity analysis for linear optimization problem with fuzzy data in the objetive funtion Stephan Dempe, Tatiana Starostina May 5, 2004 Abstrat Linear programming problems with fuzzy oeffiients in the
More informationNon-Markovian study of the relativistic magnetic-dipole spontaneous emission process of hydrogen-like atoms
NSTTUTE OF PHYSCS PUBLSHNG JOURNAL OF PHYSCS B: ATOMC, MOLECULAR AND OPTCAL PHYSCS J. Phys. B: At. Mol. Opt. Phys. 39 ) 7 85 doi:.88/953-75/39/8/ Non-Markovian study of the relativisti magneti-dipole spontaneous
More informationSupplementary Materials
Supplementary Materials Neural population partitioning and a onurrent brain-mahine interfae for sequential motor funtion Maryam M. Shanehi, Rollin C. Hu, Marissa Powers, Gregory W. Wornell, Emery N. Brown
More informationLecture 7: Sampling/Projections for Least-squares Approximation, Cont. 7 Sampling/Projections for Least-squares Approximation, Cont.
Stat60/CS94: Randomized Algorithms for Matries and Data Leture 7-09/5/013 Leture 7: Sampling/Projetions for Least-squares Approximation, Cont. Leturer: Mihael Mahoney Sribe: Mihael Mahoney Warning: these
More informationAyan Kumar Bandyopadhyay
Charaterization of radiating apertures using Multiple Multipole Method And Modeling and Optimization of a Spiral Antenna for Ground Penetrating Radar Appliations Ayan Kumar Bandyopadhyay FET-IESK, Otto-von-Guerike-University,
More informationImplementing the Law of Sines to solve SAS triangles
Implementing the Law of Sines to solve SAS triangles June 8, 009 Konstantine Zelator Dept. of Math an Computer Siene Rhoe Islan College 600 Mount Pleasant Avenue Proviene, RI 0908 U.S.A. e-mail : kzelator@ri.eu
More informationGrasp Planning: How to Choose a Suitable Task Wrench Space
Grasp Planning: How to Choose a Suitable Task Wrenh Spae Ch. Borst, M. Fisher and G. Hirzinger German Aerospae Center - DLR Institute for Robotis and Mehatronis 8223 Wessling, Germany Email: [Christoph.Borst,
More informationReview Topic 4: Cubic polynomials
Review Topi : ui polynomials Short answer Fatorise Px ( ) = x + 5x + x- 9 into linear fators. The polynomial Px ( ) = x - ax + x- leaves a remainer of when it is ivie y ( x - ) an a remainer of - when
More informationA Global Convergent Spectral Conjugate Gradient Method
Australian Journal of Basi an Applie Sienes 77: 30-309 03 ISSN 99-878 A Global Converent Spetral Conjuate Graient Metho Abbas Y. Al-Bayati an Hawraz N. Al-Khayat Collee of Basi Euation elafer Mosul University
More informationDevelopment of an efficient finite element model for the dynamic analysis of the train-bridge interaction
Development of an effiient finite element model for the dynami analysis of the train-bridge interation S. Neves, A. Azevedo & R. Calçada Faulty of Engineering, University of Porto, Porto, Portugal ABSTRACT:
More information