6.6 Modified Bessel Functions of Integer Order

Size: px
Start display at page:

Download "6.6 Modified Bessel Functions of Integer Order"

Transcription

1 6.6 Modified Bessel Functions of Integer Order 9 bjp=bjp*bigni bessj=bessj*bigni sum=sum*bigni if(jsum.ne.0)sum=sum+bj Accumulate the sum. jsum=1-jsum Change0to1orviceversa. if(j.eq.n)bessj=bjp Save the unnormalized answer. enddo 1 sum=.*sum-bj Compute (5.5.16) bessj=bessj/sum and use it to normalize the answer. if(.lt.0..and.mod(n,).eq.1)bessj=-bessj CITED REFERENCES AND FURTHER READING: Abramowitz, M., and Stegun, I.A. 1964, Handbook of Mathematical Functions, Applied Mathematics Series, Volume 55 (Washington: National Bureau of Standards; reprinted 1968 by Dover Publications, New York), Chapter 9. Hart, J.F., et al. 1968, Computer Approimations (New York: Wiley), 6.8, p [1] 6.6 Modified Bessel Functions of Integer Order The modified Bessel functions I n () and K n () are equivalent to the usual Bessel functions J n and Y n evaluated for purely imaginary arguments. In detail, the relationship is I n () =( i) n J n (i) K n () = π in+1 [J n (i)+iy n (i)] (6.6.1) The particular choice of prefactor and of the linear combination of J n and Y n to form K n are simply choices that make the functions real-valued for real arguments. For small arguments n, both I n () and K n () become, asymptotically, simple powers of their argument I n () 1 ( ) n n 0 n! K 0 () ln() K n () (n 1)! ( ) n n>0 (6.6.) These epressions are virtually identical to those for J n () and Y n () in this region, ecept for the factor of /π difference between Y n () and K n (). In the region

2 30 Chapter 6. Special Functions 4 3 modified Bessel functions 1 K 0 K 1 K I 0 I Figure Modified Bessel functions I 0 () through I 3 (), K 0 () through K (). n, however, the modified functions have quite different behavior than the Bessel functions, I I n () 1 π ep() K n () π π ep( ) I 3 (6.6.3) The modified functions evidently have eponential rather than sinusoidal behavior for large arguments (see Figure 6.6.1). The smoothness of the modified Bessel functions, once the eponential factor is removed, makes a simple polynomial approimation of a few terms quite suitable for the functions I 0, I 1, K 0, and K 1. The following routines, based on polynomial coefficients given by Abramowitz and Stegun [1], evaluate these four functions, and will provide the basis for upward recursion for n>1when >n. FUNCTION bessi0() REAL bessi0, Returns the modified Bessel function I 0 () for any real. REAL a DOUBLE PRECISION p1,p,p3,p4,p5,p6,p7,q1,q,q3,q4,q5,q6,q7, * q8,q9,y Accumulate polynomials in double precision. SAVE p1,p,p3,p4,p5,p6,p7,q1,q,q3,q4,q5,q6,q7,q8,q9 DATA p1,p,p3,p4,p5,p6,p7/1.0d0, d0, d0, d0, * d0, d-1, d-/ DATA q1,q,q3,q4,q5,q6,q7,q8,q9/ d0, d-1, * d-, d-, d-, d-1, * d-1, d-1, d-/

3 6.6 Modified Bessel Functions of Integer Order 31 if (abs().lt.3.75) then y=(/3.75)** bessi0=p1+y*(p+y*(p3+y*(p4+y*(p5+y*(p6+y*p7))))) a=abs() y=3.75/a bessi0=(ep(a)/sqrt(a))*(q1+y*(q+y*(q3+y*(q4 * +y*(q5+y*(q6+y*(q7+y*(q8+y*q9)))))))) FUNCTION bessk0() REAL bessk0, C USES bessi0 Returns the modified Bessel function K 0 () for positive real. REAL bessi0 DOUBLE PRECISION p1,p,p3,p4,p5,p6,p7,q1, * q,q3,q4,q5,q6,q7,y Accumulate polynomials in double precision. SAVE p1,p,p3,p4,p5,p6,p7,q1,q,q3,q4,q5,q6,q7 DATA p1,p,p3,p4,p5,p6,p7/ d0, d0, d0, * d-1,0.6698d-, d-3,0.74d-5/ DATA q1,q,q3,q4,q5,q6,q7/ d0, d-1, d-1, * d-1, d-, d-,0.5308d-3/ if (.le..0) then y=*/4.0 bessk0=(-log(/.0)*bessi0())+(p1+y*(p+y*(p3+ * y*(p4+y*(p5+y*(p6+y*p7)))))) y=(.0/) bessk0=(ep(-)/sqrt())*(q1+y*(q+y*(q3+ * y*(q4+y*(q5+y*(q6+y*q7)))))) FUNCTION bessi1() REAL bessi1, Returns the modified Bessel function I 1 () for any real. REAL a DOUBLE PRECISION p1,p,p3,p4,p5,p6,p7,q1,q,q3,q4,q5,q6,q7, * q8,q9,y Accumulate polynomials in double precision. SAVE p1,p,p3,p4,p5,p6,p7,q1,q,q3,q4,q5,q6,q7,q8,q9 DATA p1,p,p3,p4,p5,p6,p7/0.5d0, d0, d0, * d0, d-1, d-,0.3411d-3/ DATA q1,q,q3,q4,q5,q6,q7,q8,q9/ d0, d-1, * d-, d-, d-1,0.8967d-1, * d-1, d-1, d-/ if (abs().lt.3.75) then y=(/3.75)** bessi1=*(p1+y*(p+y*(p3+y*(p4+y*(p5+y*(p6+y*p7)))))) a=abs() y=3.75/a bessi1=(ep(a)/sqrt(a))*(q1+y*(q+y*(q3+y*(q4+ * y*(q5+y*(q6+y*(q7+y*(q8+y*q9)))))))) if(.lt.0.)bessi1=-bessi1

