Symbolic computing 1: Proofs with SymPy
|
|
- Melvyn Ball
- 5 years ago
- Views:
Transcription
1 Bachelor of Ecole Polytechnique Computational Mathematics, year 2, semester Lecturer: Lucas Gerin (send mail) Symbolic computing : Proofs with SymPy π Table of contents Introduction to SymPy Let SymPy do the proof Archimedes vs SymPy Solving equations (Bonus) Decomposition of fractions # execute this part to modify the css style from IPython.core.display import HTML def css_styling(): styles = open("./style/custom2.css").read() return HTML(styles) css_styling() ## loading python libraries # necessary to display plots inline: %matplotlib inline # load the libraries import matplotlib.pyplot as plt # 2D plotting library import numpy as np # package for scientific computing from pylab import * from math import * # package for mathematics (pi, arctan, sqrt, factorial...) import sympy as sympy # package for symbolic computation from sympy import *
2 Introduction to SymPy Remark. From now on every Lab is graded and is individual. You must upload your.pynb le each Thursday on the moodle, before :59. Using python library SymPy we can perform exact computations. For instance, compare the following scripts: print('square root of two is '+str(np.sqrt(2))+'') print('square root of two is '+str(sqrt(2))+'') square root of two is square root of two is sqrt(2) One can expand or simplify expressions: print('the square root of 40 is '+str(sqrt(40))+'') print(expand((sqrt(3)+sqrt(2))**20)) var('x') # We declare a 'symbolic' variable print(simplify((x**2-2*x + )/(x-))) the square root of 40 is 2*sqrt(0) *sqrt(6) x - Do it yourself. (A) ϕ = + 5 F = ϕ4 ϕ 2 Set. Use SymPy to simplify. +ϕ 7 phi=(+sqrt(5))/2 formula=(phi**4-phi)/(phi**7+) print("f = "+str(formula)) print("simplified F = "+str(simplify(formula))) F = (-sqrt(5)/2 - /2 + (/2 + sqrt(5)/2)**4)/( + (/2 + sqrt(5)/2)**7 ) simplified F = -4*sqrt(5)/29 + 4/29 With Sympy one can also obtain Taylor expansions of functions with series :
3 # Real variable var('x') Expression=cos(x) print('expansion of cos(x) at x=0: '+str(expression.series(x,0))) # integer variable var('n',integer=true) Expression=cos(/n) print('expansion of cos(/n) when n->+oo: '+str(expression.series(n,oo))) # oo mean Expansion of cos(x) at x=0: - x**2/2 + x**4/24 + O(x**6) Expansion of cos(/n) when n->+oo: /(24*n**4) - /(2*n**2) + + O(n** (-6), (n, oo)) SymPy can also compute with "big O's". (By default (x) is considered for x 0.) var('x') simplify((x+o(x**3))*(x+x**2+o(x**3))) x**2 + x**3 + O(x**4) Remark. A nice feature of Sympy is that you can export formulas in instance: LateX. For phi=(+sympy.sqrt(5))/2 formula=expand(phi**) print(formula) print(latex(formula)) 99/2 + 89*sqrt(5)/2 \frac{99}{2} + \frac{89 \sqrt{5}}{2} Warning: Fractions such as /4 Sympy from evaluating the expression. For example: must be introduced with Rational(,4) to keep print('(/4)^3 = '+str((/4)**3)) print('(/4)^3 = '+str(rational(,4)**3)) (/4)^3 = (/4)^3 = /64 Let SymPy do the proofs Exercise. A simple (?) recurrence
4 Do it yourself. (B) Let a, b be two real numbers, we de ne the sequence (u n ) n 0 as follows: and for u 0 = a, u = b n 2 u n + u n =. u n 2 5 u n a, b. Write a short program which returns the rst values of in terms of. The output should be something like: u_0 = a u_ = b u_2 = (b + )/a Use SymPy to make and prove a conjecture for the asymptotic behaviour of the sequence ( u n ). def InductionFormula(x,y): return (+x)/y var('a b') Sequence=[a,b] print('u_0 = a') print('u_ = b') for i in range(2,3): Sequence.append(simplify(InductionFormula(Sequence[-],Sequence[-2]))) #Sequence.append(InductionFormula(Sequence[-],Sequence[-2])) print('u_'+str(i)+' = '+str(sequence[-])) u_0 = a u_ = b u_2 = (b + )/a u_3 = (a + b + )/(a*b) u_4 = (a + )/b u_5 = a u_6 = b u_7 = (b + )/a u_8 = (a + b + )/(a*b) u_9 = (a + )/b u_0 = a u_ = b u_2 = (b + )/a Answers.. See the cell above. 2. If a, b 0, the sequence is well de ned and we observe that u 5 = u 0 and u 6 = u. Since the sequence is de ned by induction, this implies that the sequence is periodic: u n+5 = u n for every n. So if we trust Sympy the proof is done.
5 Exercise 2. A number theory theorem with Sympy Do it yourself. (C).. With Sympy simplify the following expression: 3 x x 4 8 x x x x x + + ( 3x( x 2 + 3) 2 ) ( 3x( x 2 + 3) ) ( ( x 2 + 3) 2 ) 2. Prove the following: Theorem. Every integer n can be written as the sum of three cubes of rational numbers. 3 3 var('x') A=((x**6+45*x**4-8*x**2+27)/(3*x*(x**2+3)**2))**3 B=((-x**4+30*x**2-9)/(3*x*(x**2+3)))**3 C=((-2*x**3+36*x)/((x**2+3)**2))**3 simplify(a+b+c) 8*x Answers. If we put (which are integers) then the above script shows that Therefore p = n n 4 8 n , q = 3n( n ) 2, r = n n 2 9, s = n3( n 2 + 3), t = 2 n n, u = ( n 2 + 3) 2 8n = (p/q) 3 + (r/s) 3 + (t/u) 3. n = (p/2q) 3 + (r/2s) 3 + (t/2u) 3.
6 Exercise 3. What if Archimedes had known Sympy? n n 3 2 n For, let be a regular -gon with radius. Here is : Archimedes (IIIrd century BC) used the fact that n to obtain good approximations of. 2π gets closer and closer to the unit circle Do it yourself. (D) Let L n be the length of any side of n. Compute L and use the following picture to write L n+ as a function of L n : O OC = (OB) is the bisector of ˆDOC. ˆOAC is a right angle. is the center of the circle,.
7 OB CB = BD n+ OAC A 2 = O A 2 + A C 2 = O A 2 + ( L n /2). BAC is also rectangle at A, therefore Ln+ 2 = B C 2 = A B 2 + B C 2 = ( OA ) 2 + ( L n /2). Answers. As is the bisector we have that, which both are sides of. Besides, is rectangle at. By Pythagora's theorem Putting everything together we obtain L n+ = 2 2 ( L n /2) 2. Do it yourself. (E). Write a script which computes exact expressions for the rst values L, L 2,, L n. (First try for small n 's.) 2. Find a sequence a n such that a n L n converges to π (we don't ask for a proof). Deduce some good algebraic approximations of π. Export your results in Latex in order to get nice formulas. (In order to check your formulas, you may compute numerical evaluations. With SymPy, a numerical evaluation is obtained with N(expression).) SuccessiveApproximations=[] p=8 for n in range(,p): OldValue=SuccessiveApproximations[-] NewValue=expand(sqrt(2-2*sqrt(-(OldValue**2)*Rational(,4)))) SuccessiveApproximations.append(NewValue) print(latex(simplify(3*(2**n)*newvalue))) print(n(newvalue*3*2**(n))) 6 \sqrt{- \sqrt{3} + 2} \sqrt{- \sqrt{\sqrt{3} + 2} + 2} \sqrt{- \sqrt{\sqrt{\sqrt{3} + 2} + 2} + 2} \sqrt{- \sqrt{\sqrt{\sqrt{\sqrt{3} + 2} + 2} + 2} + 2} \sqrt{- \sqrt{\sqrt{\sqrt{\sqrt{\sqrt{3} + 2} + 2} + 2} + 2} + 2} \sqrt{- \sqrt{\sqrt{\sqrt{\sqrt{\sqrt{\sqrt{3} + 2} + 2} + 2} + 2} + 2} + 2} \sqrt{- \sqrt{\sqrt{\sqrt{\sqrt{\sqrt{\sqrt{\sqrt{3} + 2} + 2} + 2} + 2} + 2} + 2} + 2}
8 Answers. When n goes large, n gets closer and closer to the unit circle. As the perimeter of is, we expect that 3 2π, 2 n L n n 3 2 n L n a n = 3 2 n n = 8 therefore we choose. For we obtain: π = Solving equations Exercise 4. Solving equations with Sympy: the easy case One can solve equations with Sympy. The following script shows how to solve x 2 = x + : var('x') # we declare variable SetOfSolutions=solve(x**2-x-,x) print(setofsolutions) [/2 + sqrt(5)/2, -sqrt(5)/2 + /2] We will use solve to handle a more complicated equation. Let n be an integer, we are interested in solving the equation x 3 + nx =. ( ) Do it yourself. (Theory) (F) n ( ) [0, ] S n x x 3 x nx 0 x. Prove that for every, Equation has a unique real solution in the interval. This solution will be denoted by. 2. On the same gure, plot, and for and for several (large) values of n. What can we conjecture for the limit S n? 3. (Bonus) Prove that S n /n, and deduce that S n = n + o(/n).
9 RangeOf_x=np.arange(0,,0.0) plt.plot(rangeof_x,rangeof_x**3,label='$x\mapsto x$**3') for n in [0, 20, 30]: f=[-n*x for x in RangeOf_x] plt.plot(rangeof_x,f,label='n ='+str(n)+' ') plt.xlim(0, ),plt.ylim(-, ) plt.xlabel('value of x') plt.legend() plt.title('location of S_n') plt.show() Answers.. The map is continuous and increasing on, since Besides, x f (x) = x 3 + nx [0, ] f (x) = 3 x 2 + n > n > 0. f (0) = 0 3 n 0 =, f () = 3 + n = n > 0. By the intermediate value theorem, this implies that there is a unique S n (0, ) such that f ( S n ) = 0, i.e. ( S n ) 3 + n S nn =. +, the solution of Equation 2. On the gure above we observe that when ( ) seems to get closer and closer to zero. We therefore conjecture lim S n n + 3. In fact, in Question ) one could be more precise: By the intermediate value theorem, this implies that S n (0, /n). If we use the fact that ( S n ) 3 + n S n =, we obtain since S n 0. Thus S n = n + o(/n). = 0. f (0) = 0 3 n 0 =, f (/n) = (/n) 3 + n /n = (/n) 3 > 0. n S n = ( S n ) 3 = + o(),
10 Do it yourself. (G). Write a script which computes the exact expression of S n. 2. In the previous questions it is proved that S n = n + o(/n). Use series to obtain the next orders in the asymptotic behaviour of S n (up to (/ n 5 ) ). var('x') var('n',integer=true) # Question. Solutions=solve(x**3+n*x-,x) Sn=Solutions[0] # The two other solutions are complex numbers print("sn = "+str(sn)) # Question 2. Taylor=series(Sn,n,oo,5) print("taylor expansion when epsilon -> 0 : "+str(taylor)) Sn = -n/(3*(sqrt(n**3/27 + /4) + /2)**(/3)) + (sqrt(n**3/27 + /4) + /2)**(/3) Taylor expansion when epsilon -> 0 : -/n**4 + /n + O(n**(-5), (n, oo) ) Answers. According to the above script, S n = + (/ n 5 ). n n 4 Exercise 5: Solving equations and "big O's" We consider the following equation: X 5 3εX 4 = 0, ( ) ε ε > 0 where is a positive parameter. A quick analysis shows that if is small enough then ( ) has a unique real solution, that we denote by S ε. The degree of this equation is too high to be solved by SymPy : var('x') var('e') solve(x**5-3*e*x**4-,x) [] Indeed, SymPy needs help to handle this equation. f (x) = x 5 3εx 4 ε In the above script we plotted the function for some small. This suggests that. lim ε 0 S ε =
11 RangeOf_x=np.arange(0,2,0.0) def ToBeZero(x,eps): return x**5+x**4*(-3*eps) - eps=0.05 plt.plot(rangeof_x,[tobezero(x,eps) for x in RangeOf_x],label='x**5+x**4*(-3*eps)-') plt.xlim(0, 2) plt.ylim(-, ) plt.plot([-2,2],[0,0]) plt.plot([,],[-2,2]) plt.xlabel('value of x') plt.title('location of S_n, with eps ='+str(eps)) plt.legend() plt.show() Do it yourself. (H). We admit that can be written as = + rε + s + ), S ε ε 2 ( ε 3 S ε r, s r, s for some real. Use SymPy to nd. (You can use any SymPy function already seen in this notebook.) var('r') var('s') var('eps') Expression=ToBeZero(+r*eps+s*eps**2+O(eps**3),eps) Simple=simplify(Expression) print(simple) solve([-3+5*r,5*s-2*r+0*r**2],[r,s]) -3*eps + 5*eps*r + 5*eps**2*s - 2*eps**2*r + 0*eps**2*r**2 + O(eps**3 ) [(3/5, 8/25)]
12 + rε + s ε 2 + ( ε 3 ) Answers. If we plug into equation ( ) we obtain (with the script): 0 = 3ε + 5rε + 5sε 2 2r ε r 2 ε 2 + O( ε 3 ). ( ) ε 0 = 3 + 5r + 5sε 2rε + 0r 2 ε + O( ε 2 ), 3 + 5r = 0 ε 0 r = 3/5 If we divide equation ( ) by we obtain which yields by letting, i.e.. ε 2 If we plug this into ( ) and divide by we obtain 0 = 5s 2r + 0 r 2 + O(ε), which gives, i.e.. 5s 2r + 0 r 2 = 0 s = 8/25 Exercise 6: Decomposition of fractions (Bonus). Do it yourself. (J) Set With SymPy nd reals (You can use solve, simplify, expand. Remember to use Rational(,) for fractions.) f (x) =. 2x 3 0x 2 30x + 54 a, b, c, α, β, γ such that a b c f (x) = + +. x α x β x γ a b c Answers. First, we observe that and must 2x 3 0x 2 30x+54 x α + + x β x γ have the same domain of de nition. Therefore {Roots of 2x 3 0x 2 30x + 54} = {α, β, γ}. In the script below we obtain: α = 3, β = 7 + 4, γ = var('x a b c') Num=2*x**3-0*x**2-30*x + 54 Roots=solve(Num,x) print(roots) [alpha, beta, gamma] = Roots [-3, -sqrt(7) + 4, sqrt(7) + 4]
13 Answers. We must have a b = x 3 0x 2 30x + 54 x + 3 x c x 4 7 x = 0,, 2 for. We obtain with the following script that x x (x + 3) # Left-hand side: def LHS(x): return Rational(,2*x**3-0*x**2-30*x + 54) # Right-hand side: def RHS(x,a,b,c): return a/(x-alpha)+b/(x-beta)+c/(x-gamma) # we identify a,b,c by solving the following system: Solutions=solve([LHS(0)-RHS(0,a,b,c),LHS()-RHS(,a,b,c),LHS(2)-RHS(2,a,b,c)],[a,b,c] print(solutions) astar=solutions[a] bstar=solutions[b] cstar=solutions[c] anum=n(astar) bnum=n(bstar) cnum=n(cstar) # to get a nice formula: var('x') print(latex(rhs(x,astar,bstar,cstar))) {c: -/68 + sqrt(7)/68, b: -sqrt(7)/68 - /68, a: /84} \frac{- \frac{\sqrt{7}}{68} - \frac{}{68}}{x \sqrt{7}} + \frac {- \frac{}{68} + \frac{\sqrt{7}}{68}}{x \sqrt{7}} + \frac{}{8 4 \left(x + 3\right)}
Symbolic computation 2: Linear recurrences
Bachelor of Ecole Polytechique Computatioal Mathematics, year 2, semester Lecturer: Lucas Geri (sed mail) (mailto:lucas.geri@polytechique.edu) Symbolic computatio 2: Liear recurreces Table of cotets Warm-up
More informationArithmetic 1: Prime numbers and factorization (with Solutions)
Bachelor of Ecole Polytechique Computatioal Mathematics, year 2, semester 1 Arithmetic 1: Prime umbers ad factorizatio (with Solutios) P. Tchebychev (Russia, 1821-1894). Amog may thigs, he was the rst
More information1. Pace yourself. Make sure you write something on every problem to get partial credit. 2. If you need more space, use the back of the previous page.
***THIS TIME I DECIDED TO WRITE A LOT OF EXTRA PROBLEMS TO GIVE MORE PRACTICE. The actual midterm will have about 6 problems. If you want to practice something with approximately the same length as the
More information1 The Real Number System
1 The Real Number System The rational numbers are beautiful, but are not big enough for various purposes, and the set R of real numbers was constructed in the late nineteenth century, as a kind of an envelope
More informationEXTRA CREDIT FOR MATH 39
EXTRA CREDIT FOR MATH 39 This is the second, theoretical, part of an extra credit homework. This homework in not compulsory. If you do it, you can get up to 6 points (3 points for each part) of extra credit
More informationWEEK 7 NOTES AND EXERCISES
WEEK 7 NOTES AND EXERCISES RATES OF CHANGE (STRAIGHT LINES) Rates of change are very important in mathematics. Take for example the speed of a car. It is a measure of how far the car travels over a certain
More informationA polynomial expression is the addition or subtraction of many algebraic terms with positive integer powers.
LEAVING CERT Honours Maths notes on Algebra. A polynomial expression is the addition or subtraction of many algebraic terms with positive integer powers. The degree is the highest power of x. 3x 2 + 2x
More informationchapter 1 vector geometry solutions V Consider the parallelogram shown alongside. Which of the following statements are true?
chapter vector geometry solutions V. Exercise A. For the shape shown, find a single vector which is equal to a)!!! " AB + BC AC b)! AD!!! " + DB AB c)! AC + CD AD d)! BC + CD!!! " + DA BA e) CD!!! " "
More informationMathematics 1 Lecture Notes Chapter 1 Algebra Review
Mathematics 1 Lecture Notes Chapter 1 Algebra Review c Trinity College 1 A note to the students from the lecturer: This course will be moving rather quickly, and it will be in your own best interests to
More informationTaylor Series. Math114. March 1, Department of Mathematics, University of Kentucky. Math114 Lecture 18 1/ 13
Taylor Series Math114 Department of Mathematics, University of Kentucky March 1, 2017 Math114 Lecture 18 1/ 13 Given a function, can we find a power series representation? Math114 Lecture 18 2/ 13 Given
More informationConditioning and Stability
Lab 17 Conditioning and Stability Lab Objective: Explore the condition of problems and the stability of algorithms. The condition number of a function measures how sensitive that function is to changes
More informationLab 1: Iterative Methods for Solving Linear Systems
Lab 1: Iterative Methods for Solving Linear Systems January 22, 2017 Introduction Many real world applications require the solution to very large and sparse linear systems where direct methods such as
More informationMathematics Revision Guide. Algebra. Grade C B
Mathematics Revision Guide Algebra Grade C B 1 y 5 x y 4 = y 9 Add powers a 3 a 4.. (1) y 10 y 7 = y 3 (y 5 ) 3 = y 15 Subtract powers Multiply powers x 4 x 9...(1) (q 3 ) 4...(1) Keep numbers without
More informationLesson 7: Algebraic Expressions The Commutative and Associative Properties
: Algebraic Expressions The Commutative and Associative Properties Four Properties of Arithmetic: The Commutative Property of Addition: If a and b are real numbers, then a + b = b + a. The Associative
More informationCS 237 Fall 2018, Homework 06 Solution
0/9/20 hw06.solution CS 237 Fall 20, Homework 06 Solution Due date: Thursday October th at :59 pm (0% off if up to 24 hours late) via Gradescope General Instructions Please complete this notebook by filling
More informationChapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27
Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation
More informationCS 237 Fall 2018, Homework 07 Solution
CS 237 Fall 2018, Homework 07 Solution Due date: Thursday November 1st at 11:59 pm (10% off if up to 24 hours late) via Gradescope General Instructions Please complete this notebook by filling in solutions
More informationLinear Transformations
Lab 4 Linear Transformations Lab Objective: Linear transformations are the most basic and essential operators in vector space theory. In this lab we visually explore how linear transformations alter points
More informationAdding and Subtracting Terms
Adding and Subtracting Terms 1.6 OBJECTIVES 1.6 1. Identify terms and like terms 2. Combine like terms 3. Add algebraic expressions 4. Subtract algebraic expressions To find the perimeter of (or the distance
More informationSTEP Support Programme. STEP 2 Complex Numbers: Solutions
STEP Support Programme STEP Complex Numbers: Solutions i Rewriting the given relationship gives arg = arg arg = α. We can then draw a picture as below: The loci is therefore a section of the circle between
More informationCHAPTER 1 POLYNOMIALS
1 CHAPTER 1 POLYNOMIALS 1.1 Removing Nested Symbols of Grouping Simplify. 1. 4x + 3( x ) + 4( x + 1). ( ) 3x + 4 5 x 3 + x 3. 3 5( y 4) + 6 y ( y + 3) 4. 3 n ( n + 5) 4 ( n + 8) 5. ( x + 5) x + 3( x 6)
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson
JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises
More informationLecture01_Geometric_view_of_linear_systems
Lecture_Geometric_view_of_linear_systems 5/3/8, : PM This notebook is part of lecture The geometry of linear equations in the OCW MIT course 8.6 by Prof Gilbert Strang [] Created by me, Dr Juan H Klopper
More informationMATH 1A, Complete Lecture Notes. Fedor Duzhin
MATH 1A, Complete Lecture Notes Fedor Duzhin 2007 Contents I Limit 6 1 Sets and Functions 7 1.1 Sets................................. 7 1.2 Functions.............................. 8 1.3 How to define a
More informationLecture # 1 - Introduction
Lecture # 1 - Introduction Mathematical vs. Nonmathematical Economics Mathematical Economics is an approach to economic analysis Purpose of any approach: derive a set of conclusions or theorems Di erences:
More informationSolutions, 2004 NCS/MAA TEAM COMPETITION
Solutions, 004 NCS/MAA TEAM COMPETITION Each problem number is followed by an 11-tuple (a 10, a 9, a 8, a 7, a 6, a 5, a 4, a 3, a, a 1, a 0 ), where a k is the number of teams that scored k points on
More informationSection 11.1: Sequences
Section 11.1: Sequences In this section, we shall study something of which is conceptually simple mathematically, but has far reaching results in so many different areas of mathematics - sequences. 1.
More informationAlgebra 1 Unit 6B Factoring
Algebra 1 Unit 6B Factoring Monday Tuesday Wednesday Thursday Friday 9 A Day 10 B Day 11 A Day 12 B Day 13 A Day Test Exponents and Polynomials Factor GCF and Trinomials box method Factoring Trinomials
More informationComplex Numbers. Visualize complex functions to estimate their zeros and poles.
Lab 1 Complex Numbers Lab Objective: Visualize complex functions to estimate their zeros and poles. Polar Representation of Complex Numbers Any complex number z = x + iy can be written in polar coordinates
More informationNotes from Andrew Reeve,U.Maine
Differentiation Summary We have now worked through the basic methods for differentiating an equation. y = x n dy dx = n xn 1 y = e x dy dx = ex y = sin(x) dy dx = cos(x) y = cos(x) dy dx = sin(x) y = ln(x)
More informationError, Accuracy and Convergence
Error, Accuracy and Convergence Error in Numerical Methods i Every result we compute in Numerical Methods is inaccurate. What is our model of that error? Approximate Result = True Value + Error. x = x
More informationFurther Mathematics Summer work booklet
Further Mathematics Summer work booklet Further Mathematics tasks 1 Skills You Should Have Below is the list of the skills you should be confident with before starting the A-Level Further Maths course:
More informationUsing Inductive and Deductive Reasoning
Big Idea 1 CHAPTER SUMMARY BIG IDEAS Using Inductive and Deductive Reasoning For Your Notebook When you make a conjecture based on a pattern, you use inductive reasoning. You use deductive reasoning to
More informationAnalysis/Calculus Review Day 3
Analysis/Calculus Review Day 3 Arvind Saibaba arvindks@stanford.edu Institute of Computational and Mathematical Engineering Stanford University September 15, 2010 Big- Oh and Little- Oh Notation We write
More informationComplex Numbers. A complex number z = x + iy can be written in polar coordinates as re i where
Lab 20 Complex Numbers Lab Objective: Create visualizations of complex functions. Visually estimate their zeros and poles, and gain intuition about their behavior in the complex plane. Representations
More informationAll numbered readings are from Beck and Geoghegan s The art of proof.
MATH 301. Assigned readings and homework All numbered readings are from Beck and Geoghegan s The art of proof. Reading Jan 30, Feb 1: Chapters 1.1 1.2 Feb 6, 8: Chapters 1.3 2.1 Feb 13, 15: Chapters 2.2
More information2. Limits at Infinity
2 Limits at Infinity To understand sequences and series fully, we will need to have a better understanding of its at infinity We begin with a few examples to motivate our discussion EXAMPLE 1 Find SOLUTION
More informationThe Geometric Mean and the AM-GM Inequality
The Geometric Mean and the AM-GM Inequality John Treuer February 27, 2017 1 Introduction: The arithmetic mean of n numbers, better known as the average of n numbers is an example of a mathematical concept
More information7.3 Singular points and the method of Frobenius
284 CHAPTER 7. POWER SERIES METHODS 7.3 Singular points and the method of Frobenius Note: or.5 lectures, 8.4 and 8.5 in [EP], 5.4 5.7 in [BD] While behaviour of ODEs at singular points is more complicated,
More informationMATH 18.01, FALL PROBLEM SET # 2
MATH 18.01, FALL 2012 - PROBLEM SET # 2 Professor: Jared Speck Due: by Thursday 4:00pm on 9-20-12 (in the boxes outside of Room 2-255 during the day; stick it under the door if the room is locked; write
More informationUpdated: January 16, 2016 Calculus II 7.4. Math 230. Calculus II. Brian Veitch Fall 2015 Northern Illinois University
Math 30 Calculus II Brian Veitch Fall 015 Northern Illinois University Integration of Rational Functions by Partial Fractions From algebra, we learned how to find common denominators so we can do something
More informationEdexcel New GCE A Level Maths workbook Circle.
Edexcel New GCE A Level Maths workbook Circle. Edited by: K V Kumaran kumarmaths.weebly.com 1 Finding the Midpoint of a Line To work out the midpoint of line we need to find the halfway point Midpoint
More informationSolutions to O Level Add Math paper
Solutions to O Level Add Math paper 04. Find the value of k for which the coefficient of x in the expansion of 6 kx x is 860. [] The question is looking for the x term in the expansion of kx and x 6 r
More informationMATH 408N PRACTICE FINAL
2/03/20 Bormashenko MATH 408N PRACTICE FINAL Show your work for all the problems. Good luck! () Let f(x) = ex e x. (a) [5 pts] State the domain and range of f(x). Name: TA session: Since e x is defined
More informationAS PURE MATHS REVISION NOTES
AS PURE MATHS REVISION NOTES 1 SURDS A root such as 3 that cannot be written exactly as a fraction is IRRATIONAL An expression that involves irrational roots is in SURD FORM e.g. 2 3 3 + 2 and 3-2 are
More informationCHAPTER 1 NUMBER SYSTEMS. 1.1 Introduction
N UMBER S YSTEMS NUMBER SYSTEMS CHAPTER. Introduction In your earlier classes, you have learnt about the number line and how to represent various types of numbers on it (see Fig..). Fig.. : The number
More informationWorksheet Week 1 Review of Chapter 5, from Definition of integral to Substitution method
Worksheet Week Review of Chapter 5, from Definition of integral to Substitution method This worksheet is for improvement of your mathematical writing skill. Writing using correct mathematical expressions
More informationfive line proofs Victor Ufnarovski, Frank Wikström 20 december 2016 Matematikcentrum, LTH
five line proofs Victor Ufnarovski, Frank Wikström 20 december 2016 Matematikcentrum, LTH projections The orthogonal projections of a body on two different planes are disks. Prove that they have the same
More informationPre RMO Exam Paper Solution:
Paper Solution:. How many positive integers less than 000 have the property that the sum of the digits of each such number is divisible by 7 and the number itself is divisible by 3? Sum of Digits Drivable
More information1 Question related to polynomials
07-08 MATH00J Lecture 6: Taylor Series Charles Li Warning: Skip the material involving the estimation of error term Reference: APEX Calculus This lecture introduced Taylor Polynomial and Taylor Series
More information0 Real Analysis - MATH20111
0 Real Analysis - MATH20111 Warmup questions Most of you have seen limits, series, continuous functions and differentiable functions in school and/or in calculus in some form You might wonder why we are
More information5615 Chapter 3. April 8, Set #3 RP Simulation 3 Working with rp1 in the Notebook Proper Numerical Integration 5
5615 Chapter 3 April 8, 2015 Contents Mixed RV and Moments 2 Simulation............................................... 2 Set #3 RP Simulation 3 Working with rp1 in the Notebook Proper.............................
More information2007 Paper 1 Solutions
27 Paper 1 Solutions 2x 2 x 19 x 2 + x + 2 1 = 2x2 x 19 (x 2 + x + 2) x 2 + x + 2 2x 2 x 19 x 2 + x + 2 > 1 2x 2 x 19 x 2 + x + 2 1 > x 2 4x 21 x 2 + x + 2 > (x + )(x 7) (x + 2)(x + 1) > = x2 4x 21 x 2
More informationUdaan School Of Mathematics Class X Chapter 10 Circles Maths
Exercise 10.1 1. Fill in the blanks (i) The common point of tangent and the circle is called point of contact. (ii) A circle may have two parallel tangents. (iii) A tangent to a circle intersects it in
More informationSample Question Paper Mathematics First Term (SA - I) Class IX. Time: 3 to 3 ½ hours
Sample Question Paper Mathematics First Term (SA - I) Class IX Time: 3 to 3 ½ hours M.M.:90 General Instructions (i) All questions are compulsory. (ii) The question paper consists of 34 questions divided
More informationAQA Level 2 Further mathematics Further algebra. Section 4: Proof and sequences
AQA Level 2 Further mathematics Further algebra Section 4: Proof and sequences Notes and Examples These notes contain subsections on Algebraic proof Sequences The limit of a sequence Algebraic proof Proof
More informationIntegrals. D. DeTurck. January 1, University of Pennsylvania. D. DeTurck Math A: Integrals 1 / 61
Integrals D. DeTurck University of Pennsylvania January 1, 2018 D. DeTurck Math 104 002 2018A: Integrals 1 / 61 Integrals Start with dx this means a little bit of x or a little change in x If we add up
More informationTest 2 - Answer Key Version A
MATH 8 Student s Printed Name: Instructor: CUID: Section: Fall 27 8., 8.2,. -.4 Instructions: You are not permitted to use a calculator on any portion of this test. You are not allowed to use any textbook,
More informationPROBLEM SET 3: PROOF TECHNIQUES
PROBLEM SET 3: PROOF TECHNIQUES CS 198-087: INTRODUCTION TO MATHEMATICAL THINKING UC BERKELEY EECS FALL 2018 This homework is due on Monday, September 24th, at 6:30PM, on Gradescope. As usual, this homework
More informationSpring Nikos Apostolakis
Spring 07 Nikos Apostolakis Review of fractions Rational expressions are fractions with numerator and denominator polynomials. We need to remember how we work with fractions (a.k.a. rational numbers) before
More informationChapter 3 - Derivatives
Chapter 3 - Derivatives Section 3.1: The derivative at a point The derivative of a function f (x) at a point x = a is equal to the rate of change in f (x) at that point or equivalently the slope of the
More informationRelevant sections from AMATH 351 Course Notes (Wainwright): Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls):
Lecture 5 Series solutions to DEs Relevant sections from AMATH 35 Course Notes (Wainwright):.4. Relevant sections from AMATH 35 Course Notes (Poulin and Ingalls): 2.-2.3 As mentioned earlier in this course,
More informationMATH 2433 Homework 1
MATH 433 Homework 1 1. The sequence (a i ) is defined recursively by a 1 = 4 a i+1 = 3a i find a closed formula for a i in terms of i.. In class we showed that the Fibonacci sequence (a i ) defined by
More informationADDITIONAL MATHEMATICS
005-CE A MATH HONG KONG CERTIFICATE OF EDUCATION EXAMINATION 005 ADDITIONAL MATHEMATICS :00 pm 5:0 pm (½ hours) This paper must be answered in English 1. Answer ALL questions in Section A and any FOUR
More informationb n x n + b n 1 x n b 1 x + b 0
Math Partial Fractions Stewart 7.4 Integrating basic rational functions. For a function f(x), we have examined several algebraic methods for finding its indefinite integral (antiderivative) F (x) = f(x)
More informationb) Rectangular box: length L, width W, height H, volume: V = LWH, cube of side s, V = s 3
Basic Math Review for PHYS 100 - Physics of Everyday Experience ----------------------------------------------------------------------------------------------------- Basic Algebra a) If x = y + z, then:
More informationWhy does pi keep popping up? Undergraduate Colloquium, October 2007 I. Definitions and Archimedes. II. Digits and some silliness (and Ramanujan)
Why does pi keep popping up? Undergraduate Colloquium, October 7 I. Definitions and Archimedes II. Digits and some silliness (and Ramanujan III. Antidote: pi is irrational. IV. Pi popping up in factorials.
More informationCLASS-IX MATHEMATICS. For. Pre-Foundation Course CAREER POINT
CLASS-IX MATHEMATICS For Pre-Foundation Course CAREER POINT CONTENTS S. No. CHAPTERS PAGE NO. 0. Number System... 0 3 0. Polynomials... 39 53 03. Co-ordinate Geometry... 54 04. Introduction to Euclid's
More information3 Fields, Elementary Matrices and Calculating Inverses
3 Fields, Elementary Matrices and Calculating Inverses 3. Fields So far we have worked with matrices whose entries are real numbers (and systems of equations whose coefficients and solutions are real numbers).
More informationMath 101 Study Session Spring 2016 Test 4 Chapter 10, Chapter 11 Chapter 12 Section 1, and Chapter 12 Section 2
Math 101 Study Session Spring 2016 Test 4 Chapter 10, Chapter 11 Chapter 12 Section 1, and Chapter 12 Section 2 April 11, 2016 Chapter 10 Section 1: Addition and Subtraction of Polynomials A monomial is
More informationHomework and Computer Problems for Math*2130 (W17).
Homework and Computer Problems for Math*2130 (W17). MARCUS R. GARVIE 1 December 21, 2016 1 Department of Mathematics & Statistics, University of Guelph NOTES: These questions are a bare minimum. You should
More informationHigher Mathematics Course Notes
Higher Mathematics Course Notes Equation of a Line (i) Collinearity: (ii) Gradient: If points are collinear then they lie on the same straight line. i.e. to show that A, B and C are collinear, show that
More informationMath Lecture 3 Notes
Math 1010 - Lecture 3 Notes Dylan Zwick Fall 2009 1 Operations with Real Numbers In our last lecture we covered some basic operations with real numbers like addition, subtraction and multiplication. This
More informationthat if a b (mod m) and c d (mod m), then ac bd (mod m) soyou aren't allowed to use this fact!) A5. (a) Show that a perfect square must leave a remain
PUTNAM PROBLEM SOLVING SEMINAR WEEK 2 The Rules. You are not allowed to try a problem that you already know how to solve. These are way too many problems to consider. Just pick a few problems in one of
More informationDISCOVERING GEOMETRY Over 6000 years ago, geometry consisted primarily of practical rules for measuring land and for
Name Period GEOMETRY Chapter One BASICS OF GEOMETRY Geometry, like much of mathematics and science, developed when people began recognizing and describing patterns. In this course, you will study many
More informationAP Calculus Chapter 9: Infinite Series
AP Calculus Chapter 9: Infinite Series 9. Sequences a, a 2, a 3, a 4, a 5,... Sequence: A function whose domain is the set of positive integers n = 2 3 4 a n = a a 2 a 3 a 4 terms of the sequence Begin
More informationUnit 3 Factors & Products
1 Unit 3 Factors & Products General Outcome: Develop algebraic reasoning and number sense. Specific Outcomes: 3.1 Demonstrate an understanding of factors of whole number by determining the: o prime factors
More informationDecember 2012 Maths HL Holiday Pack. Paper 1.2 Paper 1 from TZ2 Paper 2.2 Paper 2 from TZ2. Paper 1.1 Paper 1 from TZ1 Paper 2.
December 2012 Maths HL Holiday Pack This pack contains 4 past papers from May 2011 in the following order: Paper 1.2 Paper 1 from TZ2 Paper 2.2 Paper 2 from TZ2 Paper 1.1 Paper 1 from TZ1 Paper 2.1 Paper
More informationLimits at. x means that x gets larger and larger without a bound. Oktay Olmez and Serhan Varma Calculus Lecture 2 1 / 1
Limits at x means that x gets larger and larger without a bound. Oktay Olmez and Serhan Varma Calculus Lecture 2 1 / 1 Limits at x means that x gets larger and larger without a bound. x means that x gets
More informationp and q are two different primes greater than 25. Pass on the least possible value of p + q.
A1 p and q are two different primes greater than 25. Pass on the least possible value of p + q. A3 A circle has an area of Tπ. Pass on the area of the largest square which can be drawn inside the circle.
More informationMath Precalculus I University of Hawai i at Mānoa Spring
Math 135 - Precalculus I University of Hawai i at Mānoa Spring - 2014 Created for Math 135, Spring 2008 by Lukasz Grabarek and Michael Joyce Send comments and corrections to lukasz@math.hawaii.edu Contents
More informationMath Precalculus I University of Hawai i at Mānoa Spring
Math 135 - Precalculus I University of Hawai i at Mānoa Spring - 2013 Created for Math 135, Spring 2008 by Lukasz Grabarek and Michael Joyce Send comments and corrections to lukasz@math.hawaii.edu Contents
More informationMath 3320 Foundations of Mathematics
Math 3320 Foundations of Mathematics Chapter 1: Fundamentals Jesse Crawford Department of Mathematics Tarleton State University (Tarleton State University) Chapter 1 1 / 55 Outline 1 Section 1.1: Why Study
More informationChapter 5 Simplifying Formulas and Solving Equations
Chapter 5 Simplifying Formulas and Solving Equations Look at the geometry formula for Perimeter of a rectangle P = L W L W. Can this formula be written in a simpler way? If it is true, that we can simplify
More informationBHP BILLITON UNIVERSITY OF MELBOURNE SCHOOL MATHEMATICS COMPETITION, 2003: INTERMEDIATE DIVISION
BHP BILLITON UNIVERSITY OF MELBOURNE SCHOOL MATHEMATICS COMPETITION, 00: INTERMEDIATE DIVISION 1. A fraction processing machine takes a fraction f and produces a new fraction 1 f. If a fraction f = p is
More informationPUTNAM TRAINING EASY PUTNAM PROBLEMS
PUTNAM TRAINING EASY PUTNAM PROBLEMS (Last updated: September 24, 2018) Remark. This is a list of exercises on Easy Putnam Problems Miguel A. Lerma Exercises 1. 2017-A1. Let S be the smallest set of positive
More informationLeast squares and Eigenvalues
Lab 1 Least squares and Eigenvalues Lab Objective: Use least squares to fit curves to data and use QR decomposition to find eigenvalues. Least Squares A linear system Ax = b is overdetermined if it has
More informationNotater: INF3331. Veronika Heimsbakk December 4, Introduction 3
Notater: INF3331 Veronika Heimsbakk veronahe@student.matnat.uio.no December 4, 2013 Contents 1 Introduction 3 2 Bash 3 2.1 Variables.............................. 3 2.2 Loops...............................
More informationIntroduction Derivation General formula Example 1 List of series Convergence Applications Test SERIES 4 INU0114/514 (MATHS 1)
MACLAURIN SERIES SERIES 4 INU0114/514 (MATHS 1) Dr Adrian Jannetta MIMA CMath FRAS Maclaurin Series 1/ 19 Adrian Jannetta Background In this presentation you will be introduced to the concept of a power
More informationNOTES ON OCT 4, 2012
NOTES ON OCT 4, 0 MIKE ZABROCKI I had planned a midterm on Oct I can t be there that day I am canceling my office hours that day and I will be available on Tuesday Oct 9 from 4-pm instead I am tempted
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 informationMath Review for Exam Answer each of the following questions as either True or False. Circle the correct answer.
Math 22 - Review for Exam 3. Answer each of the following questions as either True or False. Circle the correct answer. (a) True/False: If a n > 0 and a n 0, the series a n converges. Soln: False: Let
More informationMATH 1A - FINAL EXAM DELUXE - SOLUTIONS. x x x x x 2. = lim = 1 =0. 2) Then ln(y) = x 2 ln(x) 3) ln(x)
MATH A - FINAL EXAM DELUXE - SOLUTIONS PEYAM RYAN TABRIZIAN. ( points, 5 points each) Find the following limits (a) lim x x2 + x ( ) x lim x2 + x x2 + x 2 + + x x x x2 + + x x 2 + x 2 x x2 + + x x x2 +
More informationAlgebraic Expressions
Algebraic Expressions 1. Expressions are formed from variables and constants. 2. Terms are added to form expressions. Terms themselves are formed as product of factors. 3. Expressions that contain exactly
More informationMATH 104A HW 02 SOLUTION KEY
.3 Algorithms and Convergence Solutions by Jon Lo Kim Lin - Fall 04 MATH 04A HW 0 SOLUTION KEY You are encouraged to collaborate with your classmates and utilize internet resources provided you adhere
More informationStepping stones for Number systems. 1) Concept of a number line : Marking using sticks on the floor. (1 stick length = 1 unit)
Quality for Equality Stepping stones for Number systems 1) Concept of a number line : Marking using sticks on the floor. (1 stick length = 1 unit) 2) Counting numbers: 1,2,3,... Natural numbers Represent
More informationTaylor and Maclaurin Series
Taylor and Maclaurin Series MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Background We have seen that some power series converge. When they do, we can think of them as
More informationf ( c ) = lim{x->c} (f(x)-f(c))/(x-c) = lim{x->c} (1/x - 1/c)/(x-c) = lim {x->c} ( (c - x)/( c x)) / (x-c) = lim {x->c} -1/( c x) = - 1 / x 2
There are 9 problems, most with multiple parts. The Derivative #1. Define f: R\{0} R by [f(x) = 1/x] Use the definition of derivative (page 1 of Differentiation notes, or Def. 4.1.1, Lebl) to find, the
More informationNote: An algebraic expression is made of one or more terms. The terms in an algebraic expression are connected with either addition or subtraction.
.4 Combining Like Terms Term A term is a single number or variable, or it can be a product of a number (called its coefficient) and one or more variables. Examples of terms:, x, a, xy, 4 a bc, 5, xyz coefficient
More informationGrade 8 Chapter 7: Rational and Irrational Numbers
Grade 8 Chapter 7: Rational and Irrational Numbers In this chapter we first review the real line model for numbers, as discussed in Chapter 2 of seventh grade, by recalling how the integers and then the
More information