2D1240 Numerical Methods II / André Jaun, NADA, KTH NADA
|
|
- Hector Parker
- 6 years ago
- Views:
Transcription
1 A NADA 2D24 Numericl Methods II / André Jun, NADA, KTH. NON-LINEAR LEAST SQUARE FITTING, OPTIMIZATION.. Remember: wht we sw during the lst lesson Itertive solution for sprse sstems (Jcobi, Guss Seidel, help Mtlb sprse) Non-liner sstems Fied point, Piccrd nd Newton itertions Non-liner interpoltion..2 Overview: wht ou will lern tod Non-liner fitting Optimiztion serching for n etremum with-/ out using derivtives Ordinr differentil equtions using Euler nd Runge-Kutt methods
2 >< >: + b sin!(:5 t ) ß :3 + b sin!(:8 t ) ß :3 + b sin!(: t ) ß :5 + b sin!(:2 t ) ß :9 + b sin!(:5 t ) ß :4 + b sin!(:8 t ) ß : + b sin!(2: t ) ß :5 + b sin!(2:4 t ) ß :3 >< >: + b sin!(:5 t ) :3 ß + b sin!(:8 t ) :3 ß + b sin!(: t ) :5 ß + b sin!(:2 t ) :9 ß + b sin!(:5 t ) :4 ß + b sin!(:8 t ) : ß + b sin!(2: t ) :5 ß + b sin!(2:4 t ) :3 ß.2 Non-liner lest squre fitting (NAM 6., H 6.6) F (t) = + b sin!(t t ) t Emple Fit the model function to eperimentl dt This ields n overdetermined non-liner sstem (4 unknowns, 8 equtions) 8 8 or which cn be solved b demnding tht residuls between the mesurements i nd the F (t model ) i be k Fk smll ß 2 i.e.k Fk, = f(c) ß 2. Guss-Newton solution obtined from c new = c old + ffic; J T Jffic = J T f(c old ) The coefficients c re obtined from the Newton method for non-liner sstems, with n increment ffic solving the overdetermined liner sstem Jffic ß f.inmtlb, the solution of the norml equtions cn gin be computed with dc=-f J.
3 Ji=ones(n,); Ji2=sin(u); Ji3=b*(t-t).*cos(u); Ji4=-w*b*cos(u); J=[Ji Ji2 Ji3 Ji4]; % Jcobin dc=-j f; c=c+dc; % Guss-Newton b w t f dc step Emple solution of the non-liner fit with Mtlb t=[ ]'; =[ ]'; subplot(,2,), h=stem(t,), title('dt') =.7; b=.7; w=pi; t=.2; c=[ b w t]'; %initil guess n=size(t,); iter=;dcnorm=.; while dcnorm>e-6 & iter< u=w*(t-t); f=+b*sin(u)-;.6 Dt Model F=+b sinω(t t) (t i ).8 F(t) dcnorm=norm(dc); iter=iter+; =c(); b=c(2); w=c(3); t=c(4); D=[iter c b w t norm(f) norm(dc) ] end tt=(:.5:3)'; Ft=+b*sin(w*(tt-t)); subplot(,2,2), plot(t,,'o',tt,ft) t i t A solution with n ccurc better thn 6 is obtined with 6 steps:
4 ( c ; c ;;b); cluclte the i =@ j =[ ]'; %dt =[ ]'; = c = c = b = Emple: fit n F (; ) = ellipsis c) 2 ( 2 + ( c) 2 b 2 = to seven dt points. Set-up the overdetermined non-liner sstem with 7 equtions for 4 = 2 c = 2( c) = 2 c b = 2( c) b 3 Non liner fit to n ellipsis 2 c=; c=8; =8; b=3; %initil guess p=[c c b]'; iter=; dp=; while norm(dp)> e-6 & iter< iter=iter+; d=-c; d=-c; 8 f=d.^2/^2+d.^2/b^2-; J=-2*[d/^2 d/b^2 d.^2/^3 d.^2/b^3]; 6 dp=-j f; p=p+dp; %Guss-Newton end 4 c=p(), c=p(2), =p(3), b=p(4) v=:2*pi/6:2*pi; 2 plot(,,'o',c+*cos(v),c+b*sin(v)), is equl, Note tht fit to circle cn be written s liner lest squre problem (lb.3).
5 f=inline('.3-ep(-(-).^2)-min(/4,m(./(-),))'); [,f,flg]=fminbnd(f,,3).3 Optimiztion in one dimension (NAM 7., H 6.4) In Mtlb use fminbnd for rel rguments nd m for integer rguments rgument. =:.:3; plot(,fevl(f,),,f,'o') =.27 f =.342 flg = A vriet of methods re combined to provide rpid convergence nd robust nswers; the require unimodl functions tht re strictl incresing in ll the direction from the minimum:»» = ( )2 min f() ; m :3 e 4 Methods tht use order (Jcobin) nd 2 even order derivtives (Hessin) converge fster... if defined! Check st nd to set the properties. optimset help ;.3 ep[ ( ) 2 ] min(/4,m[/( ),]) A discountinuous unimodl function Here is n emple of n optimiztion where the rgument is n integer nd the function is chep to evlute: n=; f=rndn(,n); [fmin,ind]=min(f) fmin = ind = 98438
6 fi = fi p 5 Golden section serch of 2 [; b] minimum Λ without computing derivtives Two evlutions re required to decide if f() the minimum is locted in the left or right intervl: 6 )>f( if ) 2 f( then Λ 2 [ ; b] else Λ 2 [; 2 ] end For efficienc, choose the proportions so tht = b fi = b 2 b b 2 i i 2 b 2 b b b in order to reuse the previoius evlutions. B smmetr ( 2 )=(b ) = fi 2 = fi ) b 2 b = fi b b fi ) fi 2 + fi = ) fi = ß :68 2 which reminds the golden + rtio 5=2 ß :68 from the ntiquit. Finll p fi = + (b )( fi) 2 = + (b )fi
7 f=ff(n); function -cos(n/3+9); return f=ep(sqrt(n)/) =; b=44; n2=89; n=+b-n2; F=Ff(n); F2=Ff(n2); disp([ n n2 b]) while b-> if F>F2, b=n2; n2=n; F2=F; n=+b-n2; F=Ff(n); else =n; n=n2; F=F2; n2=+b-n; F2=Ff(n2); end disp([ n n2 b]) end; Fm=[F F2] n n2 n Fiboncci serch cn be used when the rgument is n integer Write reltion between steps in golden serch L i+2 = L i + L i+ This ields recursion formul for the Fiboncci numbers i+2 L i+ L L i+2 L i+ L i L i+ i F 2 f; ; 2; 3; 5; 8; 3; 2; 34; 55; 89; 44; 233;:::g where consecutive integers pproimte the golden section 34=55 = :682 ß fi. This suggests more efficient lgorithm thn Mtlb s m to find the mimum of Ff(n) in the intervl i = :::2 strting from lrger Fiboncci number: Fm =
8 A = ψ 4 + 2sin A = (; ) f 2 (; ) f.4 Optimiztion in higher dimensions n (NAM 7.2, H 6.5) Efficient methods such s steepest descent, conjugte-grdients) do not require evluting the Jcobin: the re the subject of dvnced courses. Newton s method is however sufficient if the function cn differentited nlticll to clculte n etremum b solving using itertions rφ( ; 2 ;:::; n ) = new = old + d; Jd = f( old ) Emple: mimum of Φ(; ) = ( + sin )e (2 + 2) round (; ) = ( 4 ; 4 ). ψ! ( 2( + sin ))e ( ( + sin ))e (2 + 2 ) (cos @ To solve this sstem of the form f() =, clculte the Jcobin J ij i =@ j =! 2 cos J = sin + 2 cos
9 =/4; =/4; F=(+sin())*ep(-(^2+^2)); z=[ ]'; dznorm=; iter=; disp([iter dznorm F]) In mtlb the solution is then obtined with while dznorm>e-6 & iter < 3 s=+sin(); f=[-2**s; cos()-2**s]; J=[-4*-2*sin() -2**cos(); -2* -2*-3*sin()-2**cos()]; dz=-j f; z=z+dz; dznorm=norm(dz,inf); iter=iter+; =z(); =z(2); F=(+sin())*ep(-(^2+^2)); disp([iter dznorm F]) Φ(,) end.4.6 dz F iter
10 . ORDINARY DIFFERENTIAL EQUATIONS (ODEs).2 Euler nd Runge-Kutt methods (NAM 8.2, H 9.3) Problem: numericll solve the ordinr differentil eqution for (t) for t > t d dt = f(t; ); with (t ) = using pproimtions with smll steps in time t ;t ;:::such tht t i+ = t i + h. Euler s method directl follows from the forwrd difference (t i + h) (t i ) ß f(t; ); ) (t i+ ) ß i + hf(t i ; i ) Strting from the initil condition (t ; ), one step with Euler s method produces n pproimtion (t ; ) with locl error O(h 2 ); fter n steps to rech the finl time t n = t + nh, the solution (t n ; n ) hs globl error no(h 2 ) ο O(h). h
11 fend=feuler(h) function t=; tend=2; n=tend/h; T=t; Y=; =6; Emple solve (t) = sin(t) with () = 6 using Euler s method 6.4 Euler method with h=.2,.,.5 nd.25 i=:n for =+h*f; t=t+h; f=sin(t*); 6.3 t]; Y=[Y; ]; T=[T; end 6.2 plot(t,y,t,y,'o'); fend=; return (t) 6. h=.2; Y2=[]; for k=:3, Y2=[Y2 feuler(h/2^k)], end; 6 Y2 = t A Richrdson etrpoltion cn be used to cncel the leding O(h) (globl) error in nd (2) ß 5:873 + clculte the vlue = 5: :873 5: Q Peer Teching minutes to think, eplin to our neighbour (2 nd vote) Richrdson etrpoltion. Which vlue for Q should ou use here bove? ψ "! 2
12 = f(t i ; i ) k 2 = f(t i + h; i + hk ) k 3 = f(t i + h 2 ; i + h 2 k 2) k 4 = f(t i + h; i + hk 3 ) k Runge-Kutt (RK2) chieves better precision with globl error in O(h 2 ) ( i+ = i + h 2 (k + k 2 ); with Runge-Kutt (RK4) chieves high precision with globl error in O(h 4 ) 8 >< = f(t i ; i ) k 2 = f(t i + h 2 ; i + h 2 k ) k i+ = i + h 6 (k + 2k 2 + 2k 3 + k 4 ); with In Mtlb use help ode23, ode45, odeset for RK methods with vrible step size,where the step size h is continull djusted to chieve specified precision with minimum number of steps. >:
13 end plot(t,y); h=.2; YN=[]; for k=:, YN=[YN rk4(h/2^k)]; Y2=YN(end), end; % --- Alterntive using Mtlb's solvers tend=2; =6; options=odeset('reltol',e-6,'abstol',e-4]); [T,Y]=ode23(@fsin,[ tend],,options); [T,Y]=ode45(@fsin,[ tend],,options); Emple solve (t) = sin(t) with () = 6 using the RK4 method f=fsin(t,), function f=sin(t*); return 6.4 Runge Kutt RK4 method with fied h=.2,. fend=rk4(h) function t=; tend=2; n=tend/h; =6; 6.3 Y=; h2=h/2; T=t; i=:n for 6.2 k=fsin(t, ); k2=fsin(t+h2,+h2*k); k3=fsin(t+h2,+h2*k2); (t) 6.,+h *k3); k4=fsin(t+h =+h/6*(k+2*k2+2*k3+k4); 6 t=t+h; t]; Y=[Y; ]; T=[T; 5.9 fend=; return t Y2 = A Richrdson etrpoltion gin cncels the leding (globl) O(h error, ) here in 4 5: with (2) 4 = 5: :8639 5:86346 ß 2
Calculus 2: Integration. Differentiation. Integration
Clculus 2: Integrtion The reverse process to differentition is known s integrtion. Differentition f() f () Integrtion As it is the opposite of finding the derivtive, the function obtined b integrtion is
More informationThe Islamic University of Gaza Faculty of Engineering Civil Engineering Department. Numerical Analysis ECIV Chapter 11
The Islmic University of Gz Fculty of Engineering Civil Engineering Deprtment Numericl Anlysis ECIV 6 Chpter Specil Mtrices nd Guss-Siedel Associte Prof Mzen Abultyef Civil Engineering Deprtment, The Islmic
More informationB.Sc. in Mathematics (Ordinary)
R48/0 DUBLIN INSTITUTE OF TECHNOLOGY KEVIN STREET, DUBLIN 8 B.Sc. in Mthemtics (Ordinry) SUPPLEMENTAL EXAMINATIONS 01 Numericl Methods Dr. D. Mckey Dr. C. Hills Dr. E.A. Cox Full mrks for complete nswers
More information1 Part II: Numerical Integration
Mth 4 Lb 1 Prt II: Numericl Integrtion This section includes severl techniques for getting pproimte numericl vlues for definite integrls without using ntiderivtives. Mthemticll, ect nswers re preferble
More informationExam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1
Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution
More informationThe Trapezoidal Rule
_.qd // : PM Pge 9 SECTION. Numericl Integrtion 9 f Section. The re of the region cn e pproimted using four trpezoids. Figure. = f( ) f( ) n The re of the first trpezoid is f f n. Figure. = Numericl Integrtion
More informationCalculus II: Integrations and Series
Clculus II: Integrtions nd Series August 7, 200 Integrls Suppose we hve generl function y = f(x) For simplicity, let f(x) > 0 nd f(x) continuous Denote F (x) = re under the grph of f in the intervl [,x]
More informationMath 1431 Section M TH 4:00 PM 6:00 PM Susan Wheeler Office Hours: Wed 6:00 7:00 PM Online ***NOTE LABS ARE MON AND WED
Mth 43 Section 4839 M TH 4: PM 6: PM Susn Wheeler swheeler@mth.uh.edu Office Hours: Wed 6: 7: PM Online ***NOTE LABS ARE MON AND WED t :3 PM to 3: pm ONLINE Approimting the re under curve given the type
More informationChapter 3 Solving Nonlinear Equations
Chpter 3 Solving Nonliner Equtions 3.1 Introduction The nonliner function of unknown vrible x is in the form of where n could be non-integer. Root is the numericl vlue of x tht stisfies f ( x) 0. Grphiclly,
More informationChapter 6 Notes, Larson/Hostetler 3e
Contents 6. Antiderivtives nd the Rules of Integrtion.......................... 6. Are nd the Definite Integrl.................................. 6.. Are............................................ 6. Reimnn
More informationOrdinary Differential Equations- Boundary Value Problem
Ordinry Differentil Equtions- Boundry Vlue Problem Shooting method Runge Kutt method Computer-bsed solutions o BVPFD subroutine (Fortrn IMSL subroutine tht Solves (prmeterized) system of differentil equtions
More informationy = f(x) This means that there must be a point, c, where the Figure 1
Clculus Investigtion A Men Slope TEACHER S Prt 1: Understnding the Men Vlue Theorem The Men Vlue Theorem for differentition sttes tht if f() is defined nd continuous over the intervl [, ], nd differentile
More informationTopic 1 Notes Jeremy Orloff
Topic 1 Notes Jerem Orloff 1 Introduction to differentil equtions 1.1 Gols 1. Know the definition of differentil eqution. 2. Know our first nd second most importnt equtions nd their solutions. 3. Be ble
More informationM344 - ADVANCED ENGINEERING MATHEMATICS
M3 - ADVANCED ENGINEERING MATHEMATICS Lecture 18: Lplce s Eqution, Anltic nd Numericl Solution Our emple of n elliptic prtil differentil eqution is Lplce s eqution, lso clled the Diffusion Eqution. If
More informationPolynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230
Polynomil Approimtions for the Nturl Logrithm nd Arctngent Functions Mth 23 You recll from first semester clculus how one cn use the derivtive to find n eqution for the tngent line to function t given
More informationFundamental Theorem of Calculus
Fundmentl Theorem of Clculus Recll tht if f is nonnegtive nd continuous on [, ], then the re under its grph etween nd is the definite integrl A= f() d Now, for in the intervl [, ], let A() e the re under
More informationZ b. f(x)dx. Yet in the above two cases we know what f(x) is. Sometimes, engineers want to calculate an area by computing I, but...
Chpter 7 Numericl Methods 7. Introduction In mny cses the integrl f(x)dx cn be found by finding function F (x) such tht F 0 (x) =f(x), nd using f(x)dx = F (b) F () which is known s the nlyticl (exct) solution.
More informationP 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)
1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this
More information1 Linear Least Squares
Lest Squres Pge 1 1 Liner Lest Squres I will try to be consistent in nottion, with n being the number of dt points, nd m < n being the number of prmeters in model function. We re interested in solving
More informationCalculus - Activity 1 Rate of change of a function at a point.
Nme: Clss: p 77 Mths Helper Plus Resource Set. Copright 00 Bruce A. Vughn, Techers Choice Softwre Clculus - Activit Rte of chnge of function t point. ) Strt Mths Helper Plus, then lod the file: Clculus
More information13.3. The Area Bounded by a Curve. Introduction. Prerequisites. Learning Outcomes
The Are Bounded b Curve 3.3 Introduction One of the importnt pplictions of integrtion is to find the re bounded b curve. Often such n re cn hve phsicl significnce like the work done b motor, or the distnce
More informationOptimization Lecture 1 Review of Differential Calculus for Functions of Single Variable.
Optimiztion Lecture 1 Review of Differentil Clculus for Functions of Single Vrible http://users.encs.concordi.c/~luisrod, Jnury 14 Outline Optimiztion Problems Rel Numbers nd Rel Vectors Open, Closed nd
More informationInterpreting Integrals and the Fundamental Theorem
Interpreting Integrls nd the Fundmentl Theorem Tody, we go further in interpreting the mening of the definite integrl. Using Units to Aid Interprettion We lredy know tht if f(t) is the rte of chnge of
More information13.3. The Area Bounded by a Curve. Introduction. Prerequisites. Learning Outcomes
The Are Bounded b Curve 3.3 Introduction One of the importnt pplictions of integrtion is to find the re bounded b curve. Often such n re cn hve phsicl significnce like the work done b motor, or the distnce
More informationA LEVEL TOPIC REVIEW. factor and remainder theorems
A LEVEL TOPIC REVIEW unit C fctor nd reminder theorems. Use the Fctor Theorem to show tht: ) ( ) is fctor of +. ( mrks) ( + ) is fctor of ( ) is fctor of + 7+. ( mrks) +. ( mrks). Use lgebric division
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More informationCS514 Fall '00 Numerical Analysis (Sketched) Solution of Homework 3 1. Questions from text in Chapter 3 Problem 1: 0 f 0 1 f 8 1 8(f 1 f 0 ) 1 f 4 2 8
CS514 Fll '00 Numericl Anlysis (Sketched) Solution of Homework 3 1. Questions from text in Chpter 3 Problem 1: 0 f 0 1 f 8 1 8(f 1 f 0 ) 1 f 4 2 8(f 2 f 1 ) 32(f 2 2f 1 f 0 ) 1 32 f 2 3 4(f 3 f 2 ) (f
More informationNow, given the derivative, can we find the function back? Can we antidifferenitate it?
Fundmentl Theorem of Clculus. Prt I Connection between integrtion nd differentition. Tody we will discuss reltionship between two mjor concepts of Clculus: integrtion nd differentition. We will show tht
More informationMathematics Extension 1
04 Bored of Studies Tril Emintions Mthemtics Etension Written by Crrotsticks & Trebl. Generl Instructions Totl Mrks 70 Reding time 5 minutes. Working time hours. Write using blck or blue pen. Blck pen
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 3 Polynomial Interpolation: Numerical Differentiation and Integration.
Advnced Computtionl Fluid Dynmics AA215A Lecture 3 Polynomil Interpoltion: Numericl Differentition nd Integrtion Antony Jmeson Winter Qurter, 2016, Stnford, CA Lst revised on Jnury 7, 2016 Contents 3 Polynomil
More informationPractice Final. Name: Problem 1. Show all of your work, label your answers clearly, and do not use a calculator.
Nme: MATH 2250 Clculus Eric Perkerson Dte: December 11, 2015 Prctice Finl Show ll of your work, lbel your nswers clerly, nd do not use clcultor. Problem 1 Compute the following limits, showing pproprite
More informationChapter Direct Method of Interpolation More Examples Electrical Engineering
Chpter. Direct Method of Interpoltion More Emples Electricl Engineering Emple hermistors re used to mesure the temperture of bodies. hermistors re bsed on mterils chnge in resistnce with temperture. o
More information4.6 Numerical Integration
.6 Numericl Integrtion 5.6 Numericl Integrtion Approimte definite integrl using the Trpezoidl Rule. Approimte definite integrl using Simpson s Rule. Anlze the pproimte errors in the Trpezoidl Rule nd Simpson
More informationNumerical Methods I Orthogonal Polynomials
Numericl Methods I Orthogonl Polynomils Aleksndr Donev Cournt Institute, NYU 1 donev@cournt.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fll 2014 Nov 6th, 2014 A. Donev (Cournt Institute) Lecture IX
More information1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation
1 1.1. Liner Constnt Coefficient Equtions Section Objective(s): Overview of Differentil Equtions. Liner Differentil Equtions. Solving Liner Differentil Equtions. The Initil Vlue Problem. 1.1.1. Overview
More informationNumerical Integration. Newton Cotes Formulas. Quadrature. Newton Cotes Formulas. To approximate the integral b
Numericl Integrtion Newton Cotes Formuls Given function f : R R nd two rel numbers, b R, < b, we clculte (pproximtely) the integrl I(f,, b) = f (x) dx K. Frischmuth (IfM UR) Numerics for CSE 08/09 8 /
More informationMath 270A: Numerical Linear Algebra
Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner
More informationf(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral
Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one
More informationThe problems that follow illustrate the methods covered in class. They are typical of the types of problems that will be on the tests.
ADVANCED CALCULUS PRACTICE PROBLEMS JAMES KEESLING The problems tht follow illustrte the methods covered in clss. They re typicl of the types of problems tht will be on the tests. 1. Riemnn Integrtion
More informationSOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 5 NOVEMBER 2014
SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 5 NOVEMBER 014 Mrk Scheme: Ech prt of Question 1 is worth four mrks which re wrded solely for the correct nswer.
More informationCalculus of Variations
Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function
More informationGoals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite
Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite
More informationLecture 20: Numerical Integration III
cs4: introduction to numericl nlysis /8/0 Lecture 0: Numericl Integrtion III Instructor: Professor Amos Ron Scribes: Mrk Cowlishw, Yunpeng Li, Nthnel Fillmore For the lst few lectures we hve discussed
More informationMATH , Calculus 2, Fall 2018
MATH 36-2, 36-3 Clculus 2, Fll 28 The FUNdmentl Theorem of Clculus Sections 5.4 nd 5.5 This worksheet focuses on the most importnt theorem in clculus. In fct, the Fundmentl Theorem of Clculus (FTC is rgubly
More informationMath 113 Exam 2 Practice
Mth Em Prctice Februry, 8 Em will cover sections 6.5, 7.-7.5 nd 7.8. This sheet hs three sections. The first section will remind you bout techniques nd formuls tht you should know. The second gives number
More informationFinal Exam - Review MATH Spring 2017
Finl Exm - Review MATH 5 - Spring 7 Chpter, 3, nd Sections 5.-5.5, 5.7 Finl Exm: Tuesdy 5/9, :3-7:pm The following is list of importnt concepts from the sections which were not covered by Midterm Exm or.
More informationLog1 Contest Round 3 Theta Individual. 4 points each 1 What is the sum of the first 5 Fibonacci numbers if the first two are 1, 1?
008 009 Log1 Contest Round Thet Individul Nme: points ech 1 Wht is the sum of the first Fiboncci numbers if the first two re 1, 1? If two crds re drwn from stndrd crd deck, wht is the probbility of drwing
More informationARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac
REVIEW OF ALGEBRA Here we review the bsic rules nd procedures of lgebr tht you need to know in order to be successful in clculus. ARITHMETIC OPERATIONS The rel numbers hve the following properties: b b
More informationMAT137 Calculus! Lecture 20
officil website http://uoft.me/mat137 MAT137 Clculus! Lecture 20 Tody: 4.6 Concvity 4.7 Asypmtotes Net: 4.8 Curve Sketching 4.5 More Optimiztion Problems MVT Applictions Emple 1 Let f () = 3 27 20. 1 Find
More informationx ) dx dx x sec x over the interval (, ).
Curve on 6 For -, () Evlute the integrl, n (b) check your nswer by ifferentiting. ( ). ( ). ( ).. 6. sin cos 7. sec csccot 8. sec (sec tn ) 9. sin csc. Evlute the integrl sin by multiplying the numertor
More informationThe First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).
The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples
More informationChapter 3 Polynomials
Dr M DRAIEF As described in the introduction of Chpter 1, pplictions of solving liner equtions rise in number of different settings In prticulr, we will in this chpter focus on the problem of modelling
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More informationTotal Score Maximum
Lst Nme: Mth 8: Honours Clculus II Dr. J. Bowmn 9: : April 5, 7 Finl Em First Nme: Student ID: Question 4 5 6 7 Totl Score Mimum 6 4 8 9 4 No clcultors or formul sheets. Check tht you hve 6 pges.. Find
More information6.2 CONCEPTS FOR ADVANCED MATHEMATICS, C2 (4752) AS
6. CONCEPTS FOR ADVANCED MATHEMATICS, C (475) AS Objectives To introduce students to number of topics which re fundmentl to the dvnced study of mthemtics. Assessment Emintion (7 mrks) 1 hour 30 minutes.
More informationOperations with Polynomials
38 Chpter P Prerequisites P.4 Opertions with Polynomils Wht you should lern: How to identify the leding coefficients nd degrees of polynomils How to dd nd subtrct polynomils How to multiply polynomils
More informationAPPM 1360 Exam 2 Spring 2016
APPM 6 Em Spring 6. 8 pts, 7 pts ech For ech of the following prts, let f + nd g 4. For prts, b, nd c, set up, but do not evlute, the integrl needed to find the requested informtion. The volume of the
More informationA-Level Mathematics Transition Task (compulsory for all maths students and all further maths student)
A-Level Mthemtics Trnsition Tsk (compulsory for ll mths students nd ll further mths student) Due: st Lesson of the yer. Length: - hours work (depending on prior knowledge) This trnsition tsk provides revision
More informationDERIVATIVES NOTES HARRIS MATH CAMP Introduction
f DERIVATIVES NOTES HARRIS MATH CAMP 208. Introduction Reding: Section 2. The derivtive of function t point is the slope of the tngent line to the function t tht point. Wht does this men, nd how do we
More informationHIGHER SCHOOL CERTIFICATE EXAMINATION MATHEMATICS 3 UNIT (ADDITIONAL) AND 3/4 UNIT (COMMON) Time allowed Two hours (Plus 5 minutes reading time)
HIGHER SCHOOL CERTIFICATE EXAMINATION 998 MATHEMATICS 3 UNIT (ADDITIONAL) AND 3/4 UNIT (COMMON) Time llowed Two hours (Plus 5 minutes reding time) DIRECTIONS TO CANDIDATES Attempt ALL questions ALL questions
More informationMatrices and Determinants
Nme Chpter 8 Mtrices nd Determinnts Section 8.1 Mtrices nd Systems of Equtions Objective: In this lesson you lerned how to use mtrices, Gussin elimintion, nd Guss-Jordn elimintion to solve systems of liner
More informationSOLVING SYSTEMS OF EQUATIONS, ITERATIVE METHODS
ELM Numericl Anlysis Dr Muhrrem Mercimek SOLVING SYSTEMS OF EQUATIONS, ITERATIVE METHODS ELM Numericl Anlysis Some of the contents re dopted from Lurene V. Fusett, Applied Numericl Anlysis using MATLAB.
More informationOverview of Calculus I
Overview of Clculus I Prof. Jim Swift Northern Arizon University There re three key concepts in clculus: The limit, the derivtive, nd the integrl. You need to understnd the definitions of these three things,
More informationSpace Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space.
Clculus 3 Li Vs Spce Curves Recll the prmetric equtions of curve in xy-plne nd compre them with prmetric equtions of curve in spce. Prmetric curve in plne x = x(t) y = y(t) Prmetric curve in spce x = x(t)
More informationBIFURCATIONS IN ONE-DIMENSIONAL DISCRETE SYSTEMS
BIFRCATIONS IN ONE-DIMENSIONAL DISCRETE SYSTEMS FRANCESCA AICARDI In this lesson we will study the simplest dynmicl systems. We will see, however, tht even in this cse the scenrio of different possible
More informationBest Approximation. Chapter The General Case
Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given
More information3.4 Numerical integration
3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,
More informationImproper Integrals. Type I Improper Integrals How do we evaluate an integral such as
Improper Integrls Two different types of integrls cn qulify s improper. The first type of improper integrl (which we will refer to s Type I) involves evluting n integrl over n infinite region. In the grph
More informationDefinite integral. Mathematics FRDIS MENDELU
Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová Brno 1 Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function defined on [, b]. Wht is the re of the
More informationD 4u 2 p. u D 12 C p 144 C 4.72/ 16. D 12 C 12p r 2 D 5 C 3p 3. r D c D
pu D 4u 2 p 4u u C 2 C 36u C 54 SECTION 4.8 Newton s Method 491 Therefore the eqution becomes 8.S.0/=/ 2 D 8.S. 2 /=/ 2 8u 2 C 40u C 48 D 4u 2 4u. p u p u C 2/ C 36u C 54 4u 2 C 4u 6 D 4u. p u p u C 2/
More informationSUMMER KNOWHOW STUDY AND LEARNING CENTRE
SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18
More informationu t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx
Chpter 9: Green s functions for time-independent problems Introductory emples One-dimensionl het eqution Consider the one-dimensionl het eqution with boundry conditions nd initil condition We lredy know
More informationMCR 3U Exam Review. 1. Determine which of the following equations represent functions. Explain. Include a graph. 2. y x
MCR U MCR U Em Review Introduction to Functions. Determine which of the following equtions represent functions. Eplin. Include grph. ) b) c) d) 0. Stte the domin nd rnge for ech reltion in question.. If
More informationn f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1
The Fundmentl Theorem of Clculus As we continue to study the re problem, let s think bck to wht we know bout computing res of regions enclosed by curves. If we wnt to find the re of the region below the
More informationPartial Derivatives. Limits. For a single variable function f (x), the limit lim
Limits Prtil Derivtives For single vrible function f (x), the limit lim x f (x) exists only if the right-hnd side limit equls to the left-hnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles
More informationUNIVERSITY OF BOLTON SCHOOL OF ENGINEERING B.ENG (HONS) ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATION SEMESTER /2018
ENG005 B.ENG (HONS) ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATION SEMESTER 1-017/018 MODULE NO: EEE4001 Dte: 19Jnury 018 Time:.00 4.00 INSTRUCTIONS TO CANDIDATES: There re SIX questions. Answer ANY
More informationLINEAR ALGEBRA APPLIED
5.5 Applictions of Inner Product Spces 5.5 Applictions of Inner Product Spces 7 Find the cross product of two vectors in R. Find the liner or qudrtic lest squres pproimtion of function. Find the nth-order
More informationDefinite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30
Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová (Mendel University) Definite integrl MENDELU / Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function
More informationYear 12 Trial Examination Mathematics Extension 1. Question One 12 marks (Start on a new page) Marks
THGS Mthemtics etension Tril 00 Yer Tril Emintion Mthemtics Etension Question One mrks (Strt on new pge) Mrks ) If P is the point (-, 5) nd Q is the point (, -), find the co-ordintes of the point R which
More informationFirst Semester Review Calculus BC
First Semester Review lculus. Wht is the coordinte of the point of inflection on the grph of Multiple hoice: No lcultor y 3 3 5 4? 5 0 0 3 5 0. The grph of piecewise-liner function f, for 4, is shown below.
More informationES.182A Topic 32 Notes Jeremy Orloff
ES.8A Topic 3 Notes Jerem Orloff 3 Polr coordintes nd double integrls 3. Polr Coordintes (, ) = (r cos(θ), r sin(θ)) r θ Stndrd,, r, θ tringle Polr coordintes re just stndrd trigonometric reltions. In
More informationES.182A Topic 30 Notes Jeremy Orloff
ES82A opic 3 Notes Jerem Orloff 3 Non-independent vribles: chin rule Recll the chin rule: If w = f, ; nd = r, t, = r, t then = + r t r t r t = + t t t r nfortuntel, sometimes there re more complicted reltions
More informationTopics Covered AP Calculus AB
Topics Covered AP Clculus AB ) Elementry Functions ) Properties of Functions i) A function f is defined s set of ll ordered pirs (, y), such tht for ech element, there corresponds ectly one element y.
More informationOrthogonal Polynomials and Least-Squares Approximations to Functions
Chpter Orthogonl Polynomils nd Lest-Squres Approximtions to Functions **4/5/3 ET. Discrete Lest-Squres Approximtions Given set of dt points (x,y ), (x,y ),..., (x m,y m ), norml nd useful prctice in mny
More informationInstructor: Marios M. Fyrillas HOMEWORK ASSIGNMENT ON INTERPOLATION
AMAT 34 Numericl Methods Instructor: Mrios M. Fyrills Emil: m.fyrills@fit.c.cy Office Tel.: 34559/6 Et. 3 HOMEWORK ASSIGNMENT ON INTERPOATION QUESTION Using interpoltion by colloction-polynomil fit method
More informationCAAM 453 NUMERICAL ANALYSIS I Examination There are four questions, plus a bonus. Do not look at them until you begin the exam.
Exmintion 1 Posted 23 October 2002. Due no lter thn 5pm on Mondy, 28 October 2002. Instructions: 1. Time limit: 3 uninterrupted hours. 2. There re four questions, plus bonus. Do not look t them until you
More informationSection 4: Integration ECO4112F 2011
Reding: Ching Chpter Section : Integrtion ECOF Note: These notes do not fully cover the mteril in Ching, ut re ment to supplement your reding in Ching. Thus fr the optimistion you hve covered hs een sttic
More informationNumerical Integration
Chpter 1 Numericl Integrtion Numericl differentition methods compute pproximtions to the derivtive of function from known vlues of the function. Numericl integrtion uses the sme informtion to compute numericl
More informationa n = 1 58 a n+1 1 = 57a n + 1 a n = 56(a n 1) 57 so 0 a n+1 1, and the required result is true, by induction.
MAS221(216-17) Exm Solutions 1. (i) A is () bounded bove if there exists K R so tht K for ll A ; (b) it is bounded below if there exists L R so tht L for ll A. e.g. the set { n; n N} is bounded bove (by
More information37 Kragujevac J. Math. 23 (2001) A NOTE ON DENSITY OF THE ZEROS OF ff-orthogonal POLYNOMIALS Gradimir V. Milovanović a and Miodrag M. Spalević
37 Krgujevc J. Mth. 23 (2001) 37 43. A NOTE ON DENSITY OF THE ZEROS OF ff-orthogonal POLYNOMIALS Grdimir V. Milovnović nd Miodrg M. Splević b Fculty of Electronic Engineering, Deprtment of Mthemtics, University
More information5.5 The Substitution Rule
5.5 The Substitution Rule Given the usefulness of the Fundmentl Theorem, we wnt some helpful methods for finding ntiderivtives. At the moment, if n nti-derivtive is not esily recognizble, then we re in
More informationCalculus AB Section I Part A A CALCULATOR MAY NOT BE USED ON THIS PART OF THE EXAMINATION
lculus Section I Prt LULTOR MY NOT US ON THIS PRT OF TH XMINTION In this test: Unless otherwise specified, the domin of function f is ssumed to e the set of ll rel numers for which f () is rel numer..
More informationMA 124 January 18, Derivatives are. Integrals are.
MA 124 Jnury 18, 2018 Prof PB s one-minute introduction to clculus Derivtives re. Integrls re. In Clculus 1, we lern limits, derivtives, some pplictions of derivtives, indefinite integrls, definite integrls,
More informationBest Approximation in the 2-norm
Jim Lmbers MAT 77 Fll Semester 1-11 Lecture 1 Notes These notes correspond to Sections 9. nd 9.3 in the text. Best Approximtion in the -norm Suppose tht we wish to obtin function f n (x) tht is liner combintion
More informationMultiple Integrals. Review of Single Integrals. Planar Area. Volume of Solid of Revolution
Multiple Integrls eview of Single Integrls eding Trim 7.1 eview Appliction of Integrls: Are 7. eview Appliction of Integrls: Volumes 7.3 eview Appliction of Integrls: Lengths of Curves Assignment web pge
More informationAQA Further Pure 2. Hyperbolic Functions. Section 2: The inverse hyperbolic functions
Hperbolic Functions Section : The inverse hperbolic functions Notes nd Emples These notes contin subsections on The inverse hperbolic functions Integrtion using the inverse hperbolic functions Logrithmic
More informationHOMEWORK SOLUTIONS MATH 1910 Sections 7.9, 8.1 Fall 2016
HOMEWORK SOLUTIONS MATH 9 Sections 7.9, 8. Fll 6 Problem 7.9.33 Show tht for ny constnts M,, nd, the function yt) = )) t ) M + tnh stisfies the logistic eqution: y SOLUTION. Let Then nd Finlly, y = y M
More informationWe know that if f is a continuous nonnegative function on the interval [a, b], then b
1 Ares Between Curves c 22 Donld Kreider nd Dwight Lhr We know tht if f is continuous nonnegtive function on the intervl [, b], then f(x) dx is the re under the grph of f nd bove the intervl. We re going
More information1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ),
1. Guss-Jcobi qudrture nd Legendre polynomils Simpson s rule for evluting n integrl f(t)dt gives the correct nswer with error of bout O(n 4 ) (with constnt tht depends on f, in prticulr, it depends on
More information