Tricks of the Trade. Edited by Paul Abbott. Interpolation with Noise Vittorio G. Caffa. The Mathematica Journal
|
|
- Antony Montgomery
- 5 years ago
- Views:
Transcription
1 The Mathematica Journal Tricks of the Trade Edited by Paul Abbott This is a column of programming tricks and techniques, most of which, we hope, will be contributed by our readers, either directly as submissions to The Mathematica Journal or as an edited answer to a question posted in the Mathematica newsgroup, comp.soft-sys.math.mathematica. Interpolation with Noise Vittorio G. Caffa caffa@iabg.de In Interpolation, TMJ, 9(), 004, pp. 06 0, a method was presented for computing derivatives of a sampled function. This method only works if the data is unaffected by measuring error. If models of the data source and the measurement noise are available, you can use (for example) Kalman filtering techniques for computing an optimal estimation (see mathworld.wolfram.com/ KalmanFilter.html). Alternatively, you can achieve good results using a simple filtering technique. Consider again the exact function yhxl = x sinhxl - x ê. In[]:= yhx_l = x sinhxl - x ê ; Take x from a uniform random distribution over the interval Ha, bl = H, L, and append the endpoints. In[]:= 8a, b< = 8, <; SeedRandom@4D; xs = Join@8a<, Table@Random@Real, 8a, b<d, 8<D, 8b<D; Now sample yhxl over the interval adding Gaussian noise corresponding to measurement error. In[]:= Needs@"Statistics`ContinuousDistributions`"D; f Hx_L := 8x, yhxl + Random@NormalDistribution@0, 0.0DD<; data = f êû xs; Superimpose the sampled data over a plot of its interpolated function to visualie the 8x, y< data. It is clear that the interpolation is not smooth, so computing derivatives will be problematic. In[6]:= y int = Interpolation@dataD;
2 Tricks of the Trade 68 In[7]:= y int HxL, 8x, a, b<, Epilog Æ AbsolutePointSie@D, Point êû data<d A simple filtering technique involves computing a truncated Fourier series of the function yhxl, that is, n y n HxL = a 0 + Ha k coshk w xl + b k sinhk w xll. k= Requiring the wavelength l to satisfy l = p ê w > Hb - al avoids wrap-around effects. As only the frequencies up to n w are considered, the method acts like a low-pass filter. The value of the parameter n should be tuned to obtain the best results. This approach should be compared with that used in Spectral Analysis of Irregularly Sampled Data, TMJ, 7(), 997, pp Here is the Fourier basis set for arbitrary n and l. In[8]:= basis@n_, l_d := WithA9w = ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ p, k = Range@nD=, Join@8<, coshk w xl, sinhk w xlde N@lD The best truncated Fourier expansion, y n, for l =0 > Hb - al = and n = is computed using Fit. In[9]:= y = Function@x, Evaluate@Fit@data, basis@, 0D, xddd Out[9]= Function@x, cosh xl cosh xl cosh0.689 xl cosh xl -.8 cosh.047 xl sinh xl sinh xl -.0 sinh0.689 xl sinh xl +.6 sinh.047 xl D Here is the sampled data superimposed over a plot of y.
3 686 Edited by Paul Abbott In[0]:= y HxL, 8x, a, b<, Epilog Æ AbsolutePointSie@D, Point êû data<d Compare the first derivatives; the red curve denotes the exact value ( ), the green curve denotes the interpolated function ( ), and the blue curve denotes the truncated Fourier fit y ( ). The agreement between the exact function and the truncated Fourier fit is excellent. In[]:= PlotA8 y HxL, y int HxL, y HxL<, 8x, a, b<, PlotStyle Æ 9Hue@0D, HueA ÄÄÄÄÄ E, HueA ÄÄÄÄÄÄ E=E Using Reduce To Solve The Kuhn Tucker Equations Frank J. Kampas fkampas@msn.com The functions Maximie and Minimie give exact solutions to simple nonlinear optimiation problems with equality and inequality constraints. Another approach is to set up Lagrange multipliers for equality constraints and the Kuhn Tucker equations for inequality constraints (see mathworld.wolfram.com/ Kuhn-TuckerTheorem.html) and solve using Reduce. In some cases, this approach can solve problems which cannot directly be solved with Maximie and Minimie. The following function demonstrates this approach.
4 Tricks of the Trade 687 In[]:= KuhnTuckerHob_, cons_list, vars_, domain L := ModuleA8consconvrule = 8x_ y_ Æ y - x 0, x_ > y_ Æ y - x < 0, x_ y_ Æ x - y 0, x_ y_ Æ x - y 0<, stdcons, eqcons, ineqcons, lambdas, mus, numequals, numinequals, lagrangian, eqs, eqs, eqs<, stdcons = cons ê. consconvrule; eqcons = Cases@stdcons, x_ 0 xd; ineqcons = Cases@stdcons, x_ 0 xd; numequals = Length@eqconsD; numinequals = Length@ineqconsD; lambdas = Array@l, numequalsd; mus = Array@m, numinequalsd; lagrangian = ob + lambdas.eqcons + mus.ineqcons; eqs = i lagrangian ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ 0& y êû vars; eqs = Thread@mus 0D; k # { eqs = Table@musPiT ineqconspit 0, 8i, numinequals<d; Reduce@ Join@eqs, eqs, eqs, consd, Join@vars, lambdas, musd, domain, Backsubstitution Æ TrueDE Here is an example with an equality constraint in which both the maximum and minimum are obtained. Defining the obective function with a constraint and listing it first in the list of variables places the obective function first in the results. In[]:= Out[]= KuhnTucker@obf, 8obf ä x + y +, x + y + ä <, 8obf, x, y, <D i obf - Ì x - k Ì y - Ì - Ì y lhl - Ì lhl Î i obf Ì x { k Ì y y Ì Ì lhl - Ì lhl - Equivalent results are obtained by using Maximie and Minimie. In[]:= Out[]= In[4]:= Out[4]= Maximie@obf, 8obf ä x + y +, x + y + ä <, 8obf, x, y, <D êê Simplify 9,9obf Ø, x Ø, y Ø, Ø == Minimie@obf, 8obf ä x + y +, x + y + ä <, 8obf, x, y, <D êê Simplify 9-,9obf Ø-,x Ø-, y Ø-, Ø- == KuhnTucker can solve the problem when the righthand side of the constraint is changed from a number to a symbol. {
5 688 Edited by Paul Abbott In[]:= Out[]= 8obf ä x + y +, x + y + ä c<, 8obf, x, y, <D i obf - è!! è!! è!! c c c Ì x - k Ì y - Ì è!! c y - Ì lhl - Ì c 0 Ì lhl ÅÅÅÅ è!! Î c { i obf è!! è!! è!! è!! c c c c Ì x k Ì y Ì Ì y lhl - Ì c 0 Ì lhl - ÅÅÅÅ è!! c { Maximie and Minimie give an error message in this situation. In[6]:= Maximie@obf, 8obf ä x + y +, x + y + ä c<, 8obf, x, y, <D Maximie::consv : The constraint x + y + c contains a nonconstant expression c independent of variables 8obf, x, y, <. More Out[6]= Maximie@obf, 8obf x + y +, x + y + c<, 8obf, x, y, <D In the following example, the constraint x b is only active if b is less than. Therefore, the solution obtained covers b and b separately. In[7]:= KuhnTucker@obf, 8obf ä x - x, x b<, 8obf, x<d Out[7]= Hb Ï obf b - b Ï x b Ï lhl - Ï mhl - Hb - LL Í Hb Ï obf - Ï x Ï lhl - Ï mhl 0L When the solution is degenerate, all the results can be obtained. In[8]:= KuhnTucker@obf, 8obf ä x + y, x + y ä <, 8obf, x, y<d Out[8]= obf Ï y - x Ï lhl - Ï lhl - The following problem can be solved quite easily if the obective function is not defined as a constraint. Adding the extra constraint greatly increases the solution time. The problem cannot be solved by Minimie or NMinimie unless a change of variables is used to eliminate the square roots. In[9]:= KuhnTuckerA ÄÄÄÄÄÄ Jt h - "###### q h + t l - "##### q l N, 8t l 0 q l, t h q h, t l - 0 q l t h - 0 q h, t h - q h t l - q l <, 8q l, q h, t l, t h <E Out[9]= q l ÅÅÅ 776 Ì q h ÅÅÅÅÅÅÅÅ 6 Ì t l ÅÅÅ 444 Ì t h 8 ÅÅÅ 444 Ì mhl Ì mhl 0 Ì mhl 0 Ì mh4l ÅÅÅÅÅ Using Reduce with domain specified eliminates complex solution candidates.
6 Tricks of the Trade 689 In[0]:= 8obf == x + y, x + y ä <, 8obf, x, y<, RealsD êê FullSimplify Out[0]= obf ê Ì x Ì y Ì lhl - Ì lhl - ÅÅÅÅÅ Note, however, that using Lagrange multipliers and the Kuhn Tucker equations relies on assumptions that the code does not check. For instance, the solution set of equational constraints should be a manifold of dimension v - e, where v is the number of variables and e is the number of equations. If this is not satisfied, KuhnTucker may not find the extremum. In[]:= Out[]= KuhnTucker@x + y, 8 ä x - xy, ä y + xy, ä x + y <, 8x, y, <D False Maximie and Minimie work in this situation. In[]:= Minimie@x + y, 8 ä x - xy, ä y + xy, ä x + y <, 8x, y, <D Out[]= 80, 8x Ø 0, y Ø 0, Ø 0<< In[]:= Maximie@x + y, 8 ä x - xy, ä y + xy, ä x + y <, 8x, y, <D Out[]= 80, 8x Ø 0, y Ø 0, Ø 0<< Legendre Gauss Quadrature Gaussian quadrature gives the best estimate of an integral by picking optimal abscissas, x i, at which to evaluate the function f HxL. It is optimal because the n-point formula is exact for polynomials of degree n -. Legendre Gauss quadrature (see mathworld.wolfram.com/legendre-gaussquadrature.html ) is a Gaussian quadrature over the D with weighting function W HxL =, that is, - f HxL x º w i f Hx i L, where the abscissas, x i, for quadrature order n, satisfy P n Hx i L 0, and the weights are H - x i L w i ÅÅÅÅÅÅÅÅÅ Hn + L P n+ Hx i L. Define the abscissas as the roots of P n HxL using Root and save them as they are computed. This is an exact expression for the roots. In[]:= n i= x ê: x i_,n_ := x ê: x i,n = Root@Function@x, P n HxLD, id Define the weights. RootReduce expresses the weights as a single Root obect. Save the weights as they are computed.
7 690 Edited by Paul Abbott H - x i,n L In[]:= w ê: w i_,n_ := w ê: w i,n = ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ ÄÄÄÄÄÄÄÄÄÄÄÄ Hn + L P n+ Hx i,n L Here are the x i and w i for n =. In[]:= Out[]= With@8n = <, Table@8x i,n, w i,n <, 8i, n<dd i - "##### ÅÅÅÅ 0 "##### k ÅÅÅÅ ÅÅÅÅ 9 8 ÅÅÅÅ 9 ÅÅÅÅ 9 y { Here is x and w for n =, expressed as radicals. In[4]:= With@8n = <, 8x,n, w,n <êêtoradicalsd Out[4]= 9- ÅÅÅÅÅ $%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%!!! ÅÅÅÅÅ I - 70 M, 7!!! ÅÅÅÅÅÅÅÅÅÅÅ I + 70 M= 900 êê RootReduce Abscissas and weights for higher n are better left as Root obects (see Root versus Radicals, TMJ, 9(), 00, pp. 7). Note that, for fixed n, each w i can be expressed as a particular root of a polynomial of degree d ÅÅÅÅ n t. Here are the weights for n = 6. In[]:= With@8n = 6<, Table@w i,n, 8i, n<dd Out[]= 8Root@0000 # # + 69 # &, D, Root@0000 # # + 69 # &, D, Root@0000 # # + 69 # &, D, Root@0000 # # + 69 # &, D, Root@0000 # # + 69 # &, D, Root@0000 # # + 69 # &, D< n The weights satisfy i= In[6]:= Out[6]= Tr@%D êêrootreduce w i. Here we check this for n = 6. Now consider computing the integral of a polynomial of degree m, denoted p m HxL. In[7]:= p m_ Hx_L := a 0 + a x ; m = Compute the integral of p 7 HxL over - x. In[8]:= - p 7 HxL x Out[8]= a 0 + a + a 4 + a 6 7
8 Tricks of the Trade 69 As expected, an identical answer is obtained using Legendre Gauss quadrature with n = 4 since n - = 7. n In[9]:= WithA8n = 4<, w i,n p 7 Hx i,n LE êê FullSimplify Out[9]= a 0 + a i= + a 4 + a 6 7 Pascal Matrices Writing Pascal s triangle as a lower triangular matrix generates one type of Pascal matrix (see mathworld.wolfram.com/pascalmatrix.html). Such matrices have many interesting properties (web.mit.edu/8.06/www/pascal-work.pdf). Define the lower triangular matrix L n by 8L n < i, = i i y, where i, = 0,,, k { n -. In[]:= Here is L. L n_ := TableA i i y, 8i, 0,n - <, 8, 0,n - <E k { In[]:= L Out[]= i y k { And here is its inverse. In[]:= L - Out[]= i y k { In general, the inverse of L n has 8L - n < i, = H-L i+ i i y. k { Clearly 8L p < i, = p i- i i y. k {
9 69 Edited by Paul Abbott In[4]:= pd Out[4]= i y p p p 0 0 p p p 0 k p 4 4 p 6 p 4 p { Define v n HxL = 8x i < i=0,,, n. In[]:= v n_ Hx_L := x Range@nD- Using the binomial theorem, we can show that Ln p v n HxL = v n Hp + xl, that is, i y i y i y p x p + x p p x = Hp + xl. p p p x Hp + xl k ª ª ª ª { k ª { k ª { Prove the binomial theorem. In[6]:= SimplifyAHp + xl i i i y x p i-, i Œ E k { i =0 Out[6]= True Verify that Ln p v n HxL = v n Hp + xl for n =. In[7]:= MatrixPower@L, pd.v HxL êêfactor Out[7]= 8, p + x, Hp + xl, Hp + xl, Hp + xl 4 < In[8]:= Out[8]= % v Hp + xl True Writing Ln p = A n p, it follows that A n = loghl n L = lim pø0 HLn p - I n Lê p. Here is A. In[9]:= Out[9]= A = lim pæ0 i y k { MatrixPower@L, pd - IdentityMatrix@D ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ p In general, it can be shown that 8A n < i, = i d i, +. Check that L = A.
10 Tricks of the Trade 69 In[0]:= D Out[0]= i y k { Define the upper triangular matrix n by alternating the sign of the rows in L n T. In[]:= Here is. n_ := H-L Range@nD- L n T In[]:= Out[]= i y k { n is its own inverse (it is involutory). In[]:=. Out[]= i y k { Multiplying L n by its transpose yields a symmetric, positive-definite matrix, denoted S n, now with Pascal s triangle entries along each skew diagonal (see mathworld.wolfram.com/skewdiagonal.html). In[4]:= S n_ := L n.l n T In[]:= S 4 Out[]= i y k 4 0 0{ T Since S n = L n. L n and 8L T n < k, = 8L n <,k, we obtain that i 8S n < i, = 8L n < k=0 i,k 8L n <,k i i + y i i + y. k i { k {
11 694 Edited by Paul Abbott In[6]:= Out[6]= In[7]:= Out[7]= i k=0 i i y i y k k { k k { GHi + + L ÅÅÅÅÅÅÅ GHi + L GH + L FullSimplifyA% i i + y i i + y E k i { k { True The inverse of S n has integer entries. In[8]:= S - Out[8]= i y k { It can be shown that S n - n -. S n. n. That is S n - is similar to S n (see mathworld.wolfram.com/similarmatrices.html). Verify that S n - n -. S n. n for n =. In[9]:= S - -. S. Out[9]= True Since similar matrices have the same eigenvalues, the eigenvalues of S n must come in reciprocal pairs, l and l -. Moreover, if n is odd, then the middle eigenvalue must be unity. Here are the eigenvalues of S. In[0]:= Out[0]= Eigenvalues@S D 94 +!!!,, 4 -!!! = Confirm that the eigenvalues come in reciprocal pairs. In[]:= ê % Out[]= 9 In[]:= Out[]= 4 +!!!,, 4 -!!! = RootReduce@%D 94 -!!!,, 4 +!!! =
12 Tricks of the Trade 69 Generalied Padé Approximation It is often the case that, while an exact solution is difficult or impossible, series or asymptotic solutions to a particular problem can be obtained without too much difficulty. As is well known, Padé approximations (see mathworld.wolfram.com/ PadeApproximant.html) are usually superior to Taylor expansions when functions contain poles, because the use of rational functions allows them to be well represented. Define P N HxL and Q M HxL as polynomials of degrees N and M, with Q M H0L =. In[]:= P n_ Hx_L := p i x i i=0 In[]:= Q m_ Hx_L := + i= n m qi x i The coefficients in MD Padé approximant to f HxL are determined by requiring that the Taylor expansion of f HxL - P N HxLêQ M HxL is of order OHx N +M+ L. Suppose that you have obtained the following Taylor series up to OHx L for a function of interest. In[]:= taylor@x_d = - ÄÄÄÄÄÄÄÄÄÄ x + ÄÄÄÄÄÄÄÄÄÄÄÄÄ 4 x - ÄÄÄÄÄÄÄÄÄÄÄÄ 8 x x4 ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ 94 ; Here is D Padé approximant for this series. In[4]:= << Calculus` In[]:= pade@x_d = Pade@taylor@xD, 8x, 0,, <D Out[]= x - ÅÅÅÅÅÅ x ÅÅÅÅ 4 x ÅÅÅÅÅÅÅÅÅ + ÅÅÅÅÅÅÅ 4 x Verify that the Taylor expansion of f HxL - P HxLêQ HxL is of order OHx L. In[6]:= taylor@xd - pade@xd + O@xD Out[6]= OHx L A plot shows that the Taylor series ( ) diverges from the Padé approximant ( ) near x =.
13 696 Edited by Paul Abbott In[7]:= small = PlotA8taylor@xD, pade@xd<, 8x, 0,<, PlotStyle Æ 9Hue@D, HueA ÄÄÄÄÄ E=E Suppose now that you have also computed the following asymptotic series for the same function. In[8]:= asymptotic@x_d = ÄÄÄÄÄÄÄÄÄÄ x + ÄÄÄÄÄÄÄÄÄÄÄÄÄ 4 x + ÄÄÄÄÄÄÄÄÄÄÄÄ 8 x ; A plot shows that the Padé approximant about x = 0 ( ) and the asymptotic expansion ( ) intersect for x º.. In[9]:= << Graphics` In[0]:= large = DisplayTogetherA PlotAasymptotic@xD, 8x,,0<, PlotStyle Æ HueA ÄÄÄÄÄÄ EE, smalle
14 Tricks of the Trade 697 How can you incorporate information from expansions of a function about two or more points? One approach is to generalie the idea of computing Padé approximants []. Suppose f HxL has the asymptotic expansions f HxL ~ a i Hx - x 0 L i, x Ø x 0, i=0 f HxL ~ b i Hx - x L i, x Ø x, i=0 in the neighborhoods of distinct points x 0 and x, respectively. To compute a two-point Padé approximant to f HxL, we require that the rational function P N HxLêQ M HxL agrees with f HxL up to order J about x 0 and to order K about x, where J + K = N + M +. To use all the information at hand, note that for x 0 = 0, J =, and for x =, K = 4. In[]:= = Exponent@taylor@xD, xd + Out[]= In[]:= k = Exponent@asymptotic@xD, ê xd + Out[]= 4 The diagonal Padé approximant then has N = M = 4. In[]:= n = m = ÄÄÄÄÄ H + k - L Out[]= 4 It is straightforward to determine the coefficients p i and q i in terms of a i and b i. Here we obtain these coefficients using series arithmetic via Solve. In[4]:= SolveA 9taylor@xD - P n HxL ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ Q m HxL + O@xD ä 0, asymptotic@xd - P n HxL ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ Q m HxL + O@x, Dk ä 0=, Join@Table@p i, 8i, 0,n<D, Table@q i, 8i, m<dde Out[4]= 99p Ø ÅÅÅÅÅÅÅÅÅÅÅ 94, p Ø 8 94, p Ø ÅÅÅÅÅ 9, p 4 Ø 0, p 0 Ø, q Ø ÅÅÅÅÅÅÅÅÅ 8, q Ø ÅÅÅÅÅÅÅÅÅÅÅ, q Ø ÅÅÅÅÅ 8 9, q 4 Ø ÅÅÅÅÅÅÅÅÅÅÅ 6 94 == Here is the two-point Padé approximant to f HxL. In[]:= twopoint@x_d = P n HxL ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ ê. First@%D Q m HxL Out[]= 8 x ÅÅÅÅÅÅÅÅÅ x + ÅÅÅÅÅÅÅÅÅÅÅ + ÅÅÅÅÅÅÅ x ÅÅÅÅÅÅÅÅ 6 x 4 x + + ÅÅÅÅÅÅÅÅÅ 8 x + ÅÅÅÅÅÅÅ 8 x 94 ÅÅÅÅÅÅÅÅ ÅÅÅÅ ÅÅÅÅÅÅÅÅ + 9
15 698 Edited by Paul Abbott The two-point Padé approximant ( ) agrees with the Taylor series about x = 0 ( ), the one-point Padé approximant about x = 0 ( ), and the asymptotic expansion ( ) about x =. In[6]:= DisplayTogether@Plot@twopoint@xD, 8x, 0, 0<D, larged References [] C. M. Bender and S. A. Orsag, Advanced Mathematical Methods for Scientists and Engineers, New York: McGraw-Hill, 978, pp Paul Abbott School of Physics, M0 University of Western Australia Stirling Highway Crawley WA 6009, Australia tm@physics.uwa.edu.au physics.uwa.edu.au/~paul
Tricks of the Trade Edited by Paul Abbott
The Mathematica Journal Tricks of the Trade Edited by Paul Abbott This is a column of programming tricks and techniques, most of which, we hope, will be contributed by our readers, either directly as submissions
More informationÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅÅ
HW.nb HW #. D scattering (a) The delta function potential requires a discontinuity of the wave function, y' H+0L - y' H-0L = ÅÅÅÅÅÅÅÅÅÅÅÅ m m yh0l, while the wavefunction itself is continuous. SolveA9
More informationApproximate formulas for the distance term in far and near field interferometry Tim Cornwell, NRAO
EVLA Memo 75 Approximate formulas for the distance term in far and near field interferometry Tim Cornwell, NRAO tcornwel@nrao.edu Introduction In many applications in wide field and near field interferometry,
More informationThe Mathematica Journal Tricks of the Trade
The Mathematica Joural Tricks of the Trade Edited by Paul Abbott This is a colum of programmig tricks ad techiques, most of which, we hope, will be cotributed by our readers, either directly as submissios
More information1 Integration of Rational Functions Using Partial Fractions
MTH Fall 008 Essex County College Division of Mathematics Handout Version 4 September 8, 008 Integration of Rational Functions Using Partial Fractions In the past it was far more usual to simplify or combine
More informationSolutions to the 2005 AP Calculus AB Exam Free Response Questions
Solutions to the 00 AP Calculus AB Exam Free Response Questions Louis A. Talman Department of Mathematical & Computer Sciences Metropolitan State College of Denver Problem. FindRootA ÅÅÅÅ 4 + Sin@p xd
More informationThe Mathematica Journal p-adic Arithmetic
The Mathematica Journal p-adic Arithmetic Stany De Smedt The p-adic numbers were introduced by K. Hensel in 1908 in his book Theorie der algebraïschen Zahlen, Leipzig, 1908. In this article we present
More informationThe Pythagorean Means
Pythagorean-Mean-M6.nb The Pythagorean Mean Prof. Dr. J. Ziegenbalg Intitut für Mathematik und Informatik Univerity of Education, Karlruhe electronic mail: homepage: ziegenbalg@ph-karlruhe.de http://www.ph-karlruhe.de/wp/ziegenbalg/
More informationin Trigonometry Name Section 6.1 Law of Sines Important Vocabulary
Name Chapter 6 Additional Topics in Trigonometry Section 6.1 Law of Sines Objective: In this lesson you learned how to use the Law of Sines to solve oblique triangles and how to find the areas of oblique
More informationTime Independent Perturbation Theory
apr_0-may_5.nb: 5/4/04::9:56:8 Time Independent Perturbation Theory Note: In producing a "final" vrsion of these notes I decided to change my notation from that used in class and by Sakurai. In class,
More information8. Differentiation and regularization
43 8. Regularization 8. Differentiation and regularization 8. Regularization Regularization is the technique to make data behave well when an operator is applied to them. Such data could e.g. be functions,
More informationCHAPTER 7: Systems and Inequalities
(Exercises for Chapter 7: Systems and Inequalities) E.7.1 CHAPTER 7: Systems and Inequalities (A) means refer to Part A, (B) means refer to Part B, etc. (Calculator) means use a calculator. Otherwise,
More informationApplications of integrals
ApplicationsofIntegrals.nb Applications of integrals There are many applications of definite integrals and we cannot include all of them in a document intended to be a review. However there are some very
More informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Taylor s Theorem Can often approximate a function by a polynomial The error in the approximation
More informationHW9.nb 1. be a symmetric (parity-even) function, with a close excited state with an anti-symmetric (parity-odd) wave function.
HW9.nb HW #9. Rectngulr Double-Well Potentil Becuse the potentil is infinite for» x» > + b, the wve function should vnish t x = H + bl. For < x < + b nd - - b < x < -, the potentil vnishes nd the wve function
More informationRené Schoof s Algorithm for Determining the Order of the Group of Points on an Elliptic Curve over a Finite Field
René Schoof s Algorithm for Determining the Order of the Group of Points on an Elliptic Curve over a Finite Field John J. McGee Thesis submitted to the faculty of the Virginia Polytechnic Institute and
More information12.4 The Diagonalization Process
Chapter - More Matrix Algebra.4 The Diagonalization Process We now have the background to understand the main ideas behind the diagonalization process. Definition: Eigenvalue, Eigenvector. Let A be an
More informationMiller Objectives Alignment Math
Miller Objectives Alignment Math 1050 1 College Algebra Course Objectives Spring Semester 2016 1. Use algebraic methods to solve a variety of problems involving exponential, logarithmic, polynomial, and
More informationExamples of spontaneity in terms of increased spatial arrangements Notes on General Chemistry
Examples of spontaneity in terms of increased spatial arrangements Notes on General Chemistry http://quantum.bu.edu/notes/generalchemistry/spontaneityintermsofincreasedspatialarrangements.pdf Last updated
More informationQuadratic Programming
QuadProgram.nb 1 Revised: 5 April, 005 Quadratic Programming Quadratic Programs The general quadratic program proposes to minimie an obective function of the form: the linear constraints: A.x == b, x 0
More informationAlgebraic. techniques1
techniques Algebraic An electrician, a bank worker, a plumber and so on all have tools of their trade. Without these tools, and a good working knowledge of how to use them, it would be impossible for them
More informationchapter 12 MORE MATRIX ALGEBRA 12.1 Systems of Linear Equations
chapter MORE MATRIX ALGEBRA GOALS In Chapter we studied matrix operations and the algebra of sets and logic. We also made note of the strong resemblance of matrix algebra to elementary algebra. The reader
More informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
More informationDamped Rolling Sphere II
DapedRollingSphere(II).nb Daped Rolling Sphere II 008 October 30 printed 008 October 30 This adds daping to the case of a rolling sphere solved using the constraint transforation. Alost all of this notebook
More informationà FIELD EFFECT TRANSISTORS
Prof.M.G.Guvench à FIELD EFFECT TRANSISTORS ü FET: CONTENTS Principles of Operation Models: DC, S.S.A.C. and SPICE Applications: AC coupled S.S. Amplifiers ü FET: NAMES JFET Junction Field Effect Transistor
More informationdegree -6x 3 + 5x 3 Coefficients:
Date P3 Polynomials and Factoring leading coefficient degree -6 3 + 5 3 constant term coefficients Degree: the largest sum of eponents in a term Polynomial: a n n + a n-1 n-1 + + a 1 + a 0 where a n 0
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
More informationLagrange Multiplier Method & Karush-Kuhn-Tucker (KKT) Conditions
LagrangeFunctionKKT.nb 1 Lagrange Multiplier Method & Karush-Kuhn-Tucker (KKT) Conditions General Non-Linear Program Minimize : f@x 1, x,..., x n D # variables = n Subject to : Equality constraints : h
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationNATIONAL BOARD FOR HIGHER MATHEMATICS. M. A. and M.Sc. Scholarship Test. September 22, Time Allowed: 150 Minutes Maximum Marks: 30
NATIONAL BOARD FOR HIGHER MATHEMATICS M. A. and M.Sc. Scholarship Test September 22, 2012 Time Allowed: 150 Minutes Maximum Marks: 30 Please read, carefully, the instructions on the following page 1 INSTRUCTIONS
More informationchapter 11 ALGEBRAIC SYSTEMS GOALS
chapter 11 ALGEBRAIC SYSTEMS GOALS The primary goal of this chapter is to make the reader aware of what an algebraic system is and how algebraic systems can be studied at different levels of abstraction.
More informationSMSU Mathematics Course Content
Southwest Minnesota State University Department of Mathematics SMSU Mathematics Course Content 2012-2013 Thefollowing is a list of possibletopics and techniques to cover in your SMSU College Now course.
More informationMatrix & Linear Algebra
Matrix & Linear Algebra Jamie Monogan University of Georgia For more information: http://monogan.myweb.uga.edu/teaching/mm/ Jamie Monogan (UGA) Matrix & Linear Algebra 1 / 84 Vectors Vectors Vector: A
More informationCHAPTER 8: Matrices and Determinants
(Exercises for Chapter 8: Matrices and Determinants) E.8.1 CHAPTER 8: Matrices and Determinants (A) means refer to Part A, (B) means refer to Part B, etc. Most of these exercises can be done without a
More informationReview for Exam Find all a for which the following linear system has no solutions, one solution, and infinitely many solutions.
Review for Exam. Find all a for which the following linear system has no solutions, one solution, and infinitely many solutions. x + y z = 2 x + 2y + z = 3 x + y + (a 2 5)z = a 2 The augmented matrix for
More informationà 10. DC (DIRECT-COUPLED) AMPLIFIERS
0.DC-Amps-X.nb à 0. DC (DIRECT-COUPLED) AMPLIFIERS ü AC COUPLED SMALL SIGNAL AMPLIFIERS ADVANTAGES:. Signal, load and the amplifier bias are separate. One can work on the bias calculations stage by stage
More informationNATIONAL BOARD FOR HIGHER MATHEMATICS. M. A. and M.Sc. Scholarship Test. September 25, Time Allowed: 150 Minutes Maximum Marks: 30
NATIONAL BOARD FOR HIGHER MATHEMATICS M. A. and M.Sc. Scholarship Test September 25, 2010 Time Allowed: 150 Minutes Maximum Marks: 30 Please read, carefully, the instructions on the following page 1 INSTRUCTIONS
More information8 - Series Solutions of Differential Equations
8 - Series Solutions of Differential Equations 8.2 Power Series and Analytic Functions Homework: p. 434-436 # ü Introduction Our earlier technques allowed us to write our solutions in terms of elementary
More informationA Quick Introduction to Vectors, Vector Calculus, and Coordinate Systems
Copyright 2006, David A. Randall Revised Wed, Feb 06, :36:23 A Quick Introduction to Vectors, Vector Calculus, and Coordinate Systems David A. Randall Department of Atmospheric Science Colorado State University,
More informationTennessee s State Mathematics Standards Precalculus
Tennessee s State Mathematics Standards Precalculus Domain Cluster Standard Number Expressions (N-NE) Represent, interpret, compare, and simplify number expressions 1. Use the laws of exponents and logarithms
More information1 Solutions to selected problems
Solutions to selected problems Section., #a,c,d. a. p x = n for i = n : 0 p x = xp x + i end b. z = x, y = x for i = : n y = y + x i z = zy end c. y = (t x ), p t = a for i = : n y = y(t x i ) p t = p
More informationAppendix A Taylor Approximations and Definite Matrices
Appendix A Taylor Approximations and Definite Matrices Taylor approximations provide an easy way to approximate a function as a polynomial, using the derivatives of the function. We know, from elementary
More informationDifference Equations. Representation in "closed form" Formula solutions
Difference-Equations-3-ClosedForm.nb Difference Equations - Part 3 - Representation in "closed form" Formula solutions Prof. Dr. J. Ziegenbalg Institut für Mathemati und Informati Pädagogische Hochschule
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2018 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationPreface. 2 Linear Equations and Eigenvalue Problem 22
Contents Preface xv 1 Errors in Computation 1 1.1 Introduction 1 1.2 Floating Point Representation of Number 1 1.3 Binary Numbers 2 1.3.1 Binary number representation in computer 3 1.4 Significant Digits
More informationCheck boxes of Edited Copy of Sp Topics (was 217-pilot)
Check boxes of Edited Copy of 10024 Sp 11 213 Topics (was 217-pilot) College Algebra, 9th Ed. [open all close all] R-Basic Algebra Operations Section R.1 Integers and rational numbers Rational and irrational
More informationME 406 Assignment #9 Solutions
ME 06 Assignment #9 Solutions sysid Mathematica 6.0., DynPac.0, ê8ê009 intreset; plotreset; imsize = 50; Problem The equations are x = -x, y = x +y, z = x + z. Let F = x + y +z. To show that the surface
More informationGuideBook for Symbolics
Exercises of The Mathematica GuideBook for Symbolics This document was downloaded from http://www.mathematicaguidebooks.org. It contains extracted material from the notebooks accompanying The Mathematica
More informationax 2 + bx + c = 0 where
Chapter P Prerequisites Section P.1 Real Numbers Real numbers The set of numbers formed by joining the set of rational numbers and the set of irrational numbers. Real number line A line used to graphically
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2016 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationSTATE COUNCIL OF EDUCATIONAL RESEARCH AND TRAINING TNCF DRAFT SYLLABUS.
STATE COUNCIL OF EDUCATIONAL RESEARCH AND TRAINING TNCF 2017 - DRAFT SYLLABUS Subject :Mathematics Class : XI TOPIC CONTENT Unit 1 : Real Numbers - Revision : Rational, Irrational Numbers, Basic Algebra
More informationReference Material /Formulas for Pre-Calculus CP/ H Summer Packet
Reference Material /Formulas for Pre-Calculus CP/ H Summer Packet Week # 1 Order of Operations Step 1 Evaluate expressions inside grouping symbols. Order of Step 2 Evaluate all powers. Operations Step
More informationThe wave equation. The basic equation and it's fundamental solutions on bounded domains. u t t = c 2 u x x. An amalgam of sections 1.5 and 2.2.
The wave equation An amalgam of sections.5 and.. The basic equation and it's fundamental solutions on bounded domains Suppose that uhx, tl describes the displacement from equlibrium of a vibrating string
More informationLinear Algebra: Lecture Notes. Dr Rachel Quinlan School of Mathematics, Statistics and Applied Mathematics NUI Galway
Linear Algebra: Lecture Notes Dr Rachel Quinlan School of Mathematics, Statistics and Applied Mathematics NUI Galway November 6, 23 Contents Systems of Linear Equations 2 Introduction 2 2 Elementary Row
More informationCompletion Date: Monday February 11, 2008
MATH 4 (R) Winter 8 Intermediate Calculus I Solutions to Problem Set #4 Completion Date: Monday February, 8 Department of Mathematical and Statistical Sciences University of Alberta Question. [Sec..9,
More informationGuide for Ph.D. Area Examination in Applied Mathematics
Guide for Ph.D. Area Examination in Applied Mathematics (for graduate students in Purdue University s School of Mechanical Engineering) (revised Fall 2016) This is a 3 hour, closed book, written examination.
More informationPascal s Triangle on a Budget. Accuracy, Precision and Efficiency in Sparse Grids
Covering : Accuracy, Precision and Efficiency in Sparse Grids https://people.sc.fsu.edu/ jburkardt/presentations/auburn 2009.pdf... John Interdisciplinary Center for Applied Mathematics & Information Technology
More informationFEEDBACK, STABILITY and OSCILLATORS
FEEDBACK, STABILITY and OSCILLATORS à FEEDBACK, STABILITY and OSCILLATORS - STABILITY OF FEEDBACK SYSTEMS - Example : ANALYSIS and DESIGN OF PHASE-SHIFT-OSCILLATORS - Example 2: ANALYSIS and DESIGN OF
More informationWorkshop: Average values and uncertainties Quantum aspects of physical chemistry
Workshop: Average values and uncertainties Quantum aspects of physical chemistry http://quantum.bu.edu/pltl/2/2.pdf Last updated Tuesday, November 15, 2005 12:56:11-05:00 Copyright 2005 Dan Dill (dan@bu.edu)
More informationSolution of systems of 1st order linear ODE's: Use of eigenvector expansion and the utility of finding the diagonal matrix that goes with the problem.
Solution of systems of st order linear ODE's: Use of eigenvector expansion and the utility of finding the diagonal matrix that goes with the problem. Solution of an initial value problem for a system of
More informationSolutions. Part I. The Mathematica Notebook. Section 1. Cell structure: Input/Output. Section 2. Mixing Text and Mathematics
Solutions To derive the maximum benefit from these exercises, make an honest effort to create the Mathematica code before consulting the solution offered here. Part I. The Mathematica Notebook Section.
More informationSection 6.2 Larger Systems of Linear Equations
Section 6.2 Larger Systems of Linear Equations Gaussian Elimination In general, to solve a system of linear equations using its augmented matrix, we use elementary row operations to arrive at a matrix
More informationDerivation of the Kalman Filter
Derivation of the Kalman Filter Kai Borre Danish GPS Center, Denmark Block Matrix Identities The key formulas give the inverse of a 2 by 2 block matrix, assuming T is invertible: T U 1 L M. (1) V W N P
More informationMathematics for Chemists
Mathematics for Chemists MATHEMATICS FOR CHEMISTS D. M. Hirst Department of Molecular Sciences, university of Warwick, Coventry M D. M. Hirst 1976 All rights reserved. No part of this publication may be
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson
JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises
More informationLinear Algebra. Matrices Operations. Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0.
Matrices Operations Linear Algebra Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0 The rectangular array 1 2 1 4 3 4 2 6 1 3 2 1 in which the
More informationChapter 1A -- Real Numbers. iff. Math Symbols: Sets of Numbers
Fry Texas A&M University! Fall 2016! Math 150 Notes! Section 1A! Page 1 Chapter 1A -- Real Numbers Math Symbols: iff or Example: Let A = {2, 4, 6, 8, 10, 12, 14, 16,...} and let B = {3, 6, 9, 12, 15, 18,
More informationSolving Differential Equations under Uncertainty
Slide 1 of 8 IMS 2008 June 20-24th International Mathemati Solving Differential Equations under Uncertainty Guaranteed Numerical Integration of Nonlinear Parametric ODEs Yves Papegay Nacim Ramdani INRIA
More informationCheck boxes of Edited Copy of Sp Topics (was 261-pilot)
Check boxes of Edited Copy of 10023 Sp 11 253 Topics (was 261-pilot) Intermediate Algebra (2011), 3rd Ed. [open all close all] R-Review of Basic Algebraic Concepts Section R.2 Ordering integers Plotting
More information7. Piecewise Polynomial (Spline) Interpolation
- 64-7 Piecewise Polynomial (Spline Interpolation Single polynomial interpolation has two major disadvantages First, it is not computationally efficient when the number of data points is large When the
More informationGauss-Jordan Row Reduction and Reduced Row Echelon Form
Gauss-Jordan Row Reduction and Reduced Row Echelon Form If we put the augmented matrix of a linear system in reduced row-echelon form, then we don t need to back-substitute to solve the system. To put
More informationAP Calculus Testbank (Chapter 9) (Mr. Surowski)
AP Calculus Testbank (Chapter 9) (Mr. Surowski) Part I. Multiple-Choice Questions n 1 1. The series will converge, provided that n 1+p + n + 1 (A) p > 1 (B) p > 2 (C) p >.5 (D) p 0 2. The series
More informationB œ c " " ã B œ c 8 8. such that substituting these values for the B 3 's will make all the equations true
System of Linear Equations variables Ð unknowns Ñ B" ß B# ß ÞÞÞ ß B8 Æ Æ Æ + B + B ÞÞÞ + B œ, "" " "# # "8 8 " + B + B ÞÞÞ + B œ, #" " ## # #8 8 # ã + B + B ÞÞÞ + B œ, 3" " 3# # 38 8 3 ã + 7" B" + 7# B#
More informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
More information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More informationSolving Linear Systems Using Gaussian Elimination
Solving Linear Systems Using Gaussian Elimination DEFINITION: A linear equation in the variables x 1,..., x n is an equation that can be written in the form a 1 x 1 +...+a n x n = b, where a 1,...,a n
More information6 - Theory of Higher-Order Linear Differential Equations
6 - Theory of Higher-Order Linear Differential Equations 6.1 Basic Theory of Linear Differential Equations Homework: p. 325-326 #1, 7, 15, 19 ü Introduction A linear differential equation of order n is
More informationReview. Numerical Methods Lecture 22. Prof. Jinbo Bi CSE, UConn
Review Taylor Series and Error Analysis Roots of Equations Linear Algebraic Equations Optimization Numerical Differentiation and Integration Ordinary Differential Equations Partial Differential Equations
More informationThe Fixed String of Elementary Cellular Automata
The Fixed String of Elementary Cellular Automata Jiang Zhisong Department of Mathematics East China University of Science and Technology Shanghai 200237, China zsjiang@ecust.edu.cn Qin Dakang School of
More informationMCR3U Unit 7 Lesson Notes
7.1 Arithmetic Sequences Sequence: An ordered list of numbers identified by a pattern or rule that may stop at some number or continue indefinitely. Ex. 1, 2, 4, 8,... Ex. 3, 7, 11, 15 Term (of a sequence):
More informationGregory's quadrature method
Gregory's quadrature method Gregory's method is among the very first quadrature formulas ever described in the literature, dating back to James Gregory (638-675). It seems to have been highly regarded
More informationPRECALCULUS BISHOP KELLY HIGH SCHOOL BOISE, IDAHO. Prepared by Kristina L. Gazdik. March 2005
PRECALCULUS BISHOP KELLY HIGH SCHOOL BOISE, IDAHO Prepared by Kristina L. Gazdik March 2005 1 TABLE OF CONTENTS Course Description.3 Scope and Sequence 4 Content Outlines UNIT I: FUNCTIONS AND THEIR GRAPHS
More informationMathematics. EC / EE / IN / ME / CE. for
Mathematics for EC / EE / IN / ME / CE By www.thegateacademy.com Syllabus Syllabus for Mathematics Linear Algebra: Matrix Algebra, Systems of Linear Equations, Eigenvalues and Eigenvectors. Probability
More informationSection Taylor and Maclaurin Series
Section.0 Taylor and Maclaurin Series Ruipeng Shen Feb 5 Taylor and Maclaurin Series Main Goal: How to find a power series representation for a smooth function us assume that a smooth function has a power
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationAN IDENTIFICATION ALGORITHM FOR ARMAX SYSTEMS
AN IDENTIFICATION ALGORITHM FOR ARMAX SYSTEMS First the X, then the AR, finally the MA Jan C. Willems, K.U. Leuven Workshop on Observation and Estimation Ben Gurion University, July 3, 2004 p./2 Joint
More informationNATIONAL BOARD FOR HIGHER MATHEMATICS. M. A. and M.Sc. Scholarship Test. September 17, Time Allowed: 150 Minutes Maximum Marks: 30
NATIONAL BOARD FOR HIGHER MATHEMATICS M. A. and M.Sc. Scholarship Test September 17, 2016 Time Allowed: 150 Minutes Maximum Marks: 30 Please read, carefully, the instructions that follow INSTRUCTIONS TO
More informationchapter 5 INTRODUCTION TO MATRIX ALGEBRA GOALS 5.1 Basic Definitions
chapter 5 INTRODUCTION TO MATRIX ALGEBRA GOALS The purpose of this chapter is to introduce you to matrix algebra, which has many applications. You are already familiar with several algebras: elementary
More informationNumerical Methods with MATLAB
Numerical Methods with MATLAB A Resource for Scientists and Engineers G. J. BÖRSE Lehigh University PWS Publishing Company I(T)P AN!NTERNATIONAL THOMSON PUBLISHING COMPANY Boston Albany Bonn Cincinnati
More informationIntroduction - Motivation. Many phenomena (physical, chemical, biological, etc.) are model by differential equations. f f(x + h) f(x) (x) = lim
Introduction - Motivation Many phenomena (physical, chemical, biological, etc.) are model by differential equations. Recall the definition of the derivative of f(x) f f(x + h) f(x) (x) = lim. h 0 h Its
More informationPostulate 2 [Order Axioms] in WRW the usual rules for inequalities
Number Systems N 1,2,3,... the positive integers Z 3, 2, 1,0,1,2,3,... the integers Q p q : p,q Z with q 0 the rational numbers R {numbers expressible by finite or unending decimal expansions} makes sense
More informationContents. 2 Sequences and Series Approximation by Rational Numbers Sequences Basics on Sequences...
Contents 1 Real Numbers: The Basics... 1 1.1 Notation... 1 1.2 Natural Numbers... 4 1.3 Integers... 5 1.4 Fractions and Rational Numbers... 10 1.4.1 Introduction... 10 1.4.2 Powers and Radicals of Rational
More informationPrentice Hall CME Project, Algebra
Prentice Hall Advanced Algebra C O R R E L A T E D T O Oregon High School Standards Draft 6.0, March 2009, Advanced Algebra Advanced Algebra A.A.1 Relations and Functions: Analyze functions and relations
More informationChapter 9: Systems of Equations and Inequalities
Chapter 9: Systems of Equations and Inequalities 9. Systems of Equations Solve the system of equations below. By this we mean, find pair(s) of numbers (x, y) (if possible) that satisfy both equations.
More informationHonors Algebra II Spring 2016 Final Exam Format
Honors Algebra II Spring 2016 Final Exam Format 40 questions, all multiple choice 16 questions from Units 1 4B, 24 questions from Units 5A 6 Can use graphing calculator for entire test May 23 9:38 AM Unit
More information(3) Let Y be a totally bounded subset of a metric space X. Then the closure Y of Y
() Consider A = { q Q : q 2 2} as a subset of the metric space (Q, d), where d(x, y) = x y. Then A is A) closed but not open in Q B) open but not closed in Q C) neither open nor closed in Q D) both open
More informationNATIONAL BOARD FOR HIGHER MATHEMATICS. M. A. and M.Sc. Scholarship Test. September 24, Time Allowed: 150 Minutes Maximum Marks: 30
NATIONAL BOARD FOR HIGHER MATHEMATICS M. A. and M.Sc. Scholarship Test September 24, 2011 Time Allowed: 150 Minutes Maximum Marks: 30 Please read, carefully, the instructions on the following page 1 INSTRUCTIONS
More informationLinear Equations Matrices and Determinants
Linear Equations Matrices and Determinants Introduction The primary purpose of this topic is to develop the mathematical discipline that is necessary to solve a complex problem that requires many numerical
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More information