The Euler and trapezoidal stencils to solve d d x y x = f x, y x
|
|
- Paulina Stevens
- 5 years ago
- Views:
Transcription
1 restart; Te Euler and trapezoidal stencils to solve d d x y x = y x Te purpose of tis workseet is to derive te tree simplest numerical stencils to solve te first order d equation y x d x = y x, and study some of teir properties. We first write down te ODE using te D notation: ode := D(y)(x) = f(x,y(x)); ode := D y x = y x In tis expression, f is assumed to be a known function of te independant variable x and te function tat we are trying to solve for y x. Te simplest numerical stencils to solve tis equation will give us an approximation to y at some point x = X C given some knowledge of y at x = X. All of tese stencils are based on te Taylor series approximation for y x about x = X to linear order: eq1 := y(x) = series(y(x),x=x,); eq1 := y x = y X C D y X x K X C O x K X Let us define as te difference between x and X, and ten get rid of x in te above expression: eq := = x - X; eq3 := subs(isolate(eq,x),eq1); eq := = x K X eq3 := y C X = y X C D y X C O Now, we can remove te first derivative of y by making use of te differential equation: eq4 := subs(x=x,ode); eq5 := subs(eq4,eq3); eq4 := D y X = f X, y X eq5 := y C X = y X C f X, y X C O Te forward Euler algoritm involves discarding te terms of order and iger in tis expression and ence obtain an approximation to y X C in terms of y X. We can accompis tis by converting te above series into a polynomial using te convert/polynom command. Also, we can take X = x i as te i t point on an evenly spaced lattice x i = x 0 C i, were is te lattice spacing and i is an integer; it ten follows tat X C = x i C 1. Furtermore, we label te numeric approximation of y x at x = x i as y i z y x i. Implementing tese steps and notational canges: eq6 := convert(eq5,polynom); eq7 := [y(+x)=y[i+1],y(x)=y[i],x=x[i]]; eq8 := subs(eq7,eq6); eq6 := y C X = y X C f X, y X eq7 := y C X = y i C 1, y X = y i, X = x i (1) () (3) (4) eq8 := y i C 1 = y i C f x i, y i (5) Te last expression is te forward Euler stencil for solving ode. It is called an explicit algoritm because te quantity we wis to calculate y i C 1 is given explicitly in terms of y i. Te backward Euler stencil is
2 obtained in a similar fasion, except we identify, y i, and y i C 1 differently: eq9 := = X - x; eq10 := subs(isolate(eq9,x),eq1); eq11 := subs(eq4,eq10); eq1 := convert(eq11,polynom); eq13 := [y(-+x)=y[i],y(x)=y[i+1],x=x[i+1]]; eq14 := subs(eq13,eq1); eq9 := = X K x eq10 := y K C X = y X K D y X C O K eq11 := y K C X = y X K f X, y X C O K eq1 := y K C X = y X K f X, y X eq13 := y K C X = y i, y X = y i C 1, X = x i C 1 eq14 := y i = y i C 1 K f x i C 1 (6) Tis stencil is explicit because it will not in general be possible to alegraically isolate wat we want to calculate, i.e. y i C 1 z y x i C 1, except for very specific forms of f. We can rearrange te stencils to better demonstrate teir geometric meaning: ForwardEuler := collect(expand(isolate(eq8,f(x[i],y[i]))),); BackwardEuler := collect(expand(isolate(eq14,f(x[i+1],y[i+1]))),); ForwardEuler := f x i, y i BackwardEuler := f x i C 1 Tese forms illustrate tat te forward Euler stencil (first equation) is a simple approximation of te first derivative of y x in te interval x x i, x i C 1 using ode evaluated at te leftand side of te interval. Conversely, te backward Euler stencil uses te differential equation to evalulate te derivative at te rigtand side of te interval. Tere is anoter stencil we can derive tat is te average of tese two approaces: Trapezoidal := (ForwardEuler+BackwardEuler)/; Trapezoidal := 1 f x i, y i C 1 f x i C 1 Tis is called te trapezoidal metod, and it uses ode evaluated at bot ends of te interval to approximate y' x. Like te backward Euler stencil, it represents an implicit sceme. As an example, we can look at wat tese stencils look like for te special case of y = ly were λ is a constant: f := (x,y) - lambda*y; Stencils := [ForwardEuler,BackwardEuler,Trapezoidal]; Labels := ["Forward Euler","Backward Euler","Trapezoidal"]; f := x, y /l y Stencils := l y i Labels :=, l y i C 1, 1 l y i C 1 l y i C 1 "Forward Euler", "Backward Euler", "Trapezoidal" For tis special case, we see tat it is possible to isolate y i C 1 for eac stencil. Let's do tis using a loop: (7) (8) (9)
3 for j from 1 to 3 do: Stencils[j] := factor(isolate(stencils[j],y[i+1])): print(labels[j],stencils[j]): "Forward Euler" = y i l C 1 "Backward Euler" = K l K 1 "Trapezoidal" = K y i Now, lets go back to te general case by undefining f: f := 'f'; Stencils := [ForwardEuler,BackwardEuler,Trapezoidal]; f := f Stencils := f x i, y i y i l C l K, f x i C 1, 1 f x i, y i (10) (11) C 1 f x i C 1 We now want to determine wat te error is in eac stencil. To do tis, we need to rewrite eac of tem in te standard form y i C 1 Ky i KF x i, x i C 1, y i = 0: for j from 1 to 3 do: Stencils[j] := Stencils[j]* + y[i]; Stencils[j] := (rs-ls)(stencils[j])=0; print(stencillabels[j],stencils[j]): StencilLabels 1 K y i K f x i, y i = 0 StencilLabels K f x i C 1 K y i = 0 StencilLabels 3 K 1 f x i, y i C 1 f x i C 1 K y i = 0 Note tat tese relations define our numeric stencils in tese sense tat tey give exact relations between te approximations y i and y i C 1. But if we replace te approximations wit te exact values of y x tey are supposed to represent, te above will represent only approximate equalities. Tat is, if we make te canges x i / x, x i C 1 / x C, y i / y x / y x C, te leftand sides of te above equations will only be approximately equal to zero. We call te magnitude of te discrepancy te "one step error" or "local error" in te stencil. Hence, te one step error in te various stencils will be: eq15 := [y[i]=y(x),y[i+1]=y(x+),x[i]=x,x[i+1]=x+]; for j from 1 to 3 do: Error[Labels[j]] := subs(eq15,ls(stencils[j])); od; eq15 := y i = y x = y x C, x i = x, x i C 1 = x C Error "Forward Euler" := y x C K y x K y x Error "Backward Euler" := y x C K f x C, y x C K y x (1) (13)
4 Error "Trapezoidal" := y x C K 1 y x C 1 f x C, y x C K y x To estimate te magnitudes of tese errors, we expand eac of te above expression as a power series about = 0 (since we are implicitly assuming is a small quantity): for j from 1 to 3 do: Error[Labels[j]] := series(error[labels[j]],,4); od; Error "Forward Euler" := K y x C D y x C 1 D y x C 1 6 D 3 y x 3 C O 4 (13) Error "Backward Euler" := K y x C D y x C 1 D y x K D 1 y x K D y x D y x C 1 6 D 3 y x K 1 D 1, 1 f x, y x K D 1, y x D y x K 1 D, y x D y x 3 C O 4 y x D y x K 1 D Error "Trapezoidal" := K y x C D y x C 1 D y x K 1 D 1 (14) y x K 1 D y x D y x C 1 6 D 3 y x K 1 4 D 1, 1 y x K 1 D 1, y x D y x K 1 4 D, y x D y x K 1 4 D y x D y x 3 C O 4 Notice tat first, second and tird derivatives of y appear explicitly in tese expressions. Te first derivative terms can be removed by making use of ode. Te second derivative can also be removed by examining te derivative of ode: eq16 := diff(ode,x); eq17 := convert(ode,diff); eq18 := subs(eq17,eq16); eq16 := D y x = D 1 y x C D y x eq17 := d dx y x = y x d dx y x eq18 := D y x = D 1 y x C D y x y x Similarly, te tird derivative can be removed by differentiating eq18: eq19 := subs(eq17,diff(eq18,x)); eq19 := D 3 y x = D 1, 1 y x C D 1, y x y x (15) (16)
5 C D 1, y x C D, y x y x y x C D y x D 1 y x C D y x y x We now substitute our formulae for te derivatives of y into te error expressions: for j from 1 to 3 do: Error[Labels[j]] := subs(eq18,eq19,ode,error[labels[j]]); od; 1 Error "Forward Euler" := D y x C 1 1 D y x y x C 1 6 D 1, 1 y x C 1 6 D 1, y x y x C 1 6 D 1, y x C D, y x y x y x C 1 6 D y x D 1 y x C D y x y x 3 C O 4 Error "Backward Euler" := K 1 D 1 y x K 1 D y x y x C K 1 3 D 1, 1 y x K 5 6 D 1, y x y x C 1 6 D 1, y x C D, y x y x y x K 1 3 D y x D 1 y x C D y x y x K 1 D, y x y x 3 C O 4 Error "Trapezoidal" := K 1 1 D 1, 1 y x K 1 3 D 1, y x y x (17) C 1 6 D 1, y x C D, y x y x y x K 1 1 D y x D 1 y x C D y x y x K 1 4 D, y x y x 3 C O 4 We see tat te local error for te forward and backward scemes is O. Futermore, te leading order terms for tose stencils is te negative of te oter one. Since te trapezoidal sceme is te average of te forward and backward algoritms, tis explains wy te leading order error for te trapeziodal sceme is O 3. In oter words, te trapezoidal sceme is more accurate tan te oter two. Note tat if all we are after is te magnitude of te leading order terms in te errors, we can just expand te above in a low order series: for j from 1 to 3 do: Error[Labels[j]] := series(error[labels[j]],,0); od;
6 Error "Forward Euler" := O Error "Backward Euler" := O Error "Trapezoidal" := O 3 (18) We now turn our attention to actual numerical algoritms employing tese stencils. Since te forward Euler sceme is explicit, it is particularly easy to develop some code for it tat works for arbitrary functions f. It is useful to rewrite te stencil wit y i C 1 isolated on te LHS: eq0 := isolate(stencils[1],y[i+1]); eq0 := y i C 1 = y i C f x i, y i Now we can transcribe tis as a mapping tat takes te stepsize and te old values of y i and x i and returns te new value y i C 1 : y_new := (,x_old,y_old) - y_old + f(x_old,y_old)*; y_new :=, x_old, y_old /y_old C f x_old, y_old We could ave equivalently accomplised tis by using te unapply command, wic converts expression into mappings witout aving to re-write tem manually: y_new := unapply(rs(eq0),,x[i],y[i]); y_new :=, y, y3 /y3 C f y, y3 Here is a procedure EulerSol tat calculates te numerical solution of ode in te interval x 0, 1 once te form of y as been fixed. It's arguments are initial data in te form y 0 = y0 and te number of steps N. Te output is a list of points, wic we compare in te plot to te output generated by dsolve/numeric for te same problem. f := (x,y) - -y^+*x; ode; y0 := ; N := 15; EulerSol := proc(y0,n) local, x, y, i: end proc: := evalf(1/n): x[0] := 0: y[0] := y0: for i from 1 to N do: x[i] := x[i-1] + : y[i] := y_new(,x[i-1],y[i-1]): [seq([x[i],y[i]],i=0..n)]: data := EulerSol(y0,N); dsolvesol := rs(dsolve({ode,y(0)=y0},numeric,output=listprocedure) []); plot([data,dsolvesol(x)],x=0..1,axes=boxed,legend=[`forward Euler`, `dsolve/numeric`]); f := x, y /Ky C x (19) (0) (1)
7 D y x = Ky x C x y0 := N := 15 data := 0,, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , dsolvesol := proc x... end proc x Forward Euler dsolve/numeric Of course, tis plot is really comparing te results of two numeric approximations to te true solution of te problem. It is also useful to compare numeric answers to analytic results (if available). For te above coice of y tere is an analytic solution in terms of Airy functions: analytic_sol := dsolve({ode,y(0)=y0}): analytic_sol := rs(analytic_sol); ()
8 analytic_sol := 1 / 3 K4 p 3 5 / 6 C 3 1 / 3 G 3 4 p 3 1 / 3 C 3 1 / 3 G 3 3 / 3 AiryAi 1, 1 / 3 x 3 1 / 6 C AiryBi 1, 1 / 3 x () K4 p 3 5 / 6 C 3 1 / 3 G 3 4 p 3 1 / 3 C 3 1 / 3 G 3 3 / 3 AiryAi 1 / 3 x 3 1 / 6 C AiryBi Here is some code tat generates a movie comparing ow te numerical solution converges to te analytic solution as te number of cells is increased: for j from 1 to 5 do: N := 4*j: p[j] := plot([eulersol(y0,n),analytic_sol],x=0..1,title = cat (`number of cells N = `,N),axes=boxed,legend=[`Forward Euler solution`,`analytic solution`],style=[point,line]); plots[display](convert(p,list),insequence=true); 1 / 3 x
9 number of cells N = explicitly: x Forward Euler solution analytic solution We would also like to examine te performance of te backward and trapezoidal stencils, so let's make a simple coice of y tat admits an analytic solution allows eac stencil to be solve for y i C 1 f := (x,y) - lambda*y; y0 := 'y0'; ode; analytic_sol := rs(dsolve({ode,y(0)=y0})); for j from 1 to 3 do: Stencils[j] := isolate(stencils[j],y[i+1]); print(labels[j],stencils[j]); f := x, y /l y y0 := y0 D y x = l y x analytic_sol := y0 e l x "Forward Euler" = l y i C y i
10 "Backward Euler" = y i 1 K l "Trapezoidal" = y i C 1 l y i 1 K 1 l (3) As for te above code, it will be useful to realize eac stencil as a mapping wit arguments, l and y i. We can do tis using a loop and te unapply command: for j from 1 to 3 do: stencil[j] := unapply(rs(stencils[j]),,lambda,y[i]); od; stencil 1 :=, l, y3 /l y3 C y3 stencil :=, l, y3 / y3 1 K l y3 C 1 l y3 stencil 3 :=, l, y3 / 1 K 1 l (4) Conversely, we could ave accomplised te same ting in one line by using te map command, wic applies te same mapping onto eac element of a list (in tis case we are applying operations onto eac element of Stencils to generate a new list stencil): stencil := 'stencil': stencil := map(u-unapply(rs(u),,lambda,y[i]),stencils); stencil :=, l, y3 /l y3 C y3,, l, y3 / y3 y3 C 1 l y3,, l, y3 / 1 K l 1 K 1 l It is interesting to note tat even toug eac of te stencils give different expressions for y i C 1, tey actually agree wit one anoter to order. To see tis, let's perform some Taylor series expansions: for j from 1 to 3 do: Labels[j],y[i+1] = series(stencil[j](,lambda,y[i]),,); od; "Forward Euler" = y i C l y i "Backward Euler" = y i C l y i C O "Trapezoidal" = y i C l y i C O Now, ere is a procedure EulerSol tat calculates a numeric solution for y x for x 0, 1 using N cells and assuming y 0 = y0 and a given value of λ. Te value of coice dictates wic stencil is used: 1 for forward Euler, for backward Euler, and 3 for trapezoidal. EulerSol := proc(y0,n,lambda,coice) local, x, y, i: (5) (6)
11 := evalf(1/n): x[0] := 0: y[0] := y0: for i from 1 to N do: x[i] := x[i-1] + : y[i] := stencil[coice](,lambda,y[i-1]): [seq([x[i],y[i]],i=0..n)]: end proc: Here is a plot of te numeric output for eac stencil compared to te actually solution. (Because we are going to plot 4 curves on our grap, I needed to augment Labels into NewLabels in order to gereate te legend.) It seems as if te trapezoidal metod performs muc better tan te oter two: N := 15; y0 := ; lambda := -4; NewLabels := [op(labels),"analytic solution"]; plot([seq(eulersol(y0,n,lambda,j),j=1..3),analytic_sol],x=0..1, legend=newlabels,axes=boxed,style=[point$3,line]); N := 15 NewLabels := y0 := l := K4 "Forward Euler", "Backward Euler", "Trapezoidal", "analytic solution"
12 x Forward Euler Backward Euler Trapezoidal analytic solution We can quantify exactly ow well te various stencils are doing by comparing te numeric and actual values for y at some fixed value o say x = 1. We call tis te global error in te numeric solution at x = 1, wic we denote by ϵ: epsilon := proc(y0,n,lambda,coice) local numerical, actual: numerical := EulerSol(y0,N,lambda,coice)[N+1][]: actual := y0*exp(lambda): evalf(abs(numerical - actual)): end proc: In tis procedure, note tat te [N+1][] suffix after EulerSol(y0,N,lambda,coice) as te effect of picking out te last element of te list (wic is itself a list of quantities), and ten picking out te second element of tat sub-list (wic is te numerical answer for y 1 ). We now generate a log-log plot of te global errors for eac stencil as a function of N: data := 'data': for j from 1 to 3 do: data[j] := [seq([n,epsilon(y0,n,lambda,j)],n= ,5)];
13 plots[loglogplot](convert(data,list),axes=boxed,labels=[`number of steps`,`global error`],legend=labels); global error number of steps Forward Euler Backward Euler Trapezoidal For N T 100, te error curves look linear on tis log-log plot. Tis implies tat above some tresold number of steps, te errors obey a approximate power law e z an b were a and b are constants. We can determine te values of tese parameters by fitting a power law to our error data using Statistics/PowerFit. We first need to regenerate our data for N R 100; i.e., te range for wic we believe a power law ougt to be valid: data := 'data': for j from 1 to 3 do: data[j] := [seq([n,epsilon(y0,n,lambda,j)],n= ,5)]; We ten use te Statistics[PowerFit] command to perform te fit. Note tat tis function requires te input to be a matrix, wic is wy we use te convert/matrix structure. Te raw output of te fitting algoritms is sown for eac stencil, as well as te value of b in te power-law e z an b. for j from 1 to 3 do: fit[j] := Statistics[PowerFit](convert(data[j],Matrix),_N, output=[leastsquaresfunction,parametervalues]): print(labels[j],fit[j],b=fit[j][][]);
14 "Forward Euler", K "Backward Euler", "Trapezoidal", _N , K K _N , K K _N , K K , b =, b = K , b = K We see tat for asyptotically large numbers of steps N [ 1, te global errors are O N K1 = O for te forward and backward Euler stencils, and O N K = O for te trapezoidal metod. Tis confirms te general expectation tat for a stencil wit one step error O p, we expect te global error to obey global error w (number of steps)#(one-step error for eac step) = N # O p = O K1 # O p = O p K 1. [Recall tat we sowed in (18) above tat p = for forward and backward Euler, and p = 3 for te trapezoidal algoirtm.] NOTE: te above syntax for Statistics[PowerFit] only works in Maple 15. For previous versions, you need to do tis: for j from 1 to 3 do: M := convert(data[j],matrix): X := LinearAlgebra[Column](M,1): Y := LinearAlgebra[Column](M,): fit[j] := Statistics[PowerFit](X,Y,_N,output= [leastsquaresfunction,parametervalues]): print(labels[j],fit[j],b=fit[j][][]); "Forward Euler", _N , K , b = K K "Backward Euler", "Trapezoidal", _N , K K _N , K K , b = K , b = K (7) (8)
3.6 Directional Derivatives and the Gradient Vector
288 CHAPTER 3. FUNCTIONS OF SEVERAL VARIABLES 3.6 Directional Derivatives and te Gradient Vector 3.6.1 Functions of two Variables Directional Derivatives Let us first quickly review, one more time, te
More informationSection 2.3: Calculating Limits using the Limit Laws
Section 2.3: Calculating Limits using te Limit Laws In previous sections, we used graps and numerics to approimate te value of a it if it eists. Te problem wit tis owever is tat it does not always give
More information4.1 Tangent Lines. y 2 y 1 = y 2 y 1
41 Tangent Lines Introduction Recall tat te slope of a line tells us ow fast te line rises or falls Given distinct points (x 1, y 1 ) and (x 2, y 2 ), te slope of te line troug tese two points is cange
More information2 The Derivative. 2.0 Introduction to Derivatives. Slopes of Tangent Lines: Graphically
2 Te Derivative Te two previous capters ave laid te foundation for te study of calculus. Tey provided a review of some material you will need and started to empasize te various ways we will view and use
More informationLinear Interpolating Splines
Jim Lambers MAT 772 Fall Semester 2010-11 Lecture 17 Notes Tese notes correspond to Sections 112, 11, and 114 in te text Linear Interpolating Splines We ave seen tat ig-degree polynomial interpolation
More information2.8 The derivative as a function
CHAPTER 2. LIMITS 56 2.8 Te derivative as a function Definition. Te derivative of f(x) istefunction f (x) defined as follows f f(x + ) f(x) (x). 0 Note: tis differs from te definition in section 2.7 in
More informationPiecewise Polynomial Interpolation, cont d
Jim Lambers MAT 460/560 Fall Semester 2009-0 Lecture 2 Notes Tese notes correspond to Section 4 in te text Piecewise Polynomial Interpolation, cont d Constructing Cubic Splines, cont d Having determined
More information4.2 The Derivative. f(x + h) f(x) lim
4.2 Te Derivative Introduction In te previous section, it was sown tat if a function f as a nonvertical tangent line at a point (x, f(x)), ten its slope is given by te it f(x + ) f(x). (*) Tis is potentially
More informationBounding Tree Cover Number and Positive Semidefinite Zero Forcing Number
Bounding Tree Cover Number and Positive Semidefinite Zero Forcing Number Sofia Burille Mentor: Micael Natanson September 15, 2014 Abstract Given a grap, G, wit a set of vertices, v, and edges, various
More informationMATH 5a Spring 2018 READING ASSIGNMENTS FOR CHAPTER 2
MATH 5a Spring 2018 READING ASSIGNMENTS FOR CHAPTER 2 Note: Tere will be a very sort online reading quiz (WebWork) on eac reading assignment due one our before class on its due date. Due dates can be found
More informationHaar Transform CS 430 Denbigh Starkey
Haar Transform CS Denbig Starkey. Background. Computing te transform. Restoring te original image from te transform 7. Producing te transform matrix 8 5. Using Haar for lossless compression 6. Using Haar
More informationMore on Functions and Their Graphs
More on Functions and Teir Graps Difference Quotient ( + ) ( ) f a f a is known as te difference quotient and is used exclusively wit functions. Te objective to keep in mind is to factor te appearing in
More information, 1 1, A complex fraction is a quotient of rational expressions (including their sums) that result
RT. Complex Fractions Wen working wit algebraic expressions, sometimes we come across needing to simplify expressions like tese: xx 9 xx +, xx + xx + xx, yy xx + xx + +, aa Simplifying Complex Fractions
More informationCHAPTER 7: TRANSCENDENTAL FUNCTIONS
7.0 Introduction and One to one Functions Contemporary Calculus 1 CHAPTER 7: TRANSCENDENTAL FUNCTIONS Introduction In te previous capters we saw ow to calculate and use te derivatives and integrals of
More informationMaterials: Whiteboard, TI-Nspire classroom set, quadratic tangents program, and a computer projector.
Adam Clinc Lesson: Deriving te Derivative Grade Level: 12 t grade, Calculus I class Materials: Witeboard, TI-Nspire classroom set, quadratic tangents program, and a computer projector. Goals/Objectives:
More informationLesson 6 MA Nick Egbert
Overview From kindergarten we all know ow to find te slope of a line: rise over run, or cange in over cange in. We want to be able to determine slopes of functions wic are not lines. To do tis we use te
More informationSection 1.2 The Slope of a Tangent
Section 1.2 Te Slope of a Tangent You are familiar wit te concept of a tangent to a curve. Wat geometric interpretation can be given to a tangent to te grap of a function at a point? A tangent is te straigt
More informationFast Calculation of Thermodynamic Properties of Water and Steam in Process Modelling using Spline Interpolation
P R E P R N T CPWS XV Berlin, September 8, 008 Fast Calculation of Termodynamic Properties of Water and Steam in Process Modelling using Spline nterpolation Mattias Kunick a, Hans-Joacim Kretzscmar a,
More informationNumerical Derivatives
Lab 15 Numerical Derivatives Lab Objective: Understand and implement finite difference approximations of te derivative in single and multiple dimensions. Evaluate te accuracy of tese approximations. Ten
More information1.4 RATIONAL EXPRESSIONS
6 CHAPTER Fundamentals.4 RATIONAL EXPRESSIONS Te Domain of an Algebraic Epression Simplifying Rational Epressions Multiplying and Dividing Rational Epressions Adding and Subtracting Rational Epressions
More informationAlternating Direction Implicit Methods for FDTD Using the Dey-Mittra Embedded Boundary Method
Te Open Plasma Pysics Journal, 2010, 3, 29-35 29 Open Access Alternating Direction Implicit Metods for FDTD Using te Dey-Mittra Embedded Boundary Metod T.M. Austin *, J.R. Cary, D.N. Smite C. Nieter Tec-X
More informationOur Calibrated Model has No Predictive Value: An Example from the Petroleum Industry
Our Calibrated Model as No Predictive Value: An Example from te Petroleum Industry J.N. Carter a, P.J. Ballester a, Z. Tavassoli a and P.R. King a a Department of Eart Sciences and Engineering, Imperial
More information1 Finding Trigonometric Derivatives
MTH 121 Fall 2008 Essex County College Division of Matematics Hanout Version 8 1 October 2, 2008 1 Fining Trigonometric Derivatives 1.1 Te Derivative as a Function Te efinition of te erivative as a function
More informationChapter K. Geometric Optics. Blinn College - Physics Terry Honan
Capter K Geometric Optics Blinn College - Pysics 2426 - Terry Honan K. - Properties of Ligt Te Speed of Ligt Te speed of ligt in a vacuum is approximately c > 3.0µ0 8 mês. Because of its most fundamental
More informationFault Localization Using Tarantula
Class 20 Fault localization (cont d) Test-data generation Exam review: Nov 3, after class to :30 Responsible for all material up troug Nov 3 (troug test-data generation) Send questions beforeand so all
More informationHash-Based Indexes. Chapter 11. Comp 521 Files and Databases Spring
Has-Based Indexes Capter 11 Comp 521 Files and Databases Spring 2010 1 Introduction As for any index, 3 alternatives for data entries k*: Data record wit key value k
More informationANTENNA SPHERICAL COORDINATE SYSTEMS AND THEIR APPLICATION IN COMBINING RESULTS FROM DIFFERENT ANTENNA ORIENTATIONS
NTNN SPHRICL COORDINT SSTMS ND THIR PPLICTION IN COMBINING RSULTS FROM DIFFRNT NTNN ORINTTIONS llen C. Newell, Greg Hindman Nearfield Systems Incorporated 133. 223 rd St. Bldg. 524 Carson, C 9745 US BSTRCT
More information19.2 Surface Area of Prisms and Cylinders
Name Class Date 19 Surface Area of Prisms and Cylinders Essential Question: How can you find te surface area of a prism or cylinder? Resource Locker Explore Developing a Surface Area Formula Surface area
More informationMean Shifting Gradient Vector Flow: An Improved External Force Field for Active Surfaces in Widefield Microscopy.
Mean Sifting Gradient Vector Flow: An Improved External Force Field for Active Surfaces in Widefield Microscopy. Margret Keuper Cair of Pattern Recognition and Image Processing Computer Science Department
More informationODE IVP. An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable:
Euler s Methods ODE IVP An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable: The equation is coupled with an initial value/condition (i.e., value
More information12.2 TECHNIQUES FOR EVALUATING LIMITS
Section Tecniques for Evaluating Limits 86 TECHNIQUES FOR EVALUATING LIMITS Wat ou sould learn Use te dividing out tecnique to evaluate its of functions Use te rationalizing tecnique to evaluate its of
More informationProofs of Derivative Rules
Proos o Derivative Rules Ma 16010 June 2018 Altoug proos are not generally a part o tis course, it s never a ba iea to ave some sort o explanation o wy tings are te way tey are. Tis is especially relevant
More informationThe (, D) and (, N) problems in double-step digraphs with unilateral distance
Electronic Journal of Grap Teory and Applications () (), Te (, D) and (, N) problems in double-step digraps wit unilateral distance C Dalfó, MA Fiol Departament de Matemàtica Aplicada IV Universitat Politècnica
More informationMAC-CPTM Situations Project
raft o not use witout permission -P ituations Project ituation 20: rea of Plane Figures Prompt teacer in a geometry class introduces formulas for te areas of parallelograms, trapezoids, and romi. e removes
More informationComparison of the Efficiency of the Various Algorithms in Stratified Sampling when the Initial Solutions are Determined with Geometric Method
International Journal of Statistics and Applications 0, (): -0 DOI: 0.9/j.statistics.000.0 Comparison of te Efficiency of te Various Algoritms in Stratified Sampling wen te Initial Solutions are Determined
More information13.5 DIRECTIONAL DERIVATIVES and the GRADIENT VECTOR
13.5 Directional Derivatives and te Gradient Vector Contemporary Calculus 1 13.5 DIRECTIONAL DERIVATIVES and te GRADIENT VECTOR Directional Derivatives In Section 13.3 te partial derivatives f x and f
More informationHash-Based Indexes. Chapter 11. Comp 521 Files and Databases Fall
Has-Based Indexes Capter 11 Comp 521 Files and Databases Fall 2012 1 Introduction Hasing maps a searc key directly to te pid of te containing page/page-overflow cain Doesn t require intermediate page fetces
More informationCommunicator for Mac Quick Start Guide
Communicator for Mac Quick Start Guide 503-968-8908 sterling.net training@sterling.net Pone Support 503.968.8908, option 2 pone-support@sterling.net For te most effective support, please provide your main
More informationCSE 332: Data Structures & Parallelism Lecture 8: AVL Trees. Ruth Anderson Winter 2019
CSE 2: Data Structures & Parallelism Lecture 8: AVL Trees Rut Anderson Winter 29 Today Dictionaries AVL Trees /25/29 2 Te AVL Balance Condition: Left and rigt subtrees of every node ave eigts differing
More informationVector Processing Contours
Vector Processing Contours Andrey Kirsanov Department of Automation and Control Processes MAMI Moscow State Tecnical University Moscow, Russia AndKirsanov@yandex.ru A.Vavilin and K-H. Jo Department of
More informationInvestigating an automated method for the sensitivity analysis of functions
Investigating an automated metod for te sensitivity analysis of functions Sibel EKER s.eker@student.tudelft.nl Jill SLINGER j..slinger@tudelft.nl Delft University of Tecnology 2628 BX, Delft, te Neterlands
More informationEuler s Methods (a family of Runge- Ku9a methods)
Euler s Methods (a family of Runge- Ku9a methods) ODE IVP An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable: The equation is coupled with an
More informationCE 221 Data Structures and Algorithms
CE Data Structures and Algoritms Capter 4: Trees (AVL Trees) Text: Read Weiss, 4.4 Izmir University of Economics AVL Trees An AVL (Adelson-Velskii and Landis) tree is a binary searc tree wit a balance
More informationMulti-Stack Boundary Labeling Problems
Multi-Stack Boundary Labeling Problems Micael A. Bekos 1, Micael Kaufmann 2, Katerina Potika 1 Antonios Symvonis 1 1 National Tecnical University of Atens, Scool of Applied Matematical & Pysical Sciences,
More informationCS 234. Module 6. October 16, CS 234 Module 6 ADT Dictionary 1 / 33
CS 234 Module 6 October 16, 2018 CS 234 Module 6 ADT Dictionary 1 / 33 Idea for an ADT Te ADT Dictionary stores pairs (key, element), were keys are distinct and elements can be any data. Notes: Tis is
More informationLimits and Continuity
CHAPTER Limits and Continuit. Rates of Cange and Limits. Limits Involving Infinit.3 Continuit.4 Rates of Cange and Tangent Lines An Economic Injur Level (EIL) is a measurement of te fewest number of insect
More informationData Structures and Programming Spring 2014, Midterm Exam.
Data Structures and Programming Spring 2014, Midterm Exam. 1. (10 pts) Order te following functions 2.2 n, log(n 10 ), 2 2012, 25n log(n), 1.1 n, 2n 5.5, 4 log(n), 2 10, n 1.02, 5n 5, 76n, 8n 5 + 5n 2
More informationAnnouncements. Lilian s office hours rescheduled: Fri 2-4pm HW2 out tomorrow, due Thursday, 7/7. CSE373: Data Structures & Algorithms
Announcements Lilian s office ours resceduled: Fri 2-4pm HW2 out tomorrow, due Tursday, 7/7 CSE373: Data Structures & Algoritms Deletion in BST 2 5 5 2 9 20 7 0 7 30 Wy migt deletion be arder tan insertion?
More information12.2 Techniques for Evaluating Limits
335_qd /4/5 :5 PM Page 863 Section Tecniques for Evaluating Limits 863 Tecniques for Evaluating Limits Wat ou sould learn Use te dividing out tecnique to evaluate its of functions Use te rationalizing
More informationMTH-112 Quiz 1 - Solutions
MTH- Quiz - Solutions Words in italics are for eplanation purposes onl (not necessar to write in te tests or. Determine weter te given relation is a function. Give te domain and range of te relation. {(,
More informationPYRAMID FILTERS BASED ON BILINEAR INTERPOLATION
PYRAMID FILTERS BASED ON BILINEAR INTERPOLATION Martin Kraus Computer Grapics and Visualization Group, Tecnisce Universität Müncen, Germany krausma@in.tum.de Magnus Strengert Visualization and Interactive
More informationNOTES: A quick overview of 2-D geometry
NOTES: A quick overview of 2-D geometry Wat is 2-D geometry? Also called plane geometry, it s te geometry tat deals wit two dimensional sapes flat tings tat ave lengt and widt, suc as a piece of paper.
More informationUnsupervised Learning for Hierarchical Clustering Using Statistical Information
Unsupervised Learning for Hierarcical Clustering Using Statistical Information Masaru Okamoto, Nan Bu, and Tosio Tsuji Department of Artificial Complex System Engineering Hirosima University Kagamiyama
More informationYou Try: A. Dilate the following figure using a scale factor of 2 with center of dilation at the origin.
1 G.SRT.1-Some Tings To Know Dilations affect te size of te pre-image. Te pre-image will enlarge or reduce by te ratio given by te scale factor. A dilation wit a scale factor of 1> x >1enlarges it. A dilation
More informationMeasuring Length 11and Area
Measuring Lengt 11and Area 11.1 Areas of Triangles and Parallelograms 11.2 Areas of Trapezoids, Romuses, and Kites 11.3 Perimeter and Area of Similar Figures 11.4 Circumference and Arc Lengt 11.5 Areas
More informationNumerical solution of hybrid fuzzy differential equations by fuzzy neural network
Available online at ttp://ijimsrbiauacir/ Int J Industrial Matematics (ISSN 2008-5621) Vol 6, No 2, 2014 Article ID IJIM-00371, 15 pages Researc Article Numerical solution of ybrid fuzzy differential equations
More informationA Cost Model for Distributed Shared Memory. Using Competitive Update. Jai-Hoon Kim Nitin H. Vaidya. Department of Computer Science
A Cost Model for Distributed Sared Memory Using Competitive Update Jai-Hoon Kim Nitin H. Vaidya Department of Computer Science Texas A&M University College Station, Texas, 77843-3112, USA E-mail: fjkim,vaidyag@cs.tamu.edu
More informationCubic smoothing spline
Cubic smooting spline Menu: QCExpert Regression Cubic spline e module Cubic Spline is used to fit any functional regression curve troug data wit one independent variable x and one dependent random variable
More informationComputing geodesic paths on manifolds
Proc. Natl. Acad. Sci. USA Vol. 95, pp. 8431 8435, July 1998 Applied Matematics Computing geodesic pats on manifolds R. Kimmel* and J. A. Setian Department of Matematics and Lawrence Berkeley National
More informationDensity Estimation Over Data Stream
Density Estimation Over Data Stream Aoying Zou Dept. of Computer Science, Fudan University 22 Handan Rd. Sangai, 2433, P.R. Cina ayzou@fudan.edu.cn Ziyuan Cai Dept. of Computer Science, Fudan University
More informationAn Algorithm for Loopless Deflection in Photonic Packet-Switched Networks
An Algoritm for Loopless Deflection in Potonic Packet-Switced Networks Jason P. Jue Center for Advanced Telecommunications Systems and Services Te University of Texas at Dallas Ricardson, TX 75083-0688
More informationAll truths are easy to understand once they are discovered; the point is to discover them. Galileo
Section 7. olume All truts are easy to understand once tey are discovered; te point is to discover tem. Galileo Te main topic of tis section is volume. You will specifically look at ow to find te volume
More informationTHANK YOU FOR YOUR PURCHASE!
THANK YOU FOR YOUR PURCHASE! Te resources included in tis purcase were designed and created by me. I ope tat you find tis resource elpful in your classroom. Please feel free to contact me wit any questions
More informationParallel Simulation of Equation-Based Models on CUDA-Enabled GPUs
Parallel Simulation of Equation-Based Models on CUDA-Enabled GPUs Per Ostlund Department of Computer and Information Science Linkoping University SE-58183 Linkoping, Sweden per.ostlund@liu.se Kristian
More informationCESILA: Communication Circle External Square Intersection-Based WSN Localization Algorithm
Sensors & Transducers 2013 by IFSA ttp://www.sensorsportal.com CESILA: Communication Circle External Square Intersection-Based WSN Localization Algoritm Sun Hongyu, Fang Ziyi, Qu Guannan College of Computer
More informationAreas of Parallelograms and Triangles. To find the area of parallelograms and triangles
10-1 reas of Parallelograms and Triangles ommon ore State Standards G-MG..1 Use geometric sapes, teir measures, and teir properties to descrie ojects. G-GPE..7 Use coordinates to compute perimeters of
More information6 Computing Derivatives the Quick and Easy Way
Jay Daigle Occiental College Mat 4: Calculus Experience 6 Computing Derivatives te Quick an Easy Way In te previous section we talke about wat te erivative is, an we compute several examples, an ten we
More informationOptimal In-Network Packet Aggregation Policy for Maximum Information Freshness
1 Optimal In-etwork Packet Aggregation Policy for Maimum Information Fresness Alper Sinan Akyurek, Tajana Simunic Rosing Electrical and Computer Engineering, University of California, San Diego aakyurek@ucsd.edu,
More informationz = x 2 xy + y 2 clf // c6.1(2)contour Change to make a contour plot of z=xy.
190 Lecture 6 3D equations formatting Open Lecture 6. See Capter 3, 10 of text for details. Draw a contour grap and a 3D grap of z = 1 x 2 y 2 = an upper emispere. For Classwork 1 and 2, you will grap
More informationIntra- and Inter-Session Network Coding in Wireless Networks
Intra- and Inter-Session Network Coding in Wireless Networks Hulya Seferoglu, Member, IEEE, Atina Markopoulou, Member, IEEE, K K Ramakrisnan, Fellow, IEEE arxiv:857v [csni] 3 Feb Abstract In tis paper,
More informationCoarticulation: An Approach for Generating Concurrent Plans in Markov Decision Processes
Coarticulation: An Approac for Generating Concurrent Plans in Markov Decision Processes Kasayar Roanimanes kas@cs.umass.edu Sridar Maadevan maadeva@cs.umass.edu Department of Computer Science, University
More informationNotes: Dimensional Analysis / Conversions
Wat is a unit system? A unit system is a metod of taking a measurement. Simple as tat. We ave units for distance, time, temperature, pressure, energy, mass, and many more. Wy is it important to ave a standard?
More informationAnnouncements SORTING. Prelim 1. Announcements. A3 Comments 9/26/17. This semester s event is on Saturday, November 4 Apply to be a teacher!
Announcements 2 "Organizing is wat you do efore you do someting, so tat wen you do it, it is not all mixed up." ~ A. A. Milne SORTING Lecture 11 CS2110 Fall 2017 is a program wit a teac anyting, learn
More informationAreas of Triangles and Parallelograms. Bases of a parallelogram. Height of a parallelogram THEOREM 11.3: AREA OF A TRIANGLE. a and its corresponding.
11.1 Areas of Triangles and Parallelograms Goal p Find areas of triangles and parallelograms. Your Notes VOCABULARY Bases of a parallelogram Heigt of a parallelogram POSTULATE 4: AREA OF A SQUARE POSTULATE
More informationMinimizing Memory Access By Improving Register Usage Through High-level Transformations
Minimizing Memory Access By Improving Register Usage Troug Hig-level Transformations San Li Scool of Computer Engineering anyang Tecnological University anyang Avenue, SIGAPORE 639798 Email: p144102711@ntu.edu.sg
More informationUtilizing Call Admission Control to Derive Optimal Pricing of Multiple Service Classes in Wireless Cellular Networks
Utilizing Call Admission Control to Derive Optimal Pricing of Multiple Service Classes in Wireless Cellular Networks Okan Yilmaz and Ing-Ray Cen Computer Science Department Virginia Tec {oyilmaz, ircen}@vt.edu
More informationRedundancy Awareness in SQL Queries
Redundancy Awareness in QL Queries Bin ao and Antonio Badia omputer Engineering and omputer cience Department University of Louisville bin.cao,abadia @louisville.edu Abstract In tis paper, we study QL
More information15-122: Principles of Imperative Computation, Summer 2011 Assignment 6: Trees and Secret Codes
15-122: Principles of Imperative Computation, Summer 2011 Assignment 6: Trees and Secret Codes William Lovas (wlovas@cs) Karl Naden Out: Tuesday, Friday, June 10, 2011 Due: Monday, June 13, 2011 (Written
More informationMean Waiting Time Analysis in Finite Storage Queues for Wireless Cellular Networks
Mean Waiting Time Analysis in Finite Storage ueues for Wireless ellular Networks J. YLARINOS, S. LOUVROS, K. IOANNOU, A. IOANNOU 3 A.GARMIS 2 and S.KOTSOOULOS Wireless Telecommunication Laboratory, Department
More informationGeorge Xylomenos and George C. Polyzos. with existing protocols and their eæciency in terms of
IP MULTICASTING FOR WIRELESS MOBILE OSTS George Xylomenos and George C. Polyzos fxgeorge,polyzosg@cs.ucsd.edu Computer Systems Laboratory Department of Computer Science and Engineering University of California,
More informationImplementation of Integral based Digital Curvature Estimators in DGtal
Implementation of Integral based Digital Curvature Estimators in DGtal David Coeurjolly 1, Jacques-Olivier Lacaud 2, Jérémy Levallois 1,2 1 Université de Lyon, CNRS INSA-Lyon, LIRIS, UMR5205, F-69621,
More information2D transformations Homogeneous coordinates. Uses of Transformations
2D transformations omogeneous coordinates Uses of Transformations Modeling: position and resize parts of a complex model; Viewing: define and position te virtual camera Animation: define ow objects move/cange
More informationInterference and Diffraction of Light
Interference and Diffraction of Ligt References: [1] A.P. Frenc: Vibrations and Waves, Norton Publ. 1971, Capter 8, p. 280-297 [2] PASCO Interference and Diffraction EX-9918 guide (written by Ann Hanks)
More informationAn Explicit Formula for Generalized Arithmetic-Geometric Sum
Applied Matematical Sciences, Vol. 9, 2015, no. 114, 5687-5696 HIKARI Ltd, www.m-ikari.com ttp://dx.doi.org/10.12988/ams.2015.57481 An Explicit Formula for Generalized Aritmetic-Geometric Sum Roberto B.
More information12.2 Investigate Surface Area
Investigating g Geometry ACTIVITY Use before Lesson 12.2 12.2 Investigate Surface Area MATERIALS grap paper scissors tape Q U E S T I O N How can you find te surface area of a polyedron? A net is a pattern
More informationHW 2 Bench Table. Hash Indexes: Chap. 11. Table Bench is in tablespace setq. Loading table bench. Then a bulk load
Has Indexes: Cap. CS64 Lecture 6 HW Benc Table Table of M rows, Columns of different cardinalities CREATE TABLE BENCH ( KSEQ integer primary key, K5K integer not null, K5K integer not null, KK integer
More information( )( ) ( ) MTH 95 Practice Test 1 Key = 1+ x = f x. g. ( ) ( ) The only zero of f is 7 2. The only solution to g( x ) = 4 is 2.
Mr. Simonds MTH 95 Class MTH 95 Practice Test 1 Key 1. a. g ( ) ( ) + 4( ) 4 1 c. f ( x) 7 7 7 x 14 e. + 7 + + 4 f g 1+ g. f 4 + 4 7 + 1+ i. g ( 4) ( 4) + 4( 4) k. g( x) x 16 + 16 0 x 4 + 4 4 0 x 4x+ 4
More informationA Novel QC-LDPC Code with Flexible Construction and Low Error Floor
A Novel QC-LDPC Code wit Flexile Construction and Low Error Floor Hanxin WANG,2, Saoping CHEN,2,CuitaoZHU,2 and Kaiyou SU Department of Electronics and Information Engineering, Sout-Central University
More informationOn the Use of Radio Resource Tests in Wireless ad hoc Networks
Tecnical Report RT/29/2009 On te Use of Radio Resource Tests in Wireless ad oc Networks Diogo Mónica diogo.monica@gsd.inesc-id.pt João Leitão jleitao@gsd.inesc-id.pt Luis Rodrigues ler@ist.utl.pt Carlos
More informationSymmetric Tree Replication Protocol for Efficient Distributed Storage System*
ymmetric Tree Replication Protocol for Efficient Distributed torage ystem* ung Cune Coi 1, Hee Yong Youn 1, and Joong up Coi 2 1 cool of Information and Communications Engineering ungkyunkwan University
More informationProceedings of the 8th WSEAS International Conference on Neural Networks, Vancouver, British Columbia, Canada, June 19-21,
Proceedings of te 8t WSEAS International Conference on Neural Networks, Vancouver, Britis Columbia, Canada, June 9-2, 2007 3 Neural Network Structures wit Constant Weigts to Implement Dis-Jointly Removed
More informationover The idea is to construct an algorithm to solve the IVP ODE (8.1)
Runge- Ku(a Methods Review of Heun s Method (Deriva:on from Integra:on) The idea is to construct an algorithm to solve the IVP ODE (8.1) over To obtain the solution point we can use the fundamental theorem
More informationYou should be able to visually approximate the slope of a graph. The slope m of the graph of f at the point x, f x is given by
Section. Te Tangent Line Problem 89 87. r 5 sin, e, 88. r sin sin Parabola 9 9 Hperbola e 9 9 9 89. 7,,,, 5 7 8 5 ortogonal 9. 5, 5,, 5, 5. Not multiples of eac oter; neiter parallel nor ortogonal 9.,,,
More informationHASH ALGORITHMS: A DESIGN FOR PARALLEL CALCULATIONS
HASH ALGORITHMS: A DESIGN FOR PARALLEL CALCULATIONS N.G.Bardis Researc Associate Hellenic Ministry of te Interior, Public Administration and Decentralization 8, Dragatsaniou str., Klatmonos S. 0559, Greece
More informationFinite difference stencils
> retart; wit(pdetool): wit(plot): wit(linearalgebra): Finite difference olution of te advection equation Te purpoe of ti workeet i to dicu te variou finite difference tencil one can ue to olve te diffuion
More informationGrid Adaptation for Functional Outputs: Application to Two-Dimensional Inviscid Flows
Journal of Computational Pysics 176, 40 69 (2002) doi:10.1006/jcp.2001.6967, available online at ttp://www.idealibrary.com on Grid Adaptation for Functional Outputs: Application to Two-Dimensional Inviscid
More informationTwo Modifications of Weight Calculation of the Non-Local Means Denoising Method
Engineering, 2013, 5, 522-526 ttp://dx.doi.org/10.4236/eng.2013.510b107 Publised Online October 2013 (ttp://www.scirp.org/journal/eng) Two Modifications of Weigt Calculation of te Non-Local Means Denoising
More informationOvercomplete Steerable Pyramid Filters and Rotation Invariance
vercomplete Steerable Pyramid Filters and Rotation Invariance H. Greenspan, S. Belongie R. Goodman and P. Perona S. Raksit and C. H. Anderson Department of Electrical Engineering Department of Anatomy
More informationover The idea is to construct an algorithm to solve the IVP ODE (9.1)
Runge- Ku(a Methods Review of Heun s Method (Deriva:on from Integra:on) The idea is to construct an algorithm to solve the IVP ODE (9.1) over To obtain the solution point we can use the fundamental theorem
More information16th European Signal Processing Conference (EUSIPCO 2008), Lausanne, Switzerland, August 25-29, 2008, copyright by EURASIP
16t European Signal Processing Conference (EUSIPCO 008), Lausanne, Switzerland, August 5-9, 008, copyrigt by EURASIP ADAPTIVE WINDOW FOR LOCAL POLYNOMIAL REGRESSION FROM NOISY NONUNIFORM SAMPLES A. Sreenivasa
More information