4 3 Chapter 6. Special Functions FUNCTION bessk1() REAL bessk1, C USES bessi1 Returns the modified Bessel function K 1 () for positive real. REAL bessi1 DOUBLE PRECISION p1,p,p3,p4,p5,p6,p7,q1, * q,q3,q4,q5,q6,q7,y Accumulate polynomials in double precision. SAVE p1,p,p3,p4,p5,p6,p7,q1,q,q3,q4,q5,q6,q7 DATA p1,p,p3,p4,p5,p6,p7/1.0d0, d0, d0, * d0, d-1, d-, d-4/ DATA q1,q,q3,q4,q5,q6,q7/ d0, d0, d-1, * d-1, d-, d-, d-3/ if (.le..0) then y=*/4.0 bessk1=(log(/.0)*bessi1())+(1.0/)*(p1+y*(p+ * y*(p3+y*(p4+y*(p5+y*(p6+y*p7)))))) y=.0/ bessk1=(ep(-)/sqrt())*(q1+y*(q+y*(q3+ * y*(q4+y*(q5+y*(q6+y*q7)))))) C The recurrence relation for I n () and K n () is the same as that for J n () and Y n () provided that i is substituted for. This has the effect of changing a sign in the relation, ( ) n I n+1 () = I n ()+I n 1 () ( ) (6.6.4) n K n+1 () =+ K n ()+K n 1 () These relations are always unstable for upward recurrence. For K n, itself growing, this presents no problem. For I n, however, the strategy of downward recursion is therefore required once again, and the starting point for the recursion may be chosen in the same manner as for the routine bessj. The only fundamental difference is that the normalization formula for I n () has an alternating minus sign in successive terms, which again arises from the substitution of i for in the formula used previously for J n 1=I 0 () I ()+I 4 () I 6 ()+ (6.6.5) In fact, we prefer simply to normalize with a call to bessi0. With this simple modification, the recursion routines bessj and bessy become the new routines bessi and bessk: FUNCTION bessk(n,) INTEGER n REAL bessk, USES bessk0,bessk1 Returns the modified Bessel function Kn() for positive and n. INTEGER j REAL bk,bkm,bkp,to,bessk0,bessk1 if (n.lt.) pause bad argument n in bessk to=.0/

5 6.6 Modified Bessel Functions of Integer Order 33 C bkm=bessk0() Upward recurrence for all... bk=bessk1() do 11 j=1,n-1...and here it is. bkp=bkm+j*to*bk bkm=bk bk=bkp enddo 11 bessk=bk FUNCTION bessi(n,) INTEGER n,iacc REAL bessi,,bigno,bigni PARAMETER (IACC=40,BIGNO=1.0e10,BIGNI=1.0e-10) USES bessi0 Returns the modified Bessel function In() for any real and n. INTEGER j,m REAL bi,bim,bip,to,bessi0 if (n.lt.) pause bad argument n in bessi if (.eq.0.) then bessi=0. to=.0/abs() bip=0.0 bi=1.0 bessi=0. m=*((n+int(sqrt(float(iacc*n))))) Downward recurrence from even m. do 11 j=m,1,-1 Make IACC larger to increase accuracy. bim=bip+float(j)*to*bi The downward recurrence. bip=bi bi=bim if (abs(bi).gt.bigno) then Renormalize to prevent overflows. bessi=bessi*bigni bi=bi*bigni bip=bip*bigni if (j.eq.n) bessi=bip enddo 11 bessi=bessi*bessi0()/bi Normalize with bessi0. if (.lt.0..and.mod(n,).eq.1) bessi=-bessi CITED REFERENCES AND FURTHER READING: Abramowitz, M., and Stegun, I.A. 1964, Handbook of Mathematical Functions, Applied Mathematics Series, Volume 55 (Washington: National Bureau of Standards; reprinted 1968 by Dover Publications, New York), 9.8. [1] Carrier, G.F., Krook, M. and Pearson, C.E. 1966, Functions of a Comple Variable (New York: McGraw-Hill), pp. 0ff.

6.9 Fresnel Integrals, Cosine and Sine Integrals

