arxiv: v1 [math.pr] 28 Dec 2018
|
|
- Primrose Stokes
- 5 years ago
- Views:
Transcription
1 Approximating Sepp s constants for te Slepian process Jack Noonan a, Anatoly Zigljavsky a, a Scool of Matematics, Cardiff University, Cardiff, CF4 4AG, UK arxiv:8.0v [mat.pr] 8 Dec 08 Abstract Slepian process S(t) is a stationary Gaussian process wit zero mean and covariance ES(t)S(t ) = max{0, t t }. For any T > 0 and > 0, define F T () = Pr { max t [0,T ] S(t) < } and te constants Λ() = lim T T log F T () and λ() = exp{ Λ()}; we will call tem Sepp s constants. Te aim of te paper is construction of accurate approximations for F T () and ence for te Sepp s constants. We demonstrate tat at least some of te approximations are extremely accurate. Keywords:. Introduction Slepian process, extreme value teory, boundary crossing probability Let S(t), t [0, T ], be a Gaussian process wit mean 0 and covariance ES(t)S(t ) = max{0, t t }. Tis process is often called Slepian process. For any > 0 and x <, define { F T ( x) := Pr max S(t) < } S(0) = x ; (.) t [0,T ] if x we set F T ( x) = 0. Assuming tat x as Gaussian distribution N(0, ), and ence te stationarity of te process S(t), we average F T ( x) and tus define F T () := F T ( x)ϕ(x)dx, were ϕ(x) = (π) / exp{ x /}. Key results on te boundary crossing probabilities for te Slepian process ave been establised by L.Sepp in []. In particular, L.Sepp as derived an explicit formula for F n () wit n integer, see (.) below. As tis explicit formula is quite complicated, in (3.7) in te same paper, L.Sepp as conjectured te existence of te following constant (depending on ) Λ() = lim n n log F n() (.) and raised te question of constructing accurate approximations and bounds for tis constant. Te importance of tis constant is related to te asymptotic relation F T () const[λ()] T as T, (.3) were λ() = exp{ Λ()}. We will call Λ() and λ() Sepp s constants. In tis paper, we are interested in deriving approximations for F T () in te form (.3) and ence for te Sepp s Corresponding autor addresses: NoonanJ@cardiff.ac.uk (Jack Noonan), ZigljavskyAA@cardiff.ac.uk (Anatoly Zigljavsky) Preprint submitted to Elsevier December 3, 08
2 constants. In formulation of approximations, we offer approximations for F T () for all T > and ence approximations for Λ() and λ(). One can apply general results sown in [, 3] and formula (..3) in [4] to approximate λ() but tese results only sow tat λ() as and terefore are of no use ere. On te oter and, te following approximation, also derived from general principles, is simple but more useful. Approximation 0. Poisson Clumping Heuristic, see [5]: F T () exp( ϕ()t ), Λ (0) () = ϕ(), λ (0) () = e ϕ(). For large, tis approximation is quite good, see Tables and below. In Section we derive several new approximations for F T () and λ(). In Section 3 we provide numerical results sowing tat at least some of te approximation derived in Section are extremely accurate. Section 4 contains some minor tecnical details and Section 5 delivers conclusions.. Derivation of approximations for F T () and λ().. Sepp s explicit formula for F n ( x)... Sepp s formula Te following formula is te result (.5) in []: F n ( x) = det ϕ(y i y j+ + ) n i,j=0 dy... dy n+, (.) ϕ(x) D x were T = n is a positive integer, D x = {y,..., y T + x < y < y 3 <... < y n+ }, y 0 = 0, y = x. L.Sepp in [] as also derived explicit formulas for F T ( x) wit non-integral T > 0 but tese formulas are more complicated and are realistically applicable only for small enoug T (say, T 3). From (.) we straigtforwardly obtain F ( ϕ() x) = Φ() Φ(x), (.) ϕ(x) F () = F ( x)ϕ(x)dx = Φ () ϕ() [ Φ() + ϕ() ], were Φ(x) = x ϕ(t)dt. Derivation of explicit formulas for F T ( x) and FT () wit T is relatively easy as te process S(t) is conditionally Markovian in te interval [0, ], see [6]. Formula (.) as been first derived in [7]. In wat follows, F () and F ( 0) (in addition to F () and F ( 0)) play important roles. Te expressions for bot, F () and F ( 0), are more complicated tan expressions for F () and F ( 0). Neverteless, tese expressions can be reduced to a one-dimensional integrals and furter approximated as sown in Appendix; see Section An alternative representation of te Sepp s formula (.) Let T = n be a positive integer, y 0 = 0, y = x. For i = 0,,..., n we set s i = + y i y i+ wit s 0 = x. It follows from Sepp s proof of (.) tat s 0, s,..., s n ave te meaning of te values of te process S(t) at te times t = 0,,..., n: S(i) = s i (i = 0,,..., n). Te range of te variables s i is (, ). Te variables y,..., y n+ are expressed via s 0,..., s n by y k = k s 0 s... s k (k =,..., n+) wit y 0 = 0. Canging te variables, we obtain te following equivalent expression for te probability (.): F n ( x) =... det ϕ(s i + a i,j ) n i,j=0 ds... ds n, (.3) ϕ(x)
3 were 0 for i = j a i,j = y i+ y j+ = (i j) s j+... s i+ for i > j (i j) + s i s j for i < j. Expression (.3) for te probability F n ( x) implies tat te function p(s 0, s,... s n ) = ϕ(s 0 )F n ( s0 ) det ϕ(s i + a i,j ) n i,j=0. (.4) is te joint probability density function for te values S(0), S(),..., S(n) under te condition S(t) < for all t [0, n]. Since s n is te value of S(n), te formula (.4) also sows te transition density from s 0 = x to s n conditionally S(t) < for all t [0, n]: p (n) (x s n) =... det ϕ(s i + a i,j ) n i,j=0 ds... ds n. (.5) ϕ(x) For tis transition density, p(n) (x z)dz = F n( x)... Approximating λ() troug eigenvalues of integral operators... One-step transition In te case n = we obtain from (.5): ( ) p () ϕ(x) ϕ(x +z) (x z) = ϕ(x) det ϕ() ϕ(z) [ = ϕ(z) e ( z)( x)] (.6) wit z = s <. Let λ () be te largest eigenvalue of te te integral operator wit kernel (.6): λ ()p(z) = p(x)p () (x z)dx, z <, were eigenfunction p(x) is some probability density on (, ]. Te Ruelle-Krasnoselskii-Perron- Frobenius teory of bounded linear positive operators (see e.g. Teorem XIII.43 in [8]) implies tat te maximum eigenvalue λ of te operator wit kernel K(x, z) = p () (x z) is simple, real and positive and te eigenfunction p(x) can be cosen as a probability density. Similarly to wat we ave done below in Section.., we can suggest computing good numerical approximations to λ () using Gauss-Legendre quadrature formulas. However, we suggest to use (4.5) from [9] instead; tis elps us to obtain te following simple but rater accurate approximation to λ (): [ ] ˆλ () = Φ() + ϕ()/ ϕ()[ϕ() + Φ()]/ Φ() e / /. Approximation : F T () F () [ˆλ () ] T (T ); Λ () () = log ˆλ (), λ () () = ˆλ ().... Transition in a twice longer interval Consider now te interval [0, ]. We could ave extended te metod of Section.. and used te eigenvalue (square root of it) for te transition s 0 s wit transition density expressed in (.5) wit n =. Tis would improve Approximation but tis improvement is only marginal. Instead, we will use anoter approac: we consider te transition s s but use te interval [0, ] just for setting up te initial condition for observing S(t) at t [, ]. 3
4 For n =, te expression (.4) for te joint probability density function for te values S(0), S(), S() under te condition S(t) < for all t [0, ] as te form ϕ(s 0 ) ϕ(s 0 +s ) ϕ(s 0 +s +s ) p(s 0, s, s ) = ϕ(s 0 )F ( s0 ) det ϕ() ϕ(s ) ϕ(s +s ). ϕ( s ) ϕ() ϕ(s ) Denote by p (z), z <, te non-normalized density of S() under te condition S(t) < for all t [0, ] tat satisfies p (z)dz = F (). Using (.6), we obtain p (z) = p () (x z)ϕ(x)dx = Φ()ϕ(z) Φ(z)ϕ(). Ten te transition density from x = s to z = s under te condition S(t) < for all t [0, ] is acieved by integrating s 0 out and renormalising te joint density: ϕ(s 0 ) ϕ(s 0 +x) ϕ(s 0 +x+z) q (x z) = det ϕ() ϕ(x) ϕ(x+z ) ds 0 p (x) ϕ( x) ϕ() ϕ(z) = Φ()ϕ(x) Φ(x)ϕ() det Φ() Φ(x) Φ(x+z ) ϕ() ϕ(x) ϕ(x+z ) ϕ( x) ϕ() ϕ(z) Let λ () be te largest eigenvalue of te integral operator wit kernel q : λ ()q(z) = q(x)q (x z)dx, z <, were eigenfunction q(x) is some probability density on (, ]. Similarly to te case n =, λ () is simple, real and positive eigenvalue of te operator wit kernel K(x, z) = q (x z) and te eigenfunction q(x) can be cosen as a probability density. In numerical examples below we approximate λ () using te metodology described in [0], p.54. It is based on te Gauss-Legendre discretization of te interval [ C, ], wit some large C > 0, into an N-point set x,..., x N (te x i s are te roots of te N-t Legendre polynomial on [ C, ]), and te use of te Gauss-Legendre weigts w i associated wit points x i ; λ () and q(x) are ten approximated by te largest eigenvalue and associated eigenvector of te matrix D / AD /, were D = diag(w i ), and A i,j = q (x i x j ). If N is large enoug ten te resulting approximation ˆλ () to λ () is arbitrarily accurate.. Approximation : F T () F () [ˆλ () ] T (T ); Λ () () = log ˆλ (), λ () () = ˆλ ()...3. Quality of Approximations and Approximation is more accurate tan Approximation 0 but it is still not accurate enoug. Tis is related to te fact tat te process S(t) is not Markovian and te beaviour of S(t) on te interval [i, i + ] depends on all values of S(t) in te interval [i, i] and not only on te value s i = S(i), wic is a simplification we used for derivation of Approximation. Approximation corrects te bias of Approximation by considering twice longer intervals [i, i + ] and using te beaviour of S(t) in te first alf of te interval [i, i + ] just for setting up te initial condition at [i, i + ]. As sown in Section 3, Approximation is muc more accurate tan Approximations 0 and. Te approximations developed in te following section also carefully consider te dependence of S(t) on its past; tey could be made arbitrarily accurate (on expense of increased computational complexity). 4
5 .3. Furter approximations taking into account te non-markovianity of S(t) As mentioned above, te beaviour of S(t) on te interval [i, i + ] depends on all values of S(t) in te interval [i, i] and not only on te value s i = S(i). Te exact value of te Sepp s constant λ() can be defined as te limit (as i ) of te probability tat S(t) < for all t [i, i + ] under te condition S(t) < for all t i. Using te formula for conditional probability, we obtain λ() = lim i F i ()/F i (). (.7) Waiting a long time witout reacing is not numerically possible and is not wat is really required for computation of λ(). Wat we need is for te process S(t) to (approximately) reac te stationary beaviour in te interval [i, i] under te condition S(t) < for all t < i. Since te memory of S(t) is sort (it follows from te representation S(t) = W (t) W (t + ), were W (t) is te standard Wiener process), tis stationary beaviour of S(t) is practically acieved for very small i, as is also seen from numerical results of Section 3. Moreover, since ratios F i ()/F i () are very close to F i ( 0)/F i ( 0) for i, we can use ratios F i ( 0)/F i ( 0) in (.7) instead. For computing te approximations, it makes integration easier. Te above considerations give rise to several approximations formulated below. We start wit simpler approximations wic are easy to compute and end up wit approximations wic are extremely accurate but are arder to compute. We claim tat for all, Approximation 7 as at least seven correct decimal places as te true value of λ(). However, we would not recommend extremely accurate Approximations 6 and 7 since Approximations 4 and 5 are already very accurate, see Tables and, but are muc easier to compute. Approximation 3, te simplest in te bunc, is also quite accurate. Note tat all approximations for F T () can be applied for any T. [ T Approximation 3: F T () F () λ ()] (3), were λ (3) () = F ( 0)/F ( 0). [ T Approximation 4: F T () F () λ ()] (4), were λ (4) () = F ()/F (). [ T Approximation 5: F T () F () λ ()] (5), were λ (5) () = F 3 ( 0)/F ( 0). [ T 3 Approximation 6: F T () F 3 () λ ()] (6), were λ (6) () = F 3 ()/F (). [ T 4 Approximation 7: F T () F 4 () λ ()] (7), were λ (7) () = F 4 ()/F 3 (). Numerical complexity of tese approximation is related to te necessity of computing eiter F n ( 0) or F n () for suitable n. It follows from (.3) tat F n ( 0) is an n-dimensional integral. Consequently, F n () is an (n + )-dimensional integral. In bot cases, te dimensionality of te integral can be reduced by one, respectively to n and n, wit no furter analytical reduction possible. In view of results of Section 4, computation of Approximations 3 and 4 is easy, computation of Approximation 5 requires numerical evaluation of a one-dimensional integral (wic is not ard) but to compute Approximation 7 we need to approximate a tree-dimensional integral, wic as to be done wit ig precision as oterwise Approximation 7 is not wort using: indeed, Approximations 4 6 are almost as good but are muc easier to compute. As Approximation 7 provides us wit te values wic are practically indistinguisable from te true values of λ(), we use Approximation 7 only for te assessment of te accuracy of oter approximations and do not recommend using it in practice. 3. Numerical results In tis section we discuss te quality of approximations introduced in Sections,. and.3. In Table, we present te values of λ (i) (), i = 0,,..., 7, for a number of different. As mentioned above, λ (7) () is practically te true λ() and terefore we compare all oter approximations against 5
6 λ (7) (). In Table we present te relative errors of all oter approximations against λ (7) (); tat is, te values λ (i) ()/λ (7) () for i = 0,,..., 6. From tese two tables we see tat Approximations -7 are very accurate across all. =0.5 = =.5 = =.5 =3 =3.5 =4 =4.5 λ (0) () λ () () λ () () λ (3) () λ (4) () λ (5) () λ (6) () λ (7) () Table : λ (i) (), i = 0,,..., 7, for different. =0.5 = =.5 = =.5 =3 =3.5 =4 =4.5 λ (0) ().8e e-0.0e-0.0e-0.68e-03.5e-04.4e e e-09 λ () ().4e e-0.07e e-03.6e-03.69e-04.93e-05.88e-06.56e-07 λ () () -3.0e e e e e e e e e-0 λ (3) () -.4e e-04.6e e e e e e-06 -.e-07 λ (4) () -6.48e e e e e e-07.6e-08.35e e- λ (5) () -3.6e-04.49e-04 7.e-05.38e e e e e-.6e- λ (6) () -5.e e e e e e e-.4e- -.63e- Table : Relative errors of λ (i) (), i = 0,,..., 6, against λ (7) (). A plot of te relative errors can be seen in Figure a, were te number next to te line corresponds to te approximation. Approximations,4,6 and 7 are monotonically increasing across all and suggest very accurate lower bounds for te true λ(). Approximations 0 and appear to provide upper bounds for λ() for all. (a) Relative error of λ (i) (), i = 0,..., 6, against λ (7) () (b) λ (0) () (dotted red), λ () () (dased blue) and λ (6) () (solid green) Figure As mentioned in Section..3, Approximation is not as accurate as Approximations 7 because it does not adequately take into account te non-markovianity of S(t). In Figure b we 6
7 ave plotted λ (0) () (dotted red line), λ () () (dased red line) and λ (6) () (solid green line) for a range of interesting. Visually, all λ (i) () wit i =, 4, 5, 6, 7 would be visually indistinguisable from eac oter on te plot in Figure b and λ (3) () would be very close to tem. In Figure a we illustrate te rate of convergence of [log F n ()]/n in (.) to te Sepp s constant Λ() as n increases. Te dotted black lines correspond to simulation results obtained by computing te probability F n () wit 00,000 simulations for =.5, and 3. Te solid red lines correspond to Approximation 4 for Λ() wit te cosen. Figure a demonstrates ow accurate tis computationally ceap approximation is and also demonstrates te importance of te multiplying constant F () in Approximation 4 to correct for te non-linear beaviour seen for small n. In Figure b we investigate te rate of convergence and accuracy of convergence to Λ() using all approximations, were we ave fixed = 3. In tis figure, Approximations 4, 5, 6 and 7 produce results tat are visually indistinguisable to Approximation and ence are not plotted. Terefore, in tis figure, Approximation can be considered as giving te true Sepp s constant Λ(). Te number next to te line corresponds to wic approximation was used. (a) Rate and accuracy of convergence to te Sepp s constant Λ() using simulations (dotted black line) and using Approximation 4 (solid red) for =.5, and 3. (b) Rate and accuracy of convergence to Λ() wit = 3 using all approximations (numbers next to line correspond to wic approximation as been used). Figure 4. Appendix 4.. Simplified form of F () and its approximation Using (.) and canging te order of integration were suitable, F () can expressed troug a one-dimensional integral as follows: F () = Φ() 3 + ϕ() Φ() + ϕ() ϕ()φ()[φ() + ϕ()] [ ( )Φ() + ϕ() ] + 0 Φ(y) ϕ( y)dy Φ( y)ϕ( [ ) Φ( y) / ] dy. Using approximations for Φ(t), it is possible to approximate F () very accurately. For example, 7
8 using te approximation (see []) Φ(t) = { 0.5 exp(0.77t 0.46t ) for t exp( 0.77t 0.46t ) for t > 0, (4.) we obtain F () = Φ() 3 + ϕ() Φ() + ϕ() [ ( )Φ() + ϕ() ] ϕ()φ() [Φ() + ϕ()] + Φ() Φ() 0.5 [ e J(0.96, b, ) π J(.33, b, ) V (.46, b, ) π + V (, b, ) + π π K(.5, b, ) { K(.33, b 3, 0) } ] U(.46, b 4, 0), (4.) π were b = 0.77, b = b 0.77, b =, b 3 = b +.5, b 4 = b +.434, K(x, y, z) = ( ) πe y /(4x) xz y Φ, U(x, y, z) = [ yk(x, y, z) e z(y xz)], (4.3) x x x J(x, y, z) = K(x, y, z) K(x, y, 0) and V (x, y, z) = U(x, y, z) U(x, y, 0). Table 3 sows tat approximation (4.) is very accurate across all of interest. =0.5 = =.5 = =.5 =3 =3.5 =4 =4.5 F () (4.) Table 3: Accuracy of approximation (4.) for F () Using (.), we can express F ( 0) as follows: F ( 0) = Φ() ϕ(0) ϕ()φ() ϕ()φ() + ϕ(0) Using (4.), we can obtain te approximation ϕ(0) ϕ( y)φ( y)ϕ(y)dy Φ( y)ϕ( y)ϕ(y)dy. F ( 0) = Φ() ϕ(0) ϕ()φ() ϕ()φ() + ( ) [ πϕ Φ ( ) Φ ( )] ϕ() e.664 {e.434 [K(.46, , ) K(.46, , )] e.434 [K(.46, , ) K(.46, , )]} [ ] e K(.46, , ) K(.46, , ), (4.4) π were K(x, y, z) can be found in (4.3). Table 4 sows tat approximation (4.4) is very accurate across all of interest. 5. Conclusions In is seminal paper [], L. Sepp derived explicit formulas for F T () = Pr { max t [0,T ] S(t) < }, te distribution of maximum of te so-called Slepian process S(t). As tese explicit formulas are 8
9 =0.5 = =.5 = =.5 =3 =3.5 =4 =4.5 F ( 0) (4.4) Table 4: Accuracy of approximation (4.4) for F ( 0) complicated, in te same paper L. Sepp as introduced a constant Λ() = lim T T log F T () (wic we call Sepp s constant) measuring te rate of decrease of F T () as T grows; L. Sepp also raised te question of constructing accurate approximations and bounds for tis constant. Until now, tis question as not been addressed. To answer it, we ave constructed different approximations for F T () (and ence for Λ()). We ave sown in Section 3 tat at least some of tese approximations are extremely accurate. We ave also provided oter approximations tat are almost as good but are muc simpler to compute. References [] LA Sepp. First passage time for a particular Gaussian process. Te Annals of Matematical Statistics, pages , 97. [] HJ Landau and LA Sepp. On te supremum of a Gaussian process. Sankyā: Te Indian Journal of Statistics, Series A, pages , 970. [3] M Marcus and LA Sepp. Sample beavior of Gaussian processes. In Proc. of te Sixt Berkeley Symposium on Mat. Statist. and Prob, volume, pages 43 4, 97. [4] RJ Adler and J Taylor. Random Fields and Geometry. Springer, 007. [5] D Aldous. Probability Approximations via te Poisson Clumping Heuristic. Springer Science & Business Media, 989. [6] CB Mer and JA McFadden. Certain properties of Gaussian processes and teir first-passage times. Journal of te Royal Statistical Society. Series B (Metodological), pages 505 5, 965. [7] D Slepian. First passage time for a particular Gaussian process. Te Annals of Matematical Statistics, 3():60 6, 96. [8] M Reed and B Simon. Metods of Modern Matematical Pysics: Scattering teory Vol. 3. Academic Press, 979. [9] J Noonan and A Zigljavsky. Approximations of te boundary crossing probabilities for te maximum of moving sums. arxiv preprint arxiv:80.099, 08. [0] JL Moamed and LM Delves. Computational Metods for Integral Equations. Cambridge University Press, 985. [] JT Lin. Approximating te normal tail probability and its inverse for use on a pocket calculator. Applied Statistics, 38():69 70,
Numerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationPrecalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here!
Precalculus Test 2 Practice Questions Page Note: You can expect oter types of questions on te test tan te ones presented ere! Questions Example. Find te vertex of te quadratic f(x) = 4x 2 x. Example 2.
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More information5.1 We will begin this section with the definition of a rational expression. We
Basic Properties and Reducing to Lowest Terms 5.1 We will begin tis section wit te definition of a rational epression. We will ten state te two basic properties associated wit rational epressions and go
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationDerivation Of The Schwarzschild Radius Without General Relativity
Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More informationTime (hours) Morphine sulfate (mg)
Mat Xa Fall 2002 Review Notes Limits and Definition of Derivative Important Information: 1 According to te most recent information from te Registrar, te Xa final exam will be eld from 9:15 am to 12:15
More informationWYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
More informationThe Priestley-Chao Estimator
Te Priestley-Cao Estimator In tis section we will consider te Pristley-Cao estimator of te unknown regression function. It is assumed tat we ave a sample of observations (Y i, x i ), i = 1,..., n wic are
More informationMaterial for Difference Quotient
Material for Difference Quotient Prepared by Stepanie Quintal, graduate student and Marvin Stick, professor Dept. of Matematical Sciences, UMass Lowell Summer 05 Preface Te following difference quotient
More informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationThe Complexity of Computing the MCD-Estimator
Te Complexity of Computing te MCD-Estimator Torsten Bernolt Lerstul Informatik 2 Universität Dortmund, Germany torstenbernolt@uni-dortmundde Paul Fiscer IMM, Danisc Tecnical University Kongens Lyngby,
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationMAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016
MAT244 - Ordinary Di erential Equations - Summer 206 Assignment 2 Due: July 20, 206 Full Name: Student #: Last First Indicate wic Tutorial Section you attend by filling in te appropriate circle: Tut 0
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationAverage Rate of Change
Te Derivative Tis can be tougt of as an attempt to draw a parallel (pysically and metaporically) between a line and a curve, applying te concept of slope to someting tat isn't actually straigt. Te slope
More informationSin, Cos and All That
Sin, Cos and All Tat James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 9, 2017 Outline Sin, Cos and all tat! A New Power Rule Derivatives
More informationTopics in Generalized Differentiation
Topics in Generalized Differentiation J. Marsall As Abstract Te course will be built around tree topics: ) Prove te almost everywere equivalence of te L p n-t symmetric quantum derivative and te L p Peano
More informationPreface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.
Preface Here are my online notes for my course tat I teac ere at Lamar University. Despite te fact tat tese are my class notes, tey sould be accessible to anyone wanting to learn or needing a refreser
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationDefinition of the Derivative
Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of
More informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More informationTHE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225
THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More information2.11 That s So Derivative
2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point
More informationQuantum Numbers and Rules
OpenStax-CNX module: m42614 1 Quantum Numbers and Rules OpenStax College Tis work is produced by OpenStax-CNX and licensed under te Creative Commons Attribution License 3.0 Abstract Dene quantum number.
More information1 1. Rationalize the denominator and fully simplify the radical expression 3 3. Solution: = 1 = 3 3 = 2
MTH - Spring 04 Exam Review (Solutions) Exam : February 5t 6:00-7:0 Tis exam review contains questions similar to tose you sould expect to see on Exam. Te questions included in tis review, owever, are
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationMAT 145. Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points
MAT 15 Test #2 Name Solution Guide Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points Use te grap of a function sown ere as you respond to questions 1 to 8. 1. lim f (x) 0 2. lim
More information2.1 THE DEFINITION OF DERIVATIVE
2.1 Te Derivative Contemporary Calculus 2.1 THE DEFINITION OF DERIVATIVE 1 Te grapical idea of a slope of a tangent line is very useful, but for some uses we need a more algebraic definition of te derivative
More informationMathematics 105 Calculus I. Exam 1. February 13, Solution Guide
Matematics 05 Calculus I Exam February, 009 Your Name: Solution Guide Tere are 6 total problems in tis exam. On eac problem, you must sow all your work, or oterwise torougly explain your conclusions. Tere
More information1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist
Mat 1120 Calculus Test 2. October 18, 2001 Your name Te multiple coice problems count 4 points eac. In te multiple coice section, circle te correct coice (or coices). You must sow your work on te oter
More informationTHE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein
Worksop on Transforms and Filter Banks (WTFB),Brandenburg, Germany, Marc 999 THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS L. Trautmann, R. Rabenstein Lerstul
More information. If lim. x 2 x 1. f(x+h) f(x)
Review of Differential Calculus Wen te value of one variable y is uniquely determined by te value of anoter variable x, ten te relationsip between x and y is described by a function f tat assigns a value
More informationOutline. MS121: IT Mathematics. Limits & Continuity Rates of Change & Tangents. Is there a limit to how fast a man can run?
Outline MS11: IT Matematics Limits & Continuity & 1 Limits: Atletics Perspective Jon Carroll Scool of Matematical Sciences Dublin City University 3 Atletics Atletics Outline Is tere a limit to ow fast
More informationCubic Functions: Local Analysis
Cubic function cubing coefficient Capter 13 Cubic Functions: Local Analysis Input-Output Pairs, 378 Normalized Input-Output Rule, 380 Local I-O Rule Near, 382 Local Grap Near, 384 Types of Local Graps
More informationMATH1151 Calculus Test S1 v2a
MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums
More information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More informationDiscontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System
Discontinuous Galerkin Metods for Relativistic Vlasov-Maxwell System He Yang and Fengyan Li December 1, 16 Abstract e relativistic Vlasov-Maxwell (RVM) system is a kinetic model tat describes te dynamics
More informationQuaternion Dynamics, Part 1 Functions, Derivatives, and Integrals. Gary D. Simpson. rev 01 Aug 08, 2016.
Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals Gary D. Simpson gsim1887@aol.com rev 1 Aug 8, 216 Summary Definitions are presented for "quaternion functions" of a quaternion. Polynomial
More informationFunctions of the Complex Variable z
Capter 2 Functions of te Complex Variable z Introduction We wis to examine te notion of a function of z were z is a complex variable. To be sure, a complex variable can be viewed as noting but a pair of
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationInvestigation of Tangent Polynomials with a Computer Algebra System The AMATYC Review, Vol. 14, No. 1, Fall 1992, pp
Investigation of Tangent Polynomials wit a Computer Algebra System Te AMATYC Review, Vol. 14, No. 1, Fall 199, pp. -7. Jon H. Matews California State University Fullerton Fullerton, CA 9834 By Russell
More informationPolynomial Functions. Linear Functions. Precalculus: Linear and Quadratic Functions
Concepts: definition of polynomial functions, linear functions tree representations), transformation of y = x to get y = mx + b, quadratic functions axis of symmetry, vertex, x-intercepts), transformations
More informationSolve exponential equations in one variable using a variety of strategies. LEARN ABOUT the Math. What is the half-life of radon?
8.5 Solving Exponential Equations GOAL Solve exponential equations in one variable using a variety of strategies. LEARN ABOUT te Mat All radioactive substances decrease in mass over time. Jamie works in
More informationNotes on wavefunctions II: momentum wavefunctions
Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles
More informationNADARAYA WATSON ESTIMATE JAN 10, 2006: version 2. Y ik ( x i
NADARAYA WATSON ESTIMATE JAN 0, 2006: version 2 DATA: (x i, Y i, i =,..., n. ESTIMATE E(Y x = m(x by n i= ˆm (x = Y ik ( x i x n i= K ( x i x EXAMPLES OF K: K(u = I{ u c} (uniform or box kernel K(u = u
More informationStationary Gaussian Markov Processes As Limits of Stationary Autoregressive Time Series
Stationary Gaussian Markov Processes As Limits of Stationary Autoregressive Time Series Lawrence D. Brown, Pilip A. Ernst, Larry Sepp, and Robert Wolpert August 27, 2015 Abstract We consider te class,
More informationBrazilian Journal of Physics, vol. 29, no. 1, March, Ensemble and their Parameter Dierentiation. A. K. Rajagopal. Naval Research Laboratory,
Brazilian Journal of Pysics, vol. 29, no. 1, Marc, 1999 61 Fractional Powers of Operators of sallis Ensemble and teir Parameter Dierentiation A. K. Rajagopal Naval Researc Laboratory, Wasington D. C. 2375-532,
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationEOQ and EPQ-Partial Backordering-Approximations
USING A ONSTANT RATE TO APPROXIMATE A LINEARLY HANGING RATE FOR THE EOQ AND EPQ WITH PARTIAL BAKORDERING David W. Pentico, Palumo-Donaue Scool of Business, Duquesne University, Pittsurg, PA 158-18, pentico@duq.edu,
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationKernel Density Based Linear Regression Estimate
Kernel Density Based Linear Regression Estimate Weixin Yao and Zibiao Zao Abstract For linear regression models wit non-normally distributed errors, te least squares estimate (LSE will lose some efficiency
More informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
More informationAnalytic Functions. Differentiable Functions of a Complex Variable
Analytic Functions Differentiable Functions of a Complex Variable In tis capter, we sall generalize te ideas for polynomials power series of a complex variable we developed in te previous capter to general
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
INTODUCTION ND MTHEMTICL CONCEPTS PEVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips of sine,
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationDifferential equations. Differential equations
Differential equations A differential equation (DE) describes ow a quantity canges (as a function of time, position, ) d - A ball dropped from a building: t gt () dt d S qx - Uniformly loaded beam: wx
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
Capter 1 INTRODUCTION ND MTHEMTICL CONCEPTS PREVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationCS522 - Partial Di erential Equations
CS5 - Partial Di erential Equations Tibor Jánosi April 5, 5 Numerical Di erentiation In principle, di erentiation is a simple operation. Indeed, given a function speci ed as a closed-form formula, its
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More informationQuantum Mechanics Chapter 1.5: An illustration using measurements of particle spin.
I Introduction. Quantum Mecanics Capter.5: An illustration using measurements of particle spin. Quantum mecanics is a teory of pysics tat as been very successful in explaining and predicting many pysical
More informationSolution. Solution. f (x) = (cos x)2 cos(2x) 2 sin(2x) 2 cos x ( sin x) (cos x) 4. f (π/4) = ( 2/2) ( 2/2) ( 2/2) ( 2/2) 4.
December 09, 20 Calculus PracticeTest s Name: (4 points) Find te absolute extrema of f(x) = x 3 0 on te interval [0, 4] Te derivative of f(x) is f (x) = 3x 2, wic is zero only at x = 0 Tus we only need
More information, meant to remind us of the definition of f (x) as the limit of difference quotients: = lim
Mat 132 Differentiation Formulas Stewart 2.3 So far, we ave seen ow various real-world problems rate of cange and geometric problems tangent lines lead to derivatives. In tis section, we will see ow to
More informationPractice Problem Solutions: Exam 1
Practice Problem Solutions: Exam 1 1. (a) Algebraic Solution: Te largest term in te numerator is 3x 2, wile te largest term in te denominator is 5x 2 3x 2 + 5. Tus lim x 5x 2 2x 3x 2 x 5x 2 = 3 5 Numerical
More informationRecall from our discussion of continuity in lecture a function is continuous at a point x = a if and only if
Computational Aspects of its. Keeping te simple simple. Recall by elementary functions we mean :Polynomials (including linear and quadratic equations) Eponentials Logaritms Trig Functions Rational Functions
More informationChapter Seven The Quantum Mechanical Simple Harmonic Oscillator
Capter Seven Te Quantum Mecanical Simple Harmonic Oscillator Introduction Te potential energy function for a classical, simple armonic oscillator is given by ZÐBÑ œ 5B were 5 is te spring constant. Suc
More informationNON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE
Global and Stocastic Analysis Vol. 4 No. 1, January 2017, 1-10 NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE K. PHANEENDRA AND E. SIVA PRASAD Abstract.
More informationSFU UBC UNBC Uvic Calculus Challenge Examination June 5, 2008, 12:00 15:00
SFU UBC UNBC Uvic Calculus Callenge Eamination June 5, 008, :00 5:00 Host: SIMON FRASER UNIVERSITY First Name: Last Name: Scool: Student signature INSTRUCTIONS Sow all your work Full marks are given only
More informationMAT 1339-S14 Class 2
MAT 1339-S14 Class 2 July 07, 2014 Contents 1 Rate of Cange 1 1.5 Introduction to Derivatives....................... 1 2 Derivatives 5 2.1 Derivative of Polynomial function.................... 5 2.2 Te
More informationNew families of estimators and test statistics in log-linear models
Journal of Multivariate Analysis 99 008 1590 1609 www.elsevier.com/locate/jmva ew families of estimators and test statistics in log-linear models irian Martín a,, Leandro Pardo b a Department of Statistics
More informationLogarithmic functions
Roberto s Notes on Differential Calculus Capter 5: Derivatives of transcendental functions Section Derivatives of Logaritmic functions Wat ou need to know alread: Definition of derivative and all basic
More informationAMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.
Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error
More informationReflection Symmetries of q-bernoulli Polynomials
Journal of Nonlinear Matematical Pysics Volume 1, Supplement 1 005, 41 4 Birtday Issue Reflection Symmetries of q-bernoulli Polynomials Boris A KUPERSHMIDT Te University of Tennessee Space Institute Tullaoma,
More information