Kevin Mitchell Assignment #2 Solutions MACM 316
|
|
- Evelyn Hood
- 6 years ago
- Views:
Transcription
1 Machine arithmatic. BF, p29 #25 The binomial coefficient ( m m! = k k!(m k! describes the number of ways of choosing a subset of k objects from a set of m elements... ( point Suppose decimal machine numbers are of the form ± 0.d d 2 d 3 d 4 0 n d 9 0 d i 9 i = 2, 3, 4 n 5 ( m What is the largest value of m for which the binomial coefficient can be computed for all k directly from the k above definition without overflowing? Since the factorial (! is undefined on negative integers, we have the following restriction on k: It is therefore also true that and and therefore 0 k m k! m! (m k! m! k! (m k! k!(m k! so that ( m m! = k k!(m k! m! ( m and therefore m! is the largest number stored during the computation of regardless of the value of k. k From the specification above, the largest representable number is we therefore must find an integer m such that m! (m +! > to find this value we compute, (2, (2(3,... and check for overflow at each step. This is carried out by the very simple script fact max.m:
2 m=; mfact=m; o f l = ˆ5; while mfact<o f l m=m+; mfact=mfact m; When the while loop is completed m = 8, which is the first number to overflow. Indeed it can be checked that 7! = ! = and therefore m = 7 is the largest value that can be compute without overflowing...2 ( point It can be shown that there is a safer method of computing therefore ( m k. Note that m! (m(m... (m k + (m k(m k... ( = (m k! (m k(m k... ( = (m(m... (m k + m! (m(m... (m k + = k!(m k! (k(k... ( and furthermore the above expression has k factors in both the numerator and the denominator. They can therefore be paired in a one-to-one manner. It is clear that in order to minimise the size of the maximum number stored, it is best to pair the largest numerator with the largest denominator (i.e. m with k, m with k, etc. We are therefore dealing with a sequence of factors like..3 ( point m! k!(m k! = ( ( m k We now use the improved formula to compute ( m ( m = 3 3 ( m... k ( m 2 ( m k + ( m 2 Note that since we have already established k m (and therefore (k (m, etc. each of the factors described above is larger than one. Therefore each additional multiplication increases the result so that the largest number stored in the computation is the final result itself. We therefore seek m such that We know that therefore we bisect the function ( m ( m + 3 f(m = ( ( m on the initial interval m [, 0 6 ] with binom bisect.m > = 6 m ( m 2 3m
3 a=;b=e6 ;%guess f o r i n t e r v a l % m must always be an i n t e g e r. We round down to be on the s a f e s i d e m=floor ( ( a+b / 2 ; fun=i n l i n e ( m (mˆ2 3 m+2/ ˆ5 ; i t s =0; while b a> switch sign ( fun (m sign ( fun ( b case b=m; case a=m; case 0 break ; m=floor ( ( a+b / 2 ; i t s=i t s +; which s m = 8707 as the maximum allowable m after 20 iterations...4 (2 points While we have hitherto been investigating problems with overflow, we now investigate the effects of finite precision. In particular ( we use 4 digit chopping arithmetic to compute the number of possible 5-card hands in a 52-card deck 52 (i.e.. To implement chopping arithmatic in the least tedious manner possible, the following xchop=chop(x,t 5 function is provided function x=chop ( x, t % xchop=chop ( x, t % Loosely based on chop ( pg. 30 o f % RL Burden, JD Faires. Numerical Analysis 7 th ed. Thompson % chops the argument x ( which may be a v e c t o r o f v a l u e s to t % s i g n i f i c a n t decimal p l a c e s. idx=find ( x ; %only operate on nonzero v a l u e s o f x e=floor ( log0 ( abs ( x ( idx + ;% compute the exponent % put a l l the d i g i t s we w i l l keep to the l e f t o f the decimal p l a c e x ( idx=x ( idx. 0. ˆ ( t e ; %t r u n c a t e e v e r y t h i n g to the r i g h t o f the decimal p l a c e x ( idx= f i x ( x ( idx ; % put the exponent back how we found i t x ( idx=x ( idx. 0. ˆ ( e t ; To properly simulate 4 digit ( chopping arithmetic, we must apply this chop function to every intermediate result 52 stored in the computation of. This is implemented in the function bchop=binomchop(m,k,t 5 function [ bchop]=binomchop (m, k, t % [ r e t ]=binomchop (m, k, t % compute the binomial c o e f f i c i e n t % / m \ % \ k / % using a method c a r e f u l to avoid o v e r f l o w % and 4 d i g i t chopping a r i t h m e t i c c o u r t e s y o f the chop ( f u n c t i o n bchop =.0; for i i =0:k
4 bchop=chop ( bchop chop ( chop (m i i, t / chop ( k i i, t, t, t ; In order to assess the accuracy of the computation performed using chopping arithmetic, we compute a much more accurate result by providing binomchop( with a t argument greater than 6. In this case, our result will be dominated by limitations on machine precision, rather than the artificial limitation provided by chop(. octave :> format long octave :2> binomchop ( 5 2, 5, 4 ans = octave :3> binomchop ( 5 2, 5, 2 0 ans = octave :4> a b s e r r=abs ( binomchop (52,5,20 binomchop ( 5 2, 5, 4 a b s e r r = octave :4> r e l e r r=abs ( binomchop (52,5,20 binomchop ( 5 2, 5, 4 / binomchop ( 5 2, 5, 2 0 r e l e r r = e 04 We should not be too surprised by the fact that the relative error is the same order of magnitude as the smallest significant digit stored. 2 Rootfinding Write a function that solves f(x = 0 x (a, b in a safe way using Newton s method when stable, and otherwise resorting to bisection; a non-zero error code or print an error message if the method fails. This scheme is implemented in [soln,rc]=newtonsafe(fun,a,b,tol
5 function [ p, rc ] = newtonsafe ( f dx, a, b, t o l % [ soln, rc ] = newtonsafe ( fun, a, b, t o l % Use Newton s method to f i n d a root f ( x=0 on the i n t e r v a l xe(a, b % f a l l i n g back to b i s e c t i o n i f the s o l u t i o n g e t s taken o u t s i d e the % s p e c i f i e d i n t e r v a l. % % Loosely based on r t s a f e ( on pg.366 o f % WH Press, SA Teukolsky, WT V e t t e r l i n g, BP Flannery. Numerical Recipes % in C. 2nd Ed. Cambridge U n i v e r s i t y Press 988. % % My most n o t a b l e omission from above i s t h a t t h i s code only f a l l s back % to b i s e c t i o n i f Newton f a i l s, but not i f i t i s merely converging % s l o w l y. I found t h i s c h o i c e l e d to fewer t o t a l i t e r a t i o n s f o r the % quantum problem c o n s i d e r e d. % % INPUTS % f d x : f u n c t i o n p o i n t e r o f the form [ f, dx]= f d x ( x. % f i s the v a l u e o f the f u n c t i o n whose root we are l o o k i n g f o r % e v a l u a t e d at x. dx=f ( x / f ( x i s the newton c o r r e c t i o n at x. % Returning dx i n s t e a d o f f ( x d i r e c t l y a v o i d s the problem where % f ( x i s zero at the r o ot ( though you should p r o b a b l y r e t h i n k % your problem i f t h i s i s r e a l l y the case. % a, b : l e f t and r i g h t e x c l u s i v e p o i n t s f o r the i n t e r v a l in which to search % t o l : t o l e r a n c e d e s i r e d f o r the s o l u t i o n % % OUTPUTS % s o l n : computed v a l u e f o r root o f f ( s o l n =0 solne (a, b % rc : code i n d i c a t i n g s u c c e s s or mode o f f a i l u r e % 0 : s u c c e s s % : got NaN from t he f u n c t i o n e v a l u a t i o n ( can t b i s e c t t h a t 2 : % f u n c t i o n e v a l u a t i o n has non n e g l i g i b l e imaginary part ( can t % b i s e c t t h a t % 3 : the f u n c t i o n has the same s i g n on e i t h e r s i d e % o f the i n t e r v a l ( can t b i s e c t t h a t % 4 : v a l u e o f r o ot i s no % l o n g e r changing, but t o l e r a n c e was not met % SUCCESS =0; NOTFINITE=; COMPLEX =2; NOBRACKET=3; STALLED =4; %the i n t e r v a l i s e x c l u s i v e, but b i s e c t i o n i s i n c l u s i v e, so add or s u b t r a c t %the s m a l l e s t number p o s s i b l e from the points a=a+eps ( a ; b=b eps ( b ; r c=success;% optimism! p=(a+b / 2 ; % i n i t i a l guess i s the midpoint o f the provided i n t e r v a l % t h i s i s the only p l a c e we e x p l i c i t l y e v a l u a t e the f u n c t i o n at the % points % WARNING: t h i s i s a l s o the only p l a c e t h a t f a=f ( a. % As the a l g o r i t h m proceeds to r e f i n e the i n t e r v a l, and in p a r t i c u l a r, % change the p o s i t i o n o f a, then f a =f ( a s i n c e f a w i l l s i l l be f i x e d % the e v a l u a t i o n immediately below. Keep t h i s in mind i f f o r some
6 % reason you d e c i d e you want to l o o k at f a during i t e r a t i o n. f a=f d x ( a ; fb=f d x ( b ; i f f a==0 disp ( l e f t point i s the exact s o l u t i o n! p=a ; e l s e i f fb==0 disp ( r i g h t point i s the exact s o l u t i o n! p=b ; % Since t h e r e are d i s t i n c t f l o a t i n g p o i n t numbers I n f/ I n f % with o p p o s i t e signs, t he r oot i s s t i l l b i s e c t a b l e i f f ( a or f ( b i s i n f i n i t e. % Therefore, we only check f o r e v i l NaN which does not have a s i g n. i f any( isnan ( [ fa, fb ] disp ( got NaN e v a l u a t i n g f u n c t i o n at points! r c=notfinite; % Make sure t h a t f ( a and f ( b have zero imaginary p a r t s to w i t h i n % machine e p s i l o n o f the r e a l p a rt. e l s e i f or ( i s c o m p l e x l o o s e ( f a, i s c o m p l e x l o o s e ( f a disp ( f u n c t i o n i s complex at points! r c=complex; % In case t h e r e are imaginary p a r t s t h a t we have j u s t decided are n e g l i g i b l e % by a l l o w i n g e x e c u t i o n to proceed, we s e t them to zero now. f a=real ( f a ; fb=real ( fb ; % check t h a t the r o ot i s indeed b r a c k e t e d i f sign ( f a sign ( fb== disp ( root i s not bracketed by points! r c=nobracket; done=f a l s e ; newt=0; b i s c =0; % loop u n t i l we have marked done=t r u e ; while done [ pn, fp, absdx ]=newton ( f dx, p ; % same checks as f o r the i n i t i a l point e v a l u a t i o n s i f isnan ( fp disp ( got NaN f o r f ( p r c=notfinite; break ; e l s e i f i s c o m p l e x l o o s e ( fp disp ( got complex value f o r f ( p r c=complex; break ; % d i s c a r d n e g l i g i b l e imaginary components as b e f o r e f o r the points fp=real ( fp ; [ pb, a, b, fb ]= b i s e c t (p, a, b, fp, fb ; % i f we meet the t o l e r a n c e, t h i s i s the l a s t i t e r a t i o n
7 % but we might as w e l l go ahead and update our computed % s o l u t i o n with t he work we ve a l r e a d y done so we don t break % the loop here i f errorbound ( absdx, a, b< t o l done=true ; i f and ( i s s e n s i b l e r o o t ( pn, i s i n i n t e r v a l ( a, pn, b % Newton s o l u t i o n i s good i f and ( p==pn, done disp ( s o l u t i o n no l o n g e r changing, but t o l e r a n c e was not met r c=stalled; break ; p=real ( pn ;% d i s c a r d n e g l i g i b l e imaginary components newt=newt+;%increment number o f Newton i t e r a t i o n s else % Newton s o l u t i o n isn t good, use b i s e c t i o n s o l u t i o n i f and ( p==pb, done disp ( s o l u t i o n no l o n g e r changing, but t o l e r a n c e was not met r c=stalled; break ; p=pb ; b i s c=b i s c +;%increment number o f b i s e c t i o n i t e r a t i o n s % show s t a t i s t i c s disp ( [ newt b i s c t o t ] disp ( [ newt, bisc, newt+b i s c ] function e r r=errorbound ( absdx, a, b % Both our e s t i m a t e and t he e x a c t s o l u t i o n must be in the i n t e r v a l % [ a, b ]. This g i v e s ( b a as an upper bound on the e r r o r. abs ( dx w i l l % almost always be smaller, but j u s t to be c a r e f u l, we l l check both. e r r=min( absdx, ( b a ; function bool=i s c o m p l e x l o o s e ( x % check i f v a l u e has a s i g n i f i c a n t imaginary component bool=eps ( real ( x<abs (imag( x ; function bool=i s i n i n t e r v a l ( a, pn, b % check t h a t a<=pn<=b bool=0<=(pn a ( b pn ; function bool=i s s e n s i b l e r o o t ( pn % i s t h i s a s e n s i b l e number to c a l l a root? % I t should not be I n f or NaN and i t should have n e g l i g i b l e imaginary part. bool=and ( i s f i n i t e ( pn, i s c o m p l e x l o o s e ( pn ; function [ p, fp, dx]=newton ( f dx, p % Do a newton e v a l u a t i o n. This i s the only p l a c e where % we need to a c t u a l l y e v a l u a t e the provided f u n c t i o n
8 % i n s i d e the loop [ fp, dx]= f d x ( p ; p=p dx ; dx=abs ( dx ; %we only care about the s i z e, not the s i g n o f dx from here on function [ p, a, b, fb ] = b i s e c t (p, a, b, fp, fb % perform b i s e c t i o n with provided f u n c t i o n e v a l u a t i o n s at the % e s t i m a t e f o r the r o o t and the l e f t point only % We assume t h a t i t has a l r e a d y been checked t h a t % s i g n ( f ( a s i g n ( f ( b<0 % so we only need to check the s i g n o f % s i g n ( f ( p s i g n ( f ( b switch sign ( fp sign ( fb case a=p ; % we don t e v a l u a t e f ( a because we never a c t u a l l y use i t case b=p ; % don t r e e v a l u a t e f ( b, s i n c e we a l r e a d y have f ( p and now b=p fb=fp ; case 0 % We ve found the e x a c t r o ot. Make a zero s i z e d i n t e r v a l around the s o l u t i o n i f fp==0 a=p ; b=p ; fb=fp ; e l s e i f fb==0 % t h i s should p r o b a b l y never happen s i n c e we would have n o t i c e d a l r e a d y a=b ; else disp ( t h i s should never happen! p=(a+b / 2 ; 2. testing: f(x = x We confirm that newtonsafe( behaves as expected with the test function f(x = x. In this case f (x = 2 x /2 and since newtonsafe( requires a function that evaluate dx dx = f(x f (x x = 2 x /2 = 2 ( x x This demonstrates the advantage of having the function provided to newtonsafe dx instead of simply f (x. At x = 0 f f(x (x is undefined, however, the above expression instead automatically evaluates dx = lim x 0 f (x. The function we supply to newtonsafe is [f,dx]=sqrttest(x
9 function [ f, dx]= s q r t t e s t ( x % f u n c t i o n e v a l u a t i o n f ( x f=sqrt ( x ; i f nargout> % e v a l u a t i o n o f dx=f ( x / f ( x dx=2 (x sqrt ( x ; We can then find the root using octave :> [ p, rc ]= newtonsafe 0,. 3, eps newt b i s c t o t p = rc = 0 octave :2> p ans = 0 In this case we find the exact root after 6 Newton iterations. Note that even though Octave is warning about a division by zero, the 2.2 testing: f(x = log(x We next make sure that newtonsafe( finds the root of log(x by providing it with the function [f,dx]=logtest(x function [ f, dx]= l o g t e s t ( x % f u n c t i o n e v a l u a t i o n f ( x f=log ( x ; i f nargout> % e v a l u a t i o n o f dx=f ( x / f ( x dx=f. / x ; octave :3> [ p, rc ]= newtonsafe 0,. 3, eps newt b i s c t o t 5 6 p = rc = 0 octave :4> p ans = 0 This time we also perform six iterations to find the exact root, however one of them is a bisection rather than a Newtons step. 2.3 Quantum finite square well Now we seek the solution to the more challenging problem (z0 2 tan z = ( z
10 0 8 y=tan(z y=sqrt((z0/z 2 - roots 6 y z / π Figure : Finite square well roots at the intersections of the left and right-hand-side of Eq. (. where z 0 = a 2mV0 m = kg a = 30nm = Js V 0 = J 2.3. (2 points We approach the problem by plotting the left and right hand sides of Eq. ( as shown in. The roots are the points at which the two functions intersect. As seen in the figure, there are 5 real roots. There are infinitely many complex roots which are not shown graphically. We observe (and probably already knew that tan(z is singular at odd multiples of π/2 approaching + from ( the left and form the right. We also note that z0 2 z approaches + from both sides of z = 0. It is also purely imaginary when z > z 0. With the calculation of z 0 above, we find that 5 2 π < z 0 < 3π and that all real ( roots occur in the interval 5 2 π < z < 5 2π with exactly one in each interval between the singularities of tan(z: 5 2 π, 3 2 π, ( 3 2 π, 2 π, ( 2 π, 2 π, ( 2 π, 3 2 π, ( 3 2 π, 5 2 π (2 points We wish to recast Eq. ( in a form for which the roots solve f(z = 0. This is easily done by subtracting the right hand side (z0 2 f(z = tan z z Of course, in order to use Newton s method, we will need to know the derivative ( (z0 f (z = sec 2 z + z2 2 /2 0 z 3 z
11 z± 8.9e-6 λ(m x 0 (m e e e e e e e e e e 09 Table : Roots of the finite square well problem Eq. ( and the associated decay constants. We define the physical constants in the script quantum constants.m a=30e 9; m= e 3; h= e 34; V0=.7 e 20; z0=a/h sqrt (2 m V0 ; which is simply sourced into the function [f,dx]=quantum(x to be passed into newtonsafe function [ f, dx]=quantum ( z % Function to pass to newtonsafe f o r root f i n d i n g % %source in the quantum c o n s t a n t s quantum constan ts % the f u n c t i o n v a l u e i t s e l f f=tan ( z sqrt ( ( z0. / z.ˆ2 ; i f nargout> % i f we are c a l l e d with two output arguments, % go ahead and compute dx=f ( x / f ( x dx=f. / ( sec ( z.ˆ2+ z0 ˆ 2. / ( z. ˆ 3. sqrt ( ( z0. / z. ˆ 2 ; We could be more careful about computing dx, but in the present case we are content with letting Newton fail spectacularly at singularities so that we will fall back to bisection. All the roots shown in Fig. and tabulated in Table. are computed by the script doquantum.m. % S c r i p t to f i n d the r o o t s in t he quantum mechanics problem on assignment 2 % we know where tan has s i n g u l a r i t i e s in the r egion where the s q r t f u n c t i o n i s r e a l s i n g u l a r i t i e s =[ 5:2:5] pi / 2 ; a=s i n g u l a r i t i e s ( : ; b=s i n g u l a r i t i e s ( 2 : ; % We can t e x p e c t to g e t s o l u t i o n s more accurate than machine e p s i l o n o f % the l a r g e s t p o s s i b l e s o l u t i o n t o l=eps (max( abs ( s i n g u l a r i t i e s ; for i i =:length ( a % f i n d the r oot in t he c u r r e n t i n t e r v a l [ p ( i i, rc ( i i ]= newtonsafe (@quantum, a ( i i, b ( i i, t o l ; %save the r e s u l t s save quantum roots We find each root by restricting newtonsafe to the appropriate interval between the singularities of tan z.
12 2.3.3 ( point The roots shown in the first column Table can be used to compute the spatial decay rates λ using the equation λ = z a tan z where λ is the decay constant for the wave function of electrons in the well φ(x = exp( λx x a This suggests a more intuitive measure of the decay rate x 0 λ yielding φ(x = exp ( xx0 x a So that when x = x 0, the wave function will have decayed to φ(x 0 = /e 0.37 and when x = 2x 0 it will have decayed to φ(2x 0 = /e Since the probability density function is equal to the square of the wave function, P (2x 0 = φ(2x 0 2 = /e , so it would seem reasonable to space the wells 2x 0 apart to give 2% overlap between wells. Both λ and x 0 are tabulated in the second and third columns of Table.. Presumably the negative values do not correspond to bound states. The largest positive x 0 is 7nm. We should therefore place the wells 4nm apart. 2.4 BF p82 #0 (5 points Suppose p is a zero of multiplicity m of the function f(x, where the mth derivative of f(x is continuous on an open interval containing p. The following fixed-point method has g (p = 0, g(x = x m f(x f (x Since f(x has a continuous mth derivative, we can rewrite it as where q(p 0 then and so that taking the derivative g (x = f(x = (x p m q(x f (x = m(x p m q(x + (x p m q (x f(x f (x = (x pq(x mq(x + (x pq (x g(x = x = x + (x p q (x m q(x which we can evaluate at x = p since q(x 0 g (p = + 0 q (p m q(p = + 0 = 0 m(x pq(x mq(x + (x pq (x (x p + ( + (x p q (x m q(x + (x p + (x p m 0 ( + 0 m q (x q(x q (p q(p This therefore ensures quadratic convergence of the method of Eq. (2. 2 ( + 2 ( + 0 m (x p m q (p q(p q (x q(x (2
13 2.5 BF p74 #33 (5 points We must find p such that 2 = + p ( p 2 p + p 2 ( = ( + p p p + p We can reformulate this in terms of a zero finding problem ( 2 p f(p = ( + p p + p 2 The derivative of this function is ( 2 ( ( 20 f p p (p = p + p 2 2 ( + p p + p 2 p + p 2 + p ( p + p 2 2 ( 2p We use newtonsafe( again with the function [f,dx]=raquet(x function [ f, dx]= raquet ( x f =(+x. ( x./( x+x. ˆ 2. ˆ 2 ; i f nargout> df =2 (+x. ( x./( x+x. ˆ 2. ˆ 2 0. (. / ( x+x.ˆ2+ x./( x+x.ˆ2. ( 2 x ; df=df+(x./( x+x. ˆ 2. ˆ 2 ; dx=f. / df ; octave :> format long octave :2> newtonsafe 0,, eps newt b i s c t o t ans = BF p72 #9 (5 points The enticing alternative method for implementing the secant method p n = f(p n p n 2 f(p n 2 p n f(p n f(p n 2 is a classic example of how numerical precision is lost by subtraction. When subtracting two finite precision numbers, the number of accurate digits in the result is the number of digits that can be stored minus the number of digits that the two terms have in common. This becomes a big problem when two very nearby numbers are subtracted. as in both the numerator and the denominator of this expression as n. The standard method p n p n 2 p n = p n f(p n f(p n f(p n 2 p n p n 2 also has such a quotient of nearby subtractions ( f(p n f(p n 2, it however has two distinct factors in the numerator that are going to zero (f(p n and p n p n 2, so that the whole problematic term ts to zero as we approach the root. It s inaccuracy therefore matters less and less the closer we get.
Math 473: Practice Problems for Test 1, Fall 2011, SOLUTIONS
Math 473: Practice Problems for Test 1, Fall 011, SOLUTIONS Show your work: 1. (a) Compute the Taylor polynomials P n (x) for f(x) = sin x and x 0 = 0. Solution: Compute f(x) = sin x, f (x) = cos x, f
More informationSolution of Algebric & Transcendental Equations
Page15 Solution of Algebric & Transcendental Equations Contents: o Introduction o Evaluation of Polynomials by Horner s Method o Methods of solving non linear equations o Bracketing Methods o Bisection
More informationNumerical Methods in Informatics
Numerical Methods in Informatics Lecture 2, 30.09.2016: Nonlinear Equations in One Variable http://www.math.uzh.ch/binf4232 Tulin Kaman Institute of Mathematics, University of Zurich E-mail: tulin.kaman@math.uzh.ch
More informationRoot Finding Convergence Analysis
Root Finding Convergence Analysis Justin Ross & Matthew Kwitowski November 5, 2012 There are many different ways to calculate the root of a function. Some methods are direct and can be done by simply solving
More informationR x n. 2 R We simplify this algebraically, obtaining 2x n x n 1 x n x n
Math 42 Homework 4. page 3, #9 This is a modification of the bisection method. Write a MATLAB function similar to bisect.m. Here, given the points P a a,f a and P b b,f b with f a f b,we compute the point
More informationMath Introduction to Numerical Methods - Winter 2011 Homework 2 Assigned: Friday, January 14, Due: Thursday, January 27,
Math 371 - Introduction to Numerical Methods - Winter 2011 Homework 2 Assigned: Friday, January 14, 2011. Due: Thursday, January 27, 2011.. Include a cover page. You do not need to hand in a problem sheet.
More informationNumerical Methods Dr. Sanjeev Kumar Department of Mathematics Indian Institute of Technology Roorkee Lecture No 7 Regula Falsi and Secant Methods
Numerical Methods Dr. Sanjeev Kumar Department of Mathematics Indian Institute of Technology Roorkee Lecture No 7 Regula Falsi and Secant Methods So welcome to the next lecture of the 2 nd unit of this
More informationMidterm Review. Igor Yanovsky (Math 151A TA)
Midterm Review Igor Yanovsky (Math 5A TA) Root-Finding Methods Rootfinding methods are designed to find a zero of a function f, that is, to find a value of x such that f(x) =0 Bisection Method To apply
More informationMath Numerical Analysis Mid-Term Test Solutions
Math 400 - Numerical Analysis Mid-Term Test Solutions. Short Answers (a) A sufficient and necessary condition for the bisection method to find a root of f(x) on the interval [a,b] is f(a)f(b) < 0 or f(a)
More informationLIMITS AND DERIVATIVES
2 LIMITS AND DERIVATIVES LIMITS AND DERIVATIVES 2.2 The Limit of a Function In this section, we will learn: About limits in general and about numerical and graphical methods for computing them. THE LIMIT
More informationWhat Every Programmer Should Know About Floating-Point Arithmetic DRAFT. Last updated: November 3, Abstract
What Every Programmer Should Know About Floating-Point Arithmetic Last updated: November 3, 2014 Abstract The article provides simple answers to the common recurring questions of novice programmers about
More information1 ERROR ANALYSIS IN COMPUTATION
1 ERROR ANALYSIS IN COMPUTATION 1.2 Round-Off Errors & Computer Arithmetic (a) Computer Representation of Numbers Two types: integer mode (not used in MATLAB) floating-point mode x R ˆx F(β, t, l, u),
More informationFigure 1: Graph of y = x cos(x)
Chapter The Solution of Nonlinear Equations f(x) = 0 In this chapter we will study methods for find the solutions of functions of single variables, ie values of x such that f(x) = 0 For example, f(x) =
More informationOutline. Math Numerical Analysis. Intermediate Value Theorem. Lecture Notes Zeros and Roots. Joseph M. Mahaffy,
Outline Math 541 - Numerical Analysis Lecture Notes Zeros and Roots Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research
More informationMath Numerical Analysis
Math 541 - Numerical Analysis Lecture Notes Zeros and Roots Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center
More informationMATH ASSIGNMENT 03 SOLUTIONS
MATH444.0 ASSIGNMENT 03 SOLUTIONS 4.3 Newton s method can be used to compute reciprocals, without division. To compute /R, let fx) = x R so that fx) = 0 when x = /R. Write down the Newton iteration for
More informationNumerical Methods Lecture 3
Numerical Methods Lecture 3 Nonlinear Equations by Pavel Ludvík Introduction Definition (Root or zero of a function) A root (or a zero) of a function f is a solution of an equation f (x) = 0. We learn
More informationChapter 3: Root Finding. September 26, 2005
Chapter 3: Root Finding September 26, 2005 Outline 1 Root Finding 2 3.1 The Bisection Method 3 3.2 Newton s Method: Derivation and Examples 4 3.3 How To Stop Newton s Method 5 3.4 Application: Division
More informationCLASS NOTES Models, Algorithms and Data: Introduction to computing 2018
CLASS NOTES Models, Algorithms and Data: Introduction to computing 2018 Petros Koumoutsakos, Jens Honore Walther (Last update: April 16, 2018) IMPORTANT DISCLAIMERS 1. REFERENCES: Much of the material
More informationPolynomial and Rational Functions. Chapter 3
Polynomial and Rational Functions Chapter 3 Quadratic Functions and Models Section 3.1 Quadratic Functions Quadratic function: Function of the form f(x) = ax 2 + bx + c (a, b and c real numbers, a 0) -30
More informationMath 56 Homework 1 Michael Downs. ne n 10 + ne n (1)
. Problem (a) Yes. The following equation: ne n + ne n () holds for all n R but, since we re only concerned with the asymptotic behavior as n, let us only consider n >. Dividing both sides by n( + ne n
More informationNotes on floating point number, numerical computations and pitfalls
Notes on floating point number, numerical computations and pitfalls November 6, 212 1 Floating point numbers An n-digit floating point number in base β has the form x = ±(.d 1 d 2 d n ) β β e where.d 1
More informationROOT FINDING REVIEW MICHELLE FENG
ROOT FINDING REVIEW MICHELLE FENG 1.1. Bisection Method. 1. Root Finding Methods (1) Very naive approach based on the Intermediate Value Theorem (2) You need to be looking in an interval with only one
More informationOutline. Additional Nonlinear Systems. Abstract. Finding Equilibrium Points Numerically. Newton s Method
Outline Finding Equilibrium Points Numerically Additional Nonlinear Systems James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University June 13, 2017
More informationTwitter: @Owen134866 www.mathsfreeresourcelibrary.com Prior Knowledge Check 1) Simplify: a) 3x 2 5x 5 b) 5x3 y 2 15x 7 2) Factorise: a) x 2 2x 24 b) 3x 2 17x + 20 15x 2 y 3 3) Use long division to calculate:
More informationQueens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 2018 3 Lecture 3 3.1 General remarks March 4, 2018 This
More informationChapter 1 Mathematical Preliminaries and Error Analysis
Chapter 1 Mathematical Preliminaries and Error Analysis Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128A Numerical Analysis Limits and Continuity
More informationCLASS NOTES Computational Methods for Engineering Applications I Spring 2015
CLASS NOTES Computational Methods for Engineering Applications I Spring 2015 Petros Koumoutsakos Gerardo Tauriello (Last update: July 2, 2015) IMPORTANT DISCLAIMERS 1. REFERENCES: Much of the material
More informationLecture 5: Random numbers and Monte Carlo (Numerical Recipes, Chapter 7) Motivations for generating random numbers
Lecture 5: Random numbers and Monte Carlo (Numerical Recipes, Chapter 7) Motivations for generating random numbers To sample a function in a statistically controlled manner (i.e. for Monte Carlo integration)
More informationNumerical differentiation
Numerical differentiation Paul Seidel 1801 Lecture Notes Fall 011 Suppose that we have a function f(x) which is not given by a formula but as a result of some measurement or simulation (computer experiment)
More informationUnit 2: Solving Scalar Equations. Notes prepared by: Amos Ron, Yunpeng Li, Mark Cowlishaw, Steve Wright Instructor: Steve Wright
cs416: introduction to scientific computing 01/9/07 Unit : Solving Scalar Equations Notes prepared by: Amos Ron, Yunpeng Li, Mark Cowlishaw, Steve Wright Instructor: Steve Wright 1 Introduction We now
More informationter. on Can we get a still better result? Yes, by making the rectangles still smaller. As we make the rectangles smaller and smaller, the
Area and Tangent Problem Calculus is motivated by two main problems. The first is the area problem. It is a well known result that the area of a rectangle with length l and width w is given by A = wl.
More informationFunctions and their Graphs
Chapter One Due Monday, December 12 Functions and their Graphs Functions Domain and Range Composition and Inverses Calculator Input and Output Transformations Quadratics Functions A function yields a specific
More informationJim Lambers MAT 610 Summer Session Lecture 2 Notes
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 2 Notes These notes correspond to Sections 2.2-2.4 in the text. Vector Norms Given vectors x and y of length one, which are simply scalars x and y, the
More informationSOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS BISECTION METHOD
BISECTION METHOD If a function f(x) is continuous between a and b, and f(a) and f(b) are of opposite signs, then there exists at least one root between a and b. It is shown graphically as, Let f a be negative
More informationIntro to Scientific Computing: How long does it take to find a needle in a haystack?
Intro to Scientific Computing: How long does it take to find a needle in a haystack? Dr. David M. Goulet Intro Binary Sorting Suppose that you have a detector that can tell you if a needle is in a haystack,
More informationChapter 2 Polynomial and Rational Functions
Chapter 2 Polynomial and Rational Functions Overview: 2.2 Polynomial Functions of Higher Degree 2.3 Real Zeros of Polynomial Functions 2.4 Complex Numbers 2.5 The Fundamental Theorem of Algebra 2.6 Rational
More informationIntroduction to Numerical Analysis
Introduction to Numerical Analysis S. Baskar and S. Sivaji Ganesh Department of Mathematics Indian Institute of Technology Bombay Powai, Mumbai 400 076. Introduction to Numerical Analysis Lecture Notes
More informationBinary floating point
Binary floating point Notes for 2017-02-03 Why do we study conditioning of problems? One reason is that we may have input data contaminated by noise, resulting in a bad solution even if the intermediate
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 5 Nonlinear Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction
More informationCollege Algebra Through Problem Solving (2018 Edition)
City University of New York (CUNY) CUNY Academic Works Open Educational Resources Queensborough Community College Winter 1-25-2018 College Algebra Through Problem Solving (2018 Edition) Danielle Cifone
More informationPractical Numerical Analysis: Sheet 3 Solutions
Practical Numerical Analysis: Sheet 3 Solutions 1. We need to compute the roots of the function defined by f(x) = sin(x) + sin(x 2 ) on the interval [0, 3] using different numerical methods. First we consider
More informationChapter 2 Solutions of Equations of One Variable
Chapter 2 Solutions of Equations of One Variable 2.1 Bisection Method In this chapter we consider one of the most basic problems of numerical approximation, the root-finding problem. This process involves
More informationOutline. Scientific Computing: An Introductory Survey. Nonlinear Equations. Nonlinear Equations. Examples: Nonlinear Equations
Methods for Systems of Methods for Systems of Outline Scientific Computing: An Introductory Survey Chapter 5 1 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More information10.7 Trigonometric Equations and Inequalities
0.7 Trigonometric Equations and Inequalities 79 0.7 Trigonometric Equations and Inequalities In Sections 0., 0. and most recently 0., we solved some basic equations involving the trigonometric functions.
More informationLecture for Week 2 (Secs. 1.3 and ) Functions and Limits
Lecture for Week 2 (Secs. 1.3 and 2.2 2.3) Functions and Limits 1 First let s review what a function is. (See Sec. 1 of Review and Preview.) The best way to think of a function is as an imaginary machine,
More informationAlgebra 1 S1 Lesson Summaries. Lesson Goal: Mastery 70% or higher
Algebra 1 S1 Lesson Summaries For every lesson, you need to: Read through the LESSON REVIEW which is located below or on the last page of the lesson and 3-hole punch into your MATH BINDER. Read and work
More information1 + lim. n n+1. f(x) = x + 1, x 1. and we check that f is increasing, instead. Using the quotient rule, we easily find that. 1 (x + 1) 1 x (x + 1) 2 =
Chapter 5 Sequences and series 5. Sequences Definition 5. (Sequence). A sequence is a function which is defined on the set N of natural numbers. Since such a function is uniquely determined by its values
More informationChapter 13 - Inverse Functions
Chapter 13 - Inverse Functions In the second part of this book on Calculus, we shall be devoting our study to another type of function, the exponential function and its close relative the Sine function.
More informationMath 471. Numerical methods Root-finding algorithms for nonlinear equations
Math 471. Numerical methods Root-finding algorithms for nonlinear equations overlap Section.1.5 of Bradie Our goal in this chapter is to find the root(s) for f(x) = 0..1 Bisection Method Intermediate value
More information1. Method 1: bisection. The bisection methods starts from two points a 0 and b 0 such that
Chapter 4 Nonlinear equations 4.1 Root finding Consider the problem of solving any nonlinear relation g(x) = h(x) in the real variable x. We rephrase this problem as one of finding the zero (root) of a
More informationJim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More informationMAT 460: Numerical Analysis I. James V. Lambers
MAT 460: Numerical Analysis I James V. Lambers January 31, 2013 2 Contents 1 Mathematical Preliminaries and Error Analysis 7 1.1 Introduction............................ 7 1.1.1 Error Analysis......................
More informationChapter Five Notes N P U2C5
Chapter Five Notes N P UC5 Name Period Section 5.: Linear and Quadratic Functions with Modeling In every math class you have had since algebra you have worked with equations. Most of those equations have
More informationRational Exponents. Polynomial function of degree n: with leading coefficient,, with maximum number of turning points is given by (n-1)
Useful Fact Sheet Final Exam Interval, Set Builder Notation (a,b) = {x a
More information10.7 Trigonometric Equations and Inequalities
0.7 Trigonometric Equations and Inequalities 857 0.7 Trigonometric Equations and Inequalities In Sections 0., 0. and most recently 0., we solved some basic equations involving the trigonometric functions.
More information1 Functions, Graphs and Limits
1 Functions, Graphs and Limits 1.1 The Cartesian Plane In this course we will be dealing a lot with the Cartesian plane (also called the xy-plane), so this section should serve as a review of it and its
More informationPartial Fractions. (Do you see how to work it out? Substitute u = ax + b, so du = a dx.) For example, 1 dx = ln x 7 + C, x x (x 3)(x + 1) = a
Partial Fractions 7-9-005 Partial fractions is the opposite of adding fractions over a common denominator. It applies to integrals of the form P(x) dx, wherep(x) and Q(x) are polynomials. Q(x) The idea
More informationAlgebra II Vocabulary Word Wall Cards
Algebra II Vocabulary Word Wall Cards Mathematics vocabulary word wall cards provide a display of mathematics content words and associated visual cues to assist in vocabulary development. The cards should
More informationContinuity. MATH 161 Calculus I. J. Robert Buchanan. Fall Department of Mathematics
Continuity MATH 161 Calculus I J. Robert Buchanan Department of Mathematics Fall 2017 Intuitive Idea A process or an item can be described as continuous if it exists without interruption. The mathematical
More informationReference Material /Formulas for Pre-Calculus CP/ H Summer Packet
Reference Material /Formulas for Pre-Calculus CP/ H Summer Packet Week # 1 Order of Operations Step 1 Evaluate expressions inside grouping symbols. Order of Step 2 Evaluate all powers. Operations Step
More informationChapter 1: Introduction and mathematical preliminaries
Chapter 1: Introduction and mathematical preliminaries Evy Kersalé September 26, 2011 Motivation Most of the mathematical problems you have encountered so far can be solved analytically. However, in real-life,
More informationSolution of Nonlinear Equations
Solution of Nonlinear Equations (Com S 477/577 Notes) Yan-Bin Jia Sep 14, 017 One of the most frequently occurring problems in scientific work is to find the roots of equations of the form f(x) = 0. (1)
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 5. Nonlinear Equations
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T Heath Chapter 5 Nonlinear Equations Copyright c 2001 Reproduction permitted only for noncommercial, educational
More informationNotes for Chapter 1 of. Scientific Computing with Case Studies
Notes for Chapter 1 of Scientific Computing with Case Studies Dianne P. O Leary SIAM Press, 2008 Mathematical modeling Computer arithmetic Errors 1999-2008 Dianne P. O'Leary 1 Arithmetic and Error What
More informationChapter 2: Functions, Limits and Continuity
Chapter 2: Functions, Limits and Continuity Functions Limits Continuity Chapter 2: Functions, Limits and Continuity 1 Functions Functions are the major tools for describing the real world in mathematical
More informationIntroduction to Advanced Mathematics (C1) WEDNESDAY 9 JANUARY 2008
ADVANCED SUBSIDIARY GCE 475/0 MATHEMATICS (MEI) Introduction to Advanced Mathematics (C) WEDNESDAY 9 JANUARY 008 Additional materials: Answer Booklet (8 pages) MEI Examination Formulae and Tables (MF)
More informationMATH141: Calculus II Exam #4 review solutions 7/20/2017 Page 1
MATH4: Calculus II Exam #4 review solutions 7/0/07 Page. The limaçon r = + sin θ came up on Quiz. Find the area inside the loop of it. Solution. The loop is the section of the graph in between its two
More informationNumerical Solution of f(x) = 0
Numerical Solution of f(x) = 0 Gerald W. Recktenwald Department of Mechanical Engineering Portland State University gerry@pdx.edu ME 350: Finding roots of f(x) = 0 Overview Topics covered in these slides
More informationChapter 6. Nonlinear Equations. 6.1 The Problem of Nonlinear Root-finding. 6.2 Rate of Convergence
Chapter 6 Nonlinear Equations 6. The Problem of Nonlinear Root-finding In this module we consider the problem of using numerical techniques to find the roots of nonlinear equations, f () =. Initially we
More informationTu: 9/3/13 Math 471, Fall 2013, Section 001 Lecture 1
Tu: 9/3/13 Math 71, Fall 2013, Section 001 Lecture 1 1 Course intro Notes : Take attendance. Instructor introduction. Handout : Course description. Note the exam days (and don t be absent). Bookmark the
More information15 Nonlinear Equations and Zero-Finders
15 Nonlinear Equations and Zero-Finders This lecture describes several methods for the solution of nonlinear equations. In particular, we will discuss the computation of zeros of nonlinear functions f(x).
More informationArithmetic and Error. How does error arise? How does error arise? Notes for Part 1 of CMSC 460
Notes for Part 1 of CMSC 460 Dianne P. O Leary Preliminaries: Mathematical modeling Computer arithmetic Errors 1999-2006 Dianne P. O'Leary 1 Arithmetic and Error What we need to know about error: -- how
More informationZEROES OF INTEGER LINEAR RECURRENCES. 1. Introduction. 4 ( )( 2 1) n
ZEROES OF INTEGER LINEAR RECURRENCES DANIEL LITT Consider the integer linear recurrence 1. Introduction x n = x n 1 + 2x n 2 + 3x n 3 with x 0 = x 1 = x 2 = 1. For which n is x n = 0? Answer: x n is never
More informationUNC Charlotte 2004 Algebra with solutions
with solutions March 8, 2004 1. Let z denote the real number solution to of the digits of z? (A) 13 (B) 14 (C) 15 (D) 16 (E) 17 3 + x 1 = 5. What is the sum Solution: E. Square both sides twice to get
More information6.4 Logarithmic Equations and Inequalities
6.4 Logarithmic Equations and Inequalities 459 6.4 Logarithmic Equations and Inequalities In Section 6.3 we solved equations and inequalities involving exponential functions using one of two basic strategies.
More information+ 1 3 x2 2x x3 + 3x 2 + 0x x x2 2x + 3 4
Math 4030-001/Foundations of Algebra/Fall 2017 Polynomials at the Foundations: Rational Coefficients The rational numbers are our first field, meaning that all the laws of arithmetic hold, every number
More informationNumerical Methods 5633
Numerical Methods 5633 Lecture 4 Michaelmas Term 2017 Marina Krstic Marinkovic mmarina@maths.tcd.ie School of Mathematics Trinity College Dublin Marina Krstic Marinkovic 1 / 17 5633-Numerical Methods Root
More informationQuantum Field Theory Homework 3 Solution
Quantum Field Theory Homework 3 Solution 1 Complex Gaußian Integrals Consider the integral exp [ ix 2 /2 ] dx (1) First show that it is not absolutely convergent. Then we should define it as Next, show
More information(2) Generalize De Morgan s laws for n sets and prove the laws by induction. 1
ARS DIGITA UNIVERSITY MONTH 2: DISCRETE MATHEMATICS PROFESSOR SHAI SIMONSON PROBLEM SET 2 SOLUTIONS SET, FUNCTIONS, BIG-O, RATES OF GROWTH (1) Prove by formal logic: (a) The complement of the union of
More information10.7 Trigonometric Equations and Inequalities
0.7 Trigonometric Equations and Inequalities 857 0.7 Trigonometric Equations and Inequalities In Sections 0. 0. and most recently 0. we solved some basic equations involving the trigonometric functions.
More informationSOLVED PROBLEMS ON TAYLOR AND MACLAURIN SERIES
SOLVED PROBLEMS ON TAYLOR AND MACLAURIN SERIES TAYLOR AND MACLAURIN SERIES Taylor Series of a function f at x = a is ( f k )( a) ( x a) k k! It is a Power Series centered at a. Maclaurin Series of a function
More informationSection 6.3: Exponential Equations and Inequalities, from College Algebra: Corrected Edition by Carl Stitz, Ph.D. and Jeff Zeager, Ph.D.
Section 6.3: Exponential Equations and Inequalities, from College Algebra: Corrected Edition by Carl Stitz, Ph.D. and Jeff Zeager, Ph.D. is available under a Creative Commons Attribution-NonCommercial-
More informationTheme 1: Solving Nonlinear Equations
Theme 1: Solving Nonlinear Equations Prof. Mapundi K. Banda (Tuks) WTW263 Semester II 1 / 22 Sources of Errors Finite decimal representation (Rounding): Finite decimal representation will be used to represent
More informationNumerical Analysis and Computing
Numerical Analysis and Computing Lecture Notes #02 Calculus Review; Computer Artihmetic and Finite Precision; and Convergence; Joe Mahaffy, mahaffy@math.sdsu.edu Department of Mathematics Dynamical Systems
More informationNumerical Analysis: Solving Nonlinear Equations
Numerical Analysis: Solving Nonlinear Equations Mirko Navara http://cmp.felk.cvut.cz/ navara/ Center for Machine Perception, Department of Cybernetics, FEE, CTU Karlovo náměstí, building G, office 104a
More information1 Solutions to selected problems
Solutions to selected problems Section., #a,c,d. a. p x = n for i = n : 0 p x = xp x + i end b. z = x, y = x for i = : n y = y + x i z = zy end c. y = (t x ), p t = a for i = : n y = y(t x i ) p t = p
More informationAlgebra I Vocabulary Cards
Algebra I Vocabulary Cards Table of Contents Expressions and Operations Natural Numbers Whole Numbers Integers Rational Numbers Irrational Numbers Real Numbers Order of Operations Expression Variable Coefficient
More informationNumber Systems III MA1S1. Tristan McLoughlin. December 4, 2013
Number Systems III MA1S1 Tristan McLoughlin December 4, 2013 http://en.wikipedia.org/wiki/binary numeral system http://accu.org/index.php/articles/1558 http://www.binaryconvert.com http://en.wikipedia.org/wiki/ascii
More informationSolutions Key Exponents and Polynomials
CHAPTER 7 Solutions Key Exponents and Polynomials ARE YOU READY? PAGE 57. F. B. C. D 5. E 6. 7 7. 5 8. (-0 9. x 0. k 5. 9.. - -( - 5. 5. 5 6. 7. (- 6 (-(-(-(-(-(- 8 5 5 5 5 6 8. 0.06 9.,55 0. 5.6. 6 +
More informationUnit IV Derivatives 20 Hours Finish by Christmas
Unit IV Derivatives 20 Hours Finish by Christmas Calculus There two main streams of Calculus: Differentiation Integration Differentiation is used to find the rate of change of variables relative to one
More informationUnit IV Derivatives 20 Hours Finish by Christmas
Unit IV Derivatives 20 Hours Finish by Christmas Calculus There two main streams of Calculus: Differentiation Integration Differentiation is used to find the rate of change of variables relative to one
More informationSTEP Support Programme. Pure STEP 1 Questions
STEP Support Programme Pure STEP 1 Questions 2012 S1 Q4 1 Preparation Find the equation of the tangent to the curve y = x at the point where x = 4. Recall that x means the positive square root. Solve the
More information5.1. Integer Exponents and Scientific Notation. Objectives. Use the product rule for exponents. Define 0 and negative exponents.
Chapter 5 Section 5. Integer Exponents and Scientific Notation Objectives 2 4 5 6 Use the product rule for exponents. Define 0 and negative exponents. Use the quotient rule for exponents. Use the power
More informationRoot Finding For NonLinear Equations Bisection Method
Root Finding For NonLinear Equations Bisection Method P. Sam Johnson November 17, 2014 P. Sam Johnson (NITK) Root Finding For NonLinear Equations Bisection MethodNovember 17, 2014 1 / 26 Introduction The
More informationFinite and Infinite Sets
Chapter 9 Finite and Infinite Sets 9. Finite Sets Preview Activity (Equivalent Sets, Part ). Let A and B be sets and let f be a function from A to B..f W A! B/. Carefully complete each of the following
More informationNumerical Methods. Root Finding
Numerical Methods Solving Non Linear 1-Dimensional Equations Root Finding Given a real valued function f of one variable (say ), the idea is to find an such that: f() 0 1 Root Finding Eamples Find real
More informationNorth Carolina State University
North Carolina State University MA 141 Course Text Calculus I by Brenda Burns-Williams and Elizabeth Dempster August 7, 2014 Section1 Functions Introduction In this section, we will define the mathematical
More informationDecember Exam Summary
December Exam Summary 1 Lines and Distances 1.1 List of Concepts Distance between two numbers on the real number line or on the Cartesian Plane. Increments. If A = (a 1, a 2 ) and B = (b 1, b 2 ), then
More informationSolving Non-Linear Equations (Root Finding)
Solving Non-Linear Equations (Root Finding) Root finding Methods What are root finding methods? Methods for determining a solution of an equation. Essentially finding a root of a function, that is, a zero
More information