Scribd
Upload
7 views

Python 2

Uploaded by

pomenmew
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
7 views

Python 2

Uploaded by

pomenmew
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 5

Ton Duc Thang University

Faculty of Information
Technology

Applied Calculus for IT - 501031


Lab 04

1. Sympy
SymPy is a Python library for symbolic mathematics. It aims to become a full-
featured computer algebra system (CAS) while keeping the code as simple as
possible in order to be comprehensible and easily extensible.
1.1 Setup:
If you have Python and PIP already installed on a system, then install it using
this command:
python -m pip install sympy
1.2 Introduction
Symbolic computation systems such as SymPy are capable of computing
symbolic expressions with variables. Let us define a symbolic expression,
representing the mathematical expression x+2y.

Note that we wrote x +2∗y just as we would if x and y were ordinary Python
variables. But in this case, instead of evaluating to something, the expression
remains as just x +2∗y . Now let us play around with it:

https://docs.sympy.org/latest/tutorials/intro-tutorial/intro.html
1.3 Practice examples:
Solve x 2+ 4 x + 4=0.
import sympy as sp
Trinh Hung Cuong - Applied Calculus for IT - 501031 1/5
Ton Duc Thang University
Faculty of Information
Technology

x = sp.symbols('x')
root = sp.solve(x**2 + 4*x + 4)
print(root)

Plot functions:
import sympy as sp
import math

x = sp.symbols('x')

f5a = abs(x)**(1/2)
sp.plot(f5a, (x, -10, 10), line_color='red') #tuple

f5d = math.e**x
sp.plot(f5d, (x, -10, 10), line_color='blue')

f5e = sp.log(x)
sp.plot(f5e, (x, -10, 10))

Evaluate a certain expression:


import sympy as sp

x = sp.symbols('x')
f = x*x - x + 1

print( f.subs(x, 2) ) # f(2)

Exercise 0
Write a Python program to plot the following functions on a graph, and mark the intersection point of
f1 and f2:

Hint:
import sympy as sp
import matplotlib.pyplot as plt
import numpy as np

Trinh Hung Cuong - Applied Calculus for IT - 501031 2/5


Ton Duc Thang University
Faculty of Information
Technology
# ve f1 va f2
x = np.arange(...)
f1 = lambda x: ...
f2 = lambda x: ...

y1 = ...
y2 = ...

plt.plot(x, y1)
plt.plot(x, y2)

# ve diem giao nhau


x = sp.symbols('x')
f1 = -x + 5
f2 = ...

x_root = sp.solve(...)
y_root = f1.subs(...) #f1(x_root[0]) = f1(2)

plt.plot(..., y_root, 'ro')


plt.title('Find intersection point of f1(x) = -x + 5 and f2(x) = x/2 + 2')
plt.grid(linestyle='--')
plt.show()

Exercises 1, 3, 4, 5, 6, 7, 10 in the PDF file “Lab04.ex.pdf”


Hint:
Limit of a function:
import sympy as sp
import math
x = sp.symbols('x')
#1c
f1c = math.e ** (1/x)
lm = sp.limit(f1c, x, 1)
print("1c - The limit of f(x) at x = 1: " + str(lm) )
#1i
f = (2*x*x) / ( 3 - 3*sp.cos(x) )
lm = sp.limit(f, x, 0)
print("1i - The limit of f(x) at x = 0: " + str(lm) )

#infinity: sp.oo

#factorial: sp.factorial(...)

Trinh Hung Cuong - Applied Calculus for IT - 501031 3/5


Ton Duc Thang University
Faculty of Information
Technology

Right/Left limit:

import sympy as sp
import math

#3.1
x = sp.symbols('x')
f3_1 = 1/( 1 + 2**(1/x) )
lmRight = sp.limit(f3_1, x, 0, '+')
print ("Right limit = " , lmRight )
lmLeft = sp.limit(f3_1, x, 0, '-')
print ("Left limit = " , lmLeft )

Continuity

Trinh Hung Cuong - Applied Calculus for IT - 501031 4/5


Ton Duc Thang University
Faculty of Information
Technology

Ex6:
import numpy as np
import sympy as sp

x = sp.symbols('x')
f6a = ...

#At point x = 0
lm_x_0 = sp.limit(f6a, x, 0)
#Compare lm_x_0 and f(0)
...

#Other points x != 0
for c in np.arange(-100, 100, 1):
if c != 0:
lm_x_c = sp.limit(f6a, x, c)
#Compare lm_x_c and f(c)

#f(x) is continuous for all x # 0, 3

Homework:
Exercises 8 in the PDF file “Lab04.ex.pdf”

2. References
- Python Tutorial on the W3schools website:
https://www.w3schools.com/python/default.asp
- Python Tutorial on the Tutorials Point website:
https://www.tutorialspoint.com/python/index.htm

-- THE END --

Trinh Hung Cuong - Applied Calculus for IT - 501031 5/5

You might also like