6.9 Fresnel Integrals, Cosine and Sine Integrals 48 Chapter 6. Special Functions pll=(x*(*ll-*pmmp-(ll+m-*pmm/(ll-m pmm=pmmp pmmp=pll enddo plgndr=pll CITED REFERENCES AND FURTHER READING: Magnus, W., and Oberhettinger, F. 949, Formulas and Theorems

More information

6.9 Fresnel Integrals, Cosine and Sine Integrals

6.9 Fresnel Integrals, Cosine and Sine Integrals 48 Chapter 6. Special Functions pll=(x*(*ll-*pmmp-(ll+m-*pmm/(ll-m pmm=pmmp pmmp=pll plgndr=pll CITED REFERENCES AND FURTHER READING: Magnus, W., and Oberhettinger, F. 949, Formulas and Theorems for the

More information

If your problem requires a series of related binomial coefficients, a good idea is to use recurrence relations, for example. = n k

If your problem requires a series of related binomial coefficients, a good idea is to use recurrence relations, for example. = n k 6. Incomplete Gamma Function 9 If your problem requires a series of related binomial coefficients, a good idea is to use recurrence relations, for example ( n +1 = n +1 ( ( ( n n n = + k ( n k +1 n k +1

More information

There are also formulas of higher order for this situation, but we will refrain from giving them.

There are also formulas of higher order for this situation, but we will refrain from giving them. 136 Chapter 4. Integration of Functions N = 1 2 3 4 (total after N = 4) Figure 4.2.1. Sequential calls to the routine trapzd incorporate the information from previous calls and evaluate the integrand only

More information

Section 4.4 Rational Functions and Their Graphs. 1, the line x = 0 (y-axis) is its vertical asymptote.

Section 4.4 Rational Functions and Their Graphs. 1, the line x = 0 (y-axis) is its vertical asymptote. Section 4.4 Rational Functions and Their Graphs p( ) A rational function can be epressed as where p() and q() are q( ) 3 polynomial functions and q() is not equal to 0. For eample, 16 is a rational function.

More information

Section 4.4 Rational Functions and Their Graphs

Section 4.4 Rational Functions and Their Graphs Section 4.4 Rational Functions and Their Graphs p( ) A rational function can be epressed as where p() and q() are q( ) 3 polynomial functions and q() is not equal to 0. For eample, is a 16 rational function.

More information

16.1 Runge-Kutta Method

16.1 Runge-Kutta Method 70 Chapter 6. Integration of Ordinary Differential Equations CITED REFERENCES AND FURTHER READING: Gear, C.W. 97, Numerical Initial Value Problems in Ordinary Differential Equations (Englewood Cliffs,

More information

A Catalog of Essential Functions

A Catalog of Essential Functions Section. A Catalog of Essential Functions Kiryl Tsishchanka A Catalog of Essential Functions In this course we consider 6 groups of important functions:. Linear Functions. Polynomials 3. Power functions.

More information

Numerical Aspects of Special Functions

Numerical Aspects of Special Functions Numerical Aspects of Special Functions Nico M. Temme In collaboration with Amparo Gil and Javier Segura, Santander, Spain. Nico.Temme@cwi.nl Centrum voor Wiskunde en Informatica (CWI), Amsterdam Numerics

More information

Limits and Derivatives (Review of Math 249 or 251)

Limits and Derivatives (Review of Math 249 or 251) Chapter 3 Limits and Derivatives (Review of Math 249 or 251) 3.1 Overview This is the first of two chapters reviewing material from calculus; its and derivatives are discussed in this chapter, and integrals

More information

TABLE OF CONTENTS CHAPTER 1 LIMIT AND CONTINUITY... 26

TABLE OF CONTENTS CHAPTER 1 LIMIT AND CONTINUITY... 26 TABLE OF CONTENTS CHAPTER LIMIT AND CONTINUITY... LECTURE 0- BASIC ALGEBRAIC EXPRESSIONS AND SOLVING EQUATIONS... LECTURE 0- INTRODUCTION TO FUNCTIONS... 9 LECTURE 0- EXPONENTIAL AND LOGARITHMIC FUNCTIONS...

More information

Boundary/Contour Fitted Grid Generation for Effective Visualizations in a Digital Library of Mathematical Functions

Boundary/Contour Fitted Grid Generation for Effective Visualizations in a Digital Library of Mathematical Functions Boundary/Contour Fitted Grid Generation for Effective Visualizations in a Digital Library of Mathematical Functions Bonita Saunders Qiming Wang National Institute of Standards and Technology Bureau Drive

More information

Supplemental Material for Spatially-Balanced Designs that Incorporate Legacy Sites

Supplemental Material for Spatially-Balanced Designs that Incorporate Legacy Sites Supplemental Material for Spatially-Balanced Designs that Incorporate Legacy Sites S.D. Foster, G.R. Hosack, E. Lawrence, R. Przeslawski, P. Hedge, M.J. Caley, N.S. Barrett, A. Williams, J. Li, T. Lynch,

More information

NAG Library Function Document nag_exp_integral (s13aac)

NAG Library Function Document nag_exp_integral (s13aac) s Approimations of Special Functions NAG Library Function Document nag_ep_integral () 1 Purpose nag_ep_integral () returns the value of the eponential integral E 1 ðþ. 2 Specification #include

More information

Skills Practice Skills Practice for Lesson 7.1

Skills Practice Skills Practice for Lesson 7.1 Skills Practice Skills Practice for Lesson.1 Name Date What s the Inverse of an Eponent? Logarithmic Functions as Inverses Vocabulary Write the term that best completes each statement. 1. The of a number

More information

The partitions Package

The partitions Package The partitions Package August 22, 2006 Type Package Title Additive partitions of integers Version 1.3-2 Date 2006-03-14 Author Maintainer Additive partitions of integers. Enumerates

More information

General Terms: Algorithms, exponential integral function, multiple precision

General Terms: Algorithms, exponential integral function, multiple precision Multiple-Precision Exponential Integral and Related Functions David M. Smith Loyola Marymount University This article describes a collection of Fortran-95 routines for evaluating the exponential integral

More information

Algorithm 804: Subroutines for the Computation of Mathieu Functions of Integer Orders

Algorithm 804: Subroutines for the Computation of Mathieu Functions of Integer Orders Algorithm 84: Subroutines for the Computation of Mathieu Functions of Integer Orders FAYEZ A. ALHARGAN Computer and Electronics Research Institute Computer subroutines in C for computing Mathieu functions

More information

.2 Transformations of Exponential Functions. Math

.2 Transformations of Exponential Functions. Math .2 Transformations of Eponential Functions Math 30-1 1 Vertical Translation Given the graph of f() = 2 g() = 2 + 3 Shifts the graph up if k > 0. The graph of f() moves upward 3 units. (, y) (, y + k) (0,

More information

NAG Library Function Document nag_legendre_p (s22aac)

NAG Library Function Document nag_legendre_p (s22aac) NAG Library Function Document nag_legendre_p () 1 Purpose nag_legendre_p () returns a sequence of values for either the unnormalized or normalized Legendre functions of the first kind Pn m or P n m ð x

More information

MATH 1A MIDTERM 1 (8 AM VERSION) SOLUTION. (Last edited October 18, 2013 at 5:06pm.) lim

MATH 1A MIDTERM 1 (8 AM VERSION) SOLUTION. (Last edited October 18, 2013 at 5:06pm.) lim MATH A MIDTERM (8 AM VERSION) SOLUTION (Last edited October 8, 03 at 5:06pm.) Problem. (i) State the Squeeze Theorem. (ii) Prove the Squeeze Theorem. (iii) Using a carefully justified application of the

More information

Package Bessel. February 19, 2015

Package Bessel. February 19, 2015 Version 0.5-5 Date 2013-12-10 Package Bessel February 19, 2015 Title Bessel -- Bessel Functions Computations and Approximations Author Martin Maechler Maintainer Martin Maechler

More information

Making Graphs from Tables and Graphing Horizontal and Vertical Lines - Black Level Problems

Making Graphs from Tables and Graphing Horizontal and Vertical Lines - Black Level Problems Making Graphs from Tables and Graphing Horizontal and Vertical Lines - Black Level Problems Black Level Hperbola. Give the graph and find the range and domain for. EXPONENTIAL Functions - The following

More information

PROGRAM EFFICIENCY & COMPLEXITY ANALYSIS

PROGRAM EFFICIENCY & COMPLEXITY ANALYSIS Lecture 03-04 PROGRAM EFFICIENCY & COMPLEXITY ANALYSIS By: Dr. Zahoor Jan 1 ALGORITHM DEFINITION A finite set of statements that guarantees an optimal solution in finite interval of time 2 GOOD ALGORITHMS?

More information

Section 4.2 Graphs of Exponential Functions

Section 4.2 Graphs of Exponential Functions 238 Chapter 4 Section 4.2 Graphs of Eponential Functions Like with linear functions, the graph of an eponential function is determined by the values for the parameters in the function s formula. To get

More information

Math 1330 Section : Rational Functions Definition: A rational function is a function that can be written in the form f ( x ), where

Math 1330 Section : Rational Functions Definition: A rational function is a function that can be written in the form f ( x ), where 2.3: Rational Functions P( x ) Definition: A rational function is a function that can be written in the form f ( x ), where Q( x ) and Q are polynomials, consists of all real numbers x such that You will

More information

EECS 477: Introduction to algorithms. Lecture 7

EECS 477: Introduction to algorithms. Lecture 7 EECS 477: Introduction to algorithms. Lecture 7 Prof. Igor Guskov guskov@eecs.umich.edu September 26, 2002 1 Lecture outline Recursion issues Recurrence relations Examples 2 Recursion: factorial unsigned

More information

Relations and Functions

Relations and Functions Relations and Functions. RELATION Mathematical Concepts Any pair of elements (, y) is called an ordered pair where is the first component (abscissa) and y is the second component (ordinate). Relations

More information

Generation of Finite Difference Formulas on Arbitrarily Spaced Grids

Generation of Finite Difference Formulas on Arbitrarily Spaced Grids MATHEMATICS OF COMPUTATION VOLUME, NUMBER OCTOBER 9, PAGES 99 Generation of Finite Difference Formulas on Arbitrarily Spaced Grids By Bengt Fornberg Abstract. Simple recursions are derived for calculating

More information

3 = Advanced Math 3 Fall Final Exam Review. Unit 1: If f(x) = x 2 + 3, g(x) = 3x + 1, and h(x) = x + 1, evaluate each.

3 = Advanced Math 3 Fall Final Exam Review. Unit 1: If f(x) = x 2 + 3, g(x) = 3x + 1, and h(x) = x + 1, evaluate each. Advanced Math Fall Final Eam Review Name: Unit 1: If f() +, g() + 1, and h() + 1, evaluate each. 1. f(g()). f(h()). g(- 4) 4. Given ff() + 9, represent its inverse as a (a) graph, (b) chart, and (c) function.

More information

16.2 Adaptive Stepsize Control for Runge-Kutta

16.2 Adaptive Stepsize Control for Runge-Kutta 708 Chapter 6. Integration of Ordinary Differential Equations call rk4(v,dv,nvar,x,h,v,derivs) if(x+h.eq.x)pause stepsize not significant in rkdumb x=x+h xx(k+)=x Store intermediate steps. do 2 i=,nvar

More information

Lesson 3: Solving Equations; Floating-point Computation

Lesson 3: Solving Equations; Floating-point Computation Lesson 3: Solving Equations; Floating-point Computation restart; A hard equation Last time we were looking at this equation. eq := * sin() = Pi/2; (1.1) Maple didn't know the solutions. solve(eq,,allsolutions);

More information

NAG Library Function Document nag_bessel_j0 (s17aec)

NAG Library Function Document nag_bessel_j0 (s17aec) s Approimations of Special Functions NAG Library Function Document nag_bessel_j0 () 1 Purpose nag_bessel_j0 () returns the value of the Bessel function J 0 ðþ. Specification #include #include

More information

NAG Library Function Document nag_bessel_j1 (s17afc)

NAG Library Function Document nag_bessel_j1 (s17afc) s Approimations of Special Functions NAG Library Function Document nag_bessel_j1 () 1 Purpose nag_bessel_j1 () returns the value of the Bessel function J 1 ðþ. Specification #include #include

More information

Test Bank Ver. 5.0: Data Abstraction and Problem Solving with C++: Walls and Mirrors, 5 th edition, Frank M. Carrano

Test Bank Ver. 5.0: Data Abstraction and Problem Solving with C++: Walls and Mirrors, 5 th edition, Frank M. Carrano Chapter 2 Recursion: The Mirrors Multiple Choice Questions 1. In a recursive solution, the terminates the recursive processing. a) local environment b) pivot item c) base case d) recurrence relation 2.

More information

Lecture 6: Combinatorics Steven Skiena. skiena

Lecture 6: Combinatorics Steven Skiena.  skiena Lecture 6: Combinatorics Steven Skiena Department of Computer Science State University of New York Stony Brook, NY 11794 4400 http://www.cs.sunysb.edu/ skiena Learning to Count Combinatorics problems are

More information

main() { double gm; gm = (sqrt((double)5) - 1)/2; printf (" The Golden Mean is %11.8f \n", gm); } OUTPUT The Golden Mean is 0.

main() { double gm; gm = (sqrt((double)5) - 1)/2; printf ( The Golden Mean is %11.8f \n, gm); } OUTPUT The Golden Mean is 0. PHYSICS 115/242 Homework 1, Solutions 1. (a) SOURCE CODE printf (" Hello World \n"); ----------- Hello World (b) i. SOURCE CODE double x; x = 1.9e10; printf (" %11.1e \n", x); ----------- 1.9e+10 ii. SOURCE

More information

Chapter 3. Interpolation. 3.1 Introduction

Chapter 3. Interpolation. 3.1 Introduction Chapter 3 Interpolation 3 Introduction One of the fundamental problems in Numerical Methods is the problem of interpolation, that is given a set of data points ( k, k ) for k =,, n, how do we find a function

More information

3.6-Rational Functions & Their Graphs

3.6-Rational Functions & Their Graphs .6-Rational Functions & Their Graphs What is a Rational Function? A rational function is a function that is the ratio of two polynomial functions. This definition is similar to a rational number which

More information

Making Graphs from a Table of Values and Understanding the Graphs of Horizontal and Vertical Lines Blue Level Problems

Making Graphs from a Table of Values and Understanding the Graphs of Horizontal and Vertical Lines Blue Level Problems Making Graphs from a Table of Values and Understanding the Graphs of Horizontal and Vertical Lines Blue Level Problems. Coordinate Triangle? We have a triangle ABC, and it has an area of units^. Point

More information

2.1 Quadraticsnts.notebook. September 10, 2018

2.1 Quadraticsnts.notebook. September 10, 2018 1 A quadratic function is a polynomial function of second degree. The graph of a quadratic function is called a parabola. 2 Standard Form: Intercept Form: Vertex Form: f(x) = a(x h) 2 + k vertex: (h, k)

More information

2012/05/29 DRAFT WG5/N1921

2012/05/29 DRAFT WG5/N1921 // DRAFT WG/N Special Mathematical Functions in Fortran ISO/IEC - : x Auxiliary to ISO/IEC : Programming Language Fortran NOTE This paper is intended to suggest some special functions for which standard

More information

Sample tasks from: Algebra Assessments Through the Common Core (Grades 6-12)

Sample tasks from: Algebra Assessments Through the Common Core (Grades 6-12) Sample tasks from: Algebra Assessments Through the Common Core (Grades 6-12) A resource from The Charles A Dana Center at The University of Texas at Austin 2011 About the Dana Center Assessments More than

More information

MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.

MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Eam Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Find the intervals on which the function is continuous. 2 1) y = ( + 5)2 + 10 A) (-, ) B)

More information

Using Characteristics of a Quadratic Function to Describe Its Graph. The graphs of quadratic functions can be described using key characteristics:

Using Characteristics of a Quadratic Function to Describe Its Graph. The graphs of quadratic functions can be described using key characteristics: Chapter Summar Ke Terms standard form of a quadratic function (.1) factored form of a quadratic function (.1) verte form of a quadratic function (.1) concavit of a parabola (.1) reference points (.) transformation

More information

ON THE INVERSE RADON TRANSFORM

ON THE INVERSE RADON TRANSFORM ON THE INVERSE RADON TRANSFORM F Chvála Institute of Thermomechanics AS CR, vvi Abstract A reconstruction of an image by application of the discrete inverse Radon transform realized by the function iradon

More information

The p-sized partitioning algorithm for fast computation of factorials of numbers

The p-sized partitioning algorithm for fast computation of factorials of numbers J Supercomput (2006) 38:73 82 DOI 10.1007/s11227-006-7285-5 The p-sized partitioning algorithm for fast computation of factorials of numbers Ahmet Ugur Henry Thompson C Science + Business Media, LLC 2006

More information

Algebra 2 Notes Name: Section 8.4 Rational Functions. A function is a function whose rule can be written as a of. 1 x. =. Its graph is a, f x

Algebra 2 Notes Name: Section 8.4 Rational Functions. A function is a function whose rule can be written as a of. 1 x. =. Its graph is a, f x Algebra Notes Name: Section 8. Rational Functions DAY ONE: A function is a function whose rule can be written as a of two polynomials. The parent rational function is f. Its graph is a, which has two separate

More information

NAG Library Function Document nag_complex_bessel_j (s17dec)

NAG Library Function Document nag_complex_bessel_j (s17dec) NAG Library Function Document nag_complex_bessel_j () 1 Purpose nag_complex_bessel_j () returns a sequence of values for the Bessel functions J þn ðþ z for complex z, non-negative and n ¼ 0; 1;...;N 1,

More information

Lecture 16: Introduction to Dynamic Programming Steven Skiena. Department of Computer Science State University of New York Stony Brook, NY

Lecture 16: Introduction to Dynamic Programming Steven Skiena. Department of Computer Science State University of New York Stony Brook, NY Lecture 16: Introduction to Dynamic Programming Steven Skiena Department of Computer Science State University of New York Stony Brook, NY 11794 4400 http://www.cs.sunysb.edu/ skiena Problem of the Day

More information

Remark on Algorithm 916: Computing the Faddeyeva and Voigt functions : Efficiency Improvements and Fortran Translation

Remark on Algorithm 916: Computing the Faddeyeva and Voigt functions : Efficiency Improvements and Fortran Translation Remark on Algorithm 916: Computing the Faddeyeva and Voigt functions : Efficiency Improvements and Fortran Translation MOFREH R. ZAGHLOUL, United Arab Emirates University 1 This remark describes efficiency

More information

College Algebra Review for Test 4

College Algebra Review for Test 4 College Algebra Review for Test Evaluate the epression b hand. ) a) /2 b) 27/3 c) 32/ d) 8-/3 Use positive rational eponents to rewrite the epression. 2) 2 Use radical notation to rewrite. 3) a-/2b/7 Solve

More information

CS 157: Assignment 5

CS 157: Assignment 5 Problem : Printing Neatly CS 157: Assignment 5 Douglas R. Lanman 4 April 006 In a word processor or in L A TEX, one routinely encounters the pretty printing problem. That is, how does one transform text

More information

Instructor: Virginia Davis Course: Foundations for College Math (1)

Instructor: Virginia Davis Course: Foundations for College Math (1) 5/19/01 Final Eam Review Ch 10,11-Virginia Davis Student: Date: Instructor: Virginia Davis Course: Foundations for College Math (1) Assignment: Final Eam Review Ch 10,11 1. Simplif b factoring. Assume

More information

Exact Calculation of Beta Inequalities

Exact Calculation of Beta Inequalities Exact Calculation of Beta Inequalities John Cook Department of Biostatistics, Box 447 The University of Texas, M. D. Anderson Cancer Center 1515 Holcombe Blvd., Houston, Texas 77030, USA cook@mdanderson.org

More information

Dyadic Interpolation Schemes

Dyadic Interpolation Schemes Dyadic Interpolation Schemes Mie Day, January 8 Given a set of equally-spaced samples of some function, we d lie to be able to mae educated guesses for the values which the function would tae on at points

More information

Cumulative Review Problems Packet # 1

Cumulative Review Problems Packet # 1 April 15, 009 Cumulative Review Problems Packet #1 page 1 Cumulative Review Problems Packet # 1 This set of review problems will help you prepare for the cumulative test on Friday, April 17. The test will

More information

Using Templates to Introduce Time Efficiency Analysis in an Algorithms Course

Using Templates to Introduce Time Efficiency Analysis in an Algorithms Course Using Templates to Introduce Time Efficiency Analysis in an Algorithms Course Irena Pevac Department of Computer Science Central Connecticut State University, New Britain, CT, USA Abstract: We propose

More information

NAG Library Function Document nag_complex_bessel_i (s18dec)

NAG Library Function Document nag_complex_bessel_i (s18dec) NAG Library Function Document nag_complex_bessel_i () 1 Purpose nag_complex_bessel_i () returns a sequence of values for the modified Bessel functions I þn ðþ z for complex z, non-negative and n ¼ 0; 1;...;N

More information

2 Computation with Floating-Point Numbers

2 Computation with Floating-Point Numbers 2 Computation with Floating-Point Numbers 2.1 Floating-Point Representation The notion of real numbers in mathematics is convenient for hand computations and formula manipulations. However, real numbers

More information

Domain: The domain of f is all real numbers except those values for which Q(x) =0.

Domain: The domain of f is all real numbers except those values for which Q(x) =0. Math 1330 Section.3.3: Rational Functions Definition: A rational function is a function that can be written in the form P() f(), where f and g are polynomials. Q() The domain of the rational function such

More information

Core Mathematics 3 Functions

Core Mathematics 3 Functions http://kumarmaths.weebly.com/ Core Mathematics 3 Functions Core Maths 3 Functions Page 1 Functions C3 The specifications suggest that you should be able to do the following: Understand the definition of

More information

SHORT ANSWER. Write the word or phrase that best completes each statement or answers the question.

SHORT ANSWER. Write the word or phrase that best completes each statement or answers the question. Eam Name SHORT ANSWER. Write the word or phrase that best completes each statement or answers the question. Perform the indicated operation on functions f and g. 1) Given f() = 3-4 and g() = -7 + 9, find

More information

x 16 d( x) 16 n( x) 36 d( x) zeros: x 2 36 = 0 x 2 = 36 x = ±6 Section Yes. Since 1 is a polynomial (of degree 0), P(x) =

x 16 d( x) 16 n( x) 36 d( x) zeros: x 2 36 = 0 x 2 = 36 x = ±6 Section Yes. Since 1 is a polynomial (of degree 0), P(x) = 9 CHAPTER POLYNOMIAL AND RATIONAL FUNCTIONS Section -. Yes. Since is a polynomial (of degree 0), P() P( ) is a rational function if P() is a polynomial.. A vertical asymptote is a vertical line a that

More information

correlated to the Michigan High School Mathematics Content Expectations

correlated to the Michigan High School Mathematics Content Expectations correlated to the Michigan High School Mathematics Content Expectations McDougal Littell Algebra 1 Geometry Algebra 2 2007 correlated to the STRAND 1: QUANTITATIVE LITERACY AND LOGIC (L) STANDARD L1: REASONING

More information

NAG Library Function Document nag_complex_bessel_y (s17dcc)

NAG Library Function Document nag_complex_bessel_y (s17dcc) NAG Library Function Document nag_complex_bessel_y () 1 Purpose nag_complex_bessel_y () returns a sequence of values for the Bessel functions Y þn ðþfor z complex z, non-negative and n ¼ 0; 1;...;N 1,

More information

Chapter 1. Math review. 1.1 Some sets

Chapter 1. Math review. 1.1 Some sets Chapter 1 Math review This book assumes that you understood precalculus when you took it. So you used to know how to do things like factoring polynomials, solving high school geometry problems, using trigonometric

More information

Easy way to Find Multivariate Interpolation

Easy way to Find Multivariate Interpolation International Journal of Emerging Trends in Science Technology IC Value: 7689 (Index Copernicus) Impact Factor: 4219 DOI: https://dxdoiorg/1018535/ijetst/v4i511 Easy way to Find Multivariate Interpolation

More information

Continued Fractions for Special Functions: Handbook and Software

Continued Fractions for Special Functions: Handbook and Software Continued Fractions for Special Functions: Handbook and Software Annie Cuyt, Franky Backeljauw, Stefan Becuwe, Michel Colman, Tom Docx and Joris Van Deun Abstract The revived interest in continued fractions

More information

You Try: Find the x-intercepts of f ( x) Find the roots (zeros, x-intercepts) of 2. x x. x 2. x 8 x x 2 8x 4 4x 32

You Try: Find the x-intercepts of f ( x) Find the roots (zeros, x-intercepts) of 2. x x. x 2. x 8 x x 2 8x 4 4x 32 1 Find the roots (zeros, -intercepts) of 1 3. 1a Find the -intercepts of 5 1 We are looking for the solutions to Method 1: Factoring Using Guess & heck Using an rea Model, fill in a Generic Rectangle with

More information

Chapter 2 Exercises and Answers

Chapter 2 Exercises and Answers Chapter 2 Exercises and nswers nswers are in blue. For Exercises -5, match the following numbers with their definition.. Number. Natural number C. Integer number D. Negative number E. Rational number unit

More information

F.BF.B.3: Graphing Polynomial Functions

F.BF.B.3: Graphing Polynomial Functions F.BF.B.3: Graphing Polynomial Functions 1 Given the graph of the line represented by the equation f(x) = 2x + b, if b is increased by 4 units, the graph of the new line would be shifted 4 units 1) right

More information

2 Computation with Floating-Point Numbers

2 Computation with Floating-Point Numbers 2 Computation with Floating-Point Numbers 2.1 Floating-Point Representation The notion of real numbers in mathematics is convenient for hand computations and formula manipulations. However, real numbers

More information

Exponential Functions. Christopher Thomas

Exponential Functions. Christopher Thomas Mathematics Learning Centre Eponential Functions Christopher Thomas c 1998 Universit of Sdne Mathematics Learning Centre, Universit of Sdne 1 1 Eponential Functions 1.1 The functions =2 and =2 From our

More information

Numerical Analysis Fall. Numerical Differentiation

Numerical Analysis Fall. Numerical Differentiation Numerical Analysis 5 Fall Numerical Differentiation Differentiation The mathematical definition of a derivative begins with a difference approimation: and as is allowed to approach zero, the difference

More information

Analysis of Algorithms Part I: Analyzing a pseudo-code

Analysis of Algorithms Part I: Analyzing a pseudo-code Analysis of Algorithms Part I: Analyzing a pseudo-code Introduction Pseudo-code representation of an algorithm Analyzing algorithms Measuring the running time and memory size of an algorithm Calculating

More information

A Cumulative Averaging Method for Piecewise Polynomial Approximation to Discrete Data

A Cumulative Averaging Method for Piecewise Polynomial Approximation to Discrete Data Applied Mathematical Sciences, Vol. 1, 16, no. 7, 331-343 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/1.1988/ams.16.5177 A Cumulative Averaging Method for Piecewise Polynomial Approximation to Discrete

More information

Divide-and-Conquer. Divide-and conquer is a general algorithm design paradigm:

Divide-and-Conquer. Divide-and conquer is a general algorithm design paradigm: Presentation for use with the textbook Data Structures and Algorithms in Java, 6 th edition, by M. T. Goodrich, R. Tamassia, and M. H. Goldwasser, Wiley, 2014 Merge Sort 7 2 9 4 2 4 7 9 7 2 2 7 9 4 4 9

More information

Pre-Calculus Summer Assignment

Pre-Calculus Summer Assignment Name: Pre-Calculus Summer Assignment Due Date: The beginning of class on September 8, 017. The purpose of this assignment is to have you practice the mathematical skills necessary to be successful in Pre-Calculus.

More information

Data Structures and Algorithms

Data Structures and Algorithms Data Structures and Algorithms About the course (objectives, outline, recommended reading) Problem solving Notions of Algorithmics (growth of functions, efficiency, programming model, example analysis)

More information

Module 3.5: nag fresnel intg Fresnel Integrals. Contents

Module 3.5: nag fresnel intg Fresnel Integrals. Contents Special Functions Module Contents Module 3.5: nag fresnel intg Fresnel Integrals nag fresnel intg contains procedures for approximating the Fresnel integrals S(x) and C(x). Contents Introduction... 3.5.3

More information

x 2 + 8x - 12 = 0 April 18, 2016 Aim: To review for Quadratic Function Exam #1 Homework: Study Review Materials

x 2 + 8x - 12 = 0 April 18, 2016 Aim: To review for Quadratic Function Exam #1 Homework: Study Review Materials im: To review for Quadratic Function Exam #1 Homework: Study Review Materials o Now - Solve using any strategy. If irrational, express in simplest radical form x 2 + 8x - 12 = 0 Review Topic Index 1. Transformations

More information

Algebra II Quadratic Functions

Algebra II Quadratic Functions 1 Algebra II Quadratic Functions 2014-10-14 www.njctl.org 2 Ta b le o f C o n te n t Key Terms click on the topic to go to that section Explain Characteristics of Quadratic Functions Combining Transformations

More information

lim x c x 2 x +2. Suppose that, instead of calculating all the values in the above tables, you simply . What do you find? x +2

lim x c x 2 x +2. Suppose that, instead of calculating all the values in the above tables, you simply . What do you find? x +2 MA123, Chapter 3: The idea of its (pp. 47-67, Gootman) Chapter Goals: Evaluate its. Evaluate one-sided its. Understand the concepts of continuity and differentiability and their relationship. Assignments:

More information

TeleBench 1.1. software benchmark data book.

TeleBench 1.1. software benchmark data book. TeleBench 1.1 software benchmark data book Table of Contents Autocorrelation...2 Bit Allocation...4 Convolutional Encoder...6 Fast Fourier Transform (FFT)...8 Viterbi Decoder... 11 1 TeleBench Version

More information

Finding Asymptotes KEY

Finding Asymptotes KEY Unit: 0 Lesson: 0 Discontinuities Rational functions of the form f ( are undefined at values of that make 0. Wherever a rational function is undefined, a break occurs in its graph. Each such break is called

More information

Floating-point representation

Floating-point representation Lecture 3-4: Floating-point representation and arithmetic Floating-point representation The notion of real numbers in mathematics is convenient for hand computations and formula manipulations. However,

More information

Divide and Conquer Algorithms

Divide and Conquer Algorithms Divide and Conquer Algorithms T. M. Murali February 19, 2009 Divide and Conquer Break up a problem into several parts. Solve each part recursively. Solve base cases by brute force. Efficiently combine

More information

7.7 Quasi- (that is, Sub-) Random Sequences

7.7 Quasi- (that is, Sub-) Random Sequences 77 Quasi- (that is, Sub-) Random Sequences 299 If you think for a minute, you will realize that equation (767) was useful only because the part of the integrand that we wanted to eliminate (e 5z ) was

More information

Algorithm Analysis and Design

Algorithm Analysis and Design Algorithm Analysis and Design Dr. Truong Tuan Anh Faculty of Computer Science and Engineering Ho Chi Minh City University of Technology VNU- Ho Chi Minh City 1 References [1] Cormen, T. H., Leiserson,

More information

Errors in Computation

Errors in Computation Theory of Errors Content Errors in computation Absolute Error Relative Error Roundoff Errors Truncation Errors Floating Point Numbers Normalized Floating Point Numbers Roundoff Error in Floating Point

More information

Math 1314 Lesson 13 Analyzing Other Types of Functions

Math 1314 Lesson 13 Analyzing Other Types of Functions Math 1314 Lesson 13 Analyzing Other Types of Functions Asymptotes We will need to identify any vertical or horizontal asymptotes of the graph of a function. A vertical asymptote is a vertical line x =

More information

Outline. First Quiz Results. Exercise Five Goals. Question Three. Questions One and Two. Exercise five if statements February 28, 2006

Outline. First Quiz Results. Exercise Five Goals. Question Three. Questions One and Two. Exercise five if statements February 28, 2006 Eercise five if statements February 8, 6 Laboratory V Program Control Using if Statements Larry Caretto Computer Science 6 Computing in Engineering and Science February 8, 6 Outline Review first quiz Summarize

More information

Computational Methods of Scientific Programming Fall 2007

Computational Methods of Scientific Programming Fall 2007 MIT OpenCourseWare http://ocw.mit.edu 12.010 Computational Methods of Scientific Programming Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.

More information

CS6702 GRAPH THEORY AND APPLICATIONS QUESTION BANK

CS6702 GRAPH THEORY AND APPLICATIONS QUESTION BANK CS6702 GRAPH THEORY AND APPLICATIONS 2 MARKS QUESTIONS AND ANSWERS 1 UNIT I INTRODUCTION CS6702 GRAPH THEORY AND APPLICATIONS QUESTION BANK 1. Define Graph. 2. Define Simple graph. 3. Write few problems

More information

Figure 1: The evaluation window. ab a b \a.b (\.y)((\.)(\.)) Epressions with nested abstractions such as \.(\y.(\.w)) can be abbreviated as \y.w. v al

Figure 1: The evaluation window. ab a b \a.b (\.y)((\.)(\.)) Epressions with nested abstractions such as \.(\y.(\.w)) can be abbreviated as \y.w. v al v: An Interactive -Calculus Tool Doug Zongker CSE 505, Autumn 1996 December 11, 1996 \Computers are better than humans at doing these things." { Gary Leavens, CSE 505 lecture 1 Introduction The -calculus

More information

Graph Algorithms. Chromatic Polynomials. Graph Algorithms

Graph Algorithms. Chromatic Polynomials. Graph Algorithms Graph Algorithms Chromatic Polynomials Graph Algorithms Chromatic Polynomials Definition G a simple labelled graph with n vertices and m edges. k a positive integer. P G (k) number of different ways of

More information

Examples of Fourier series

Examples of Fourier series Examples of Fourier series Preliminaries In[]:= Below is the definition of a periodic extension of a function defined on L, L. This definition takes a function as a variable. The function has to be inputted

More information

1. Given each equation below, graph it on the grid provided. Label all max. & min. points: Asymptotes: D: R: Asymptotes: D: R: Asymptotes: D: R:

1. Given each equation below, graph it on the grid provided. Label all max. & min. points: Asymptotes: D: R: Asymptotes: D: R: Asymptotes: D: R: Name: Date: Math 10/11 Enriched: Section 5. Graphing Rational Functions and Applications 1. Given each equation below, graph it on the grid provided. Label all ma. & min. points: a) 4 b) 4 c) 5 d) 9 3

More information

Exploring Hyperplane Arrangements and their Regions

Exploring Hyperplane Arrangements and their Regions Exploring Hyperplane Arrangements and their Regions Jonathan Haapala Abstract This paper studies the Linial, Coxeter, Shi, and Nil hyperplane arrangements and explores ideas of how to count their regions.

More information