BÔ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC THƯƠNG MẠI KHOA HỆ THỐNG THÔNG TIN KINH TẾ & THƯƠNG MẠI ĐIỆN TỬBẢN BÁO CÁO TỔNG HỢP HỌC PHẦN TOÁN CAO CẤP 2 Nhiệm vụ: “Sử dụng Python để giải các bà
Trang 1BÔ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC THƯƠNG MẠI KHOA HỆ THỐNG THÔNG TIN KINH TẾ & THƯƠNG MẠI ĐIỆN TỬ
BẢN BÁO CÁO TỔNG HỢP
HỌC PHẦN TOÁN CAO CẤP 2
Nhiệm vụ: “Sử dụng Python để giải các bài toán”.
Lớp học phần: 2192FMAT0211 Nhóm thực hiện: 03
Giảng viên giảng dạy: Lê Văn Tuấn
HÀ NỘI – 2021
Trang 2DANH SÁCH THÀNH VIÊN
Nhóm trưởng: Lê Thị Giang
ST
21 Đinh Thùy Dương Làm nhiệm vụ 1 và kiểm tra các đáp án
23 Lê Thị Giang Làm nhiệm vụ 6 và tổng hợp đáp án
24 Nguyễn Thị Hương Giang Làm nhiệm vụ 2
25 Nguyễn Thu Hà Làm nhiệm vụ 3, word
28 Nguyễn Thị Thúy Hằng Làm nhiệm vụ 4
29 Nguyễn Thúy Hằng Làm nhiệm vụ 1
30 Trần Thị Thu Hiền Làm nhiệm vụ 5
Trang 3NHIỆM VỤ CỦA NHÓMNV1 Tính đạo hàm riêng cấp 1 & 2 của các hàm số trong bài 7.2 tại (3, 12).
NV2 Tính các tích phân trong bài 8.1 với cận 3 đến 12.
NV3 Vẽ đồ thị nghiệm riêng của các PTVP trong bài 9.4 với điều kiện ban đầu y(3) = 2
Trang 4from sympy import *
x, y = symbols ('x y', real=True)
import numpy as np
from sympy import *
x, y = symbols ('x y', real=True)
from sympy import *
x, y = symbols ('x y', real=True)
Trang 5x, y = symbols('x y', real=True)
from sympy import *
x, y = symbols('x y', real=True)
from sympy import *
x, y = symbols ('x y', real=True)
import numpy as np
from sympy import *
x, y = symbols ('x y', real=True)
Trang 6import numpy as np
from sympy import *
x, y = symbols ('x y', real=True)
from sympy import *
x, y = symbols('x y', real=True)
from sympy import *
x, y = symbols('x y', real=True)
from sympy import *
x, y = symbols ('x y', real=True)
Trang 7z_x = lambdify ((x, y), z_x)
print (z_x(3,12)
= - 0.295314730609362
import numpy as np
from sympy import *
x, y = symbols ('x y', real=True)
from sympy import *
x, y = symbols ('x y', real=True)
from sympy import *
x, y = symbols('x y', real=True)
from sympy import *
x, y = symbols('x y', real=True)
Trang 8from sympy import *
x, y = symbols ('x y', real=True)
from sympy import *
x, y = symbols ('x y', real=True)
from sympy import *
x, y = symbols ('x y', real=True)
Trang 9z_xx = lambdify((x, y), z_xx)
print(z_xx(3,12))
import numpy as np
from sympy import *
x, y = symbols('x y', real=True)
from sympy import *
x, y = symbols('x y', real=True)
from sympy import *
x, y = symbols ('x y', real=True)
Trang 10from sympy import *
x, y = symbols ('x y', real=True)
from sympy import *
x, y = symbols ('x y', real=True)
from sympy import *
x, y = symbols('x y', real=True)
from sympy import *
x, y = symbols('x y', real=True)
Trang 11from sympy import *
x, y = symbols ('x y', real=True)
from sympy import *
x, y = symbols ('x y', real=True)
from sympy import *
x, y = symbols ('x y', real=True)
from sympy import *
x, y = symbols('x y', real=True)
Trang 12import numpy as np
from sympy import *
x, y = symbols('x y', real=True)
from sympy import *
x, y = symbols ('x y', real=True)
15 ¿.15sin3 sin 12
import numpy as np
from sympy import *
x, y = symbols ('x y', real=True)
15 ¿ 15sin 3.sin 12
import numpy as np
from sympy import *
x, y = symbols ('x y', real=True)
Trang 13from sympy import *
x, y = symbols('x y', real=True)
from sympy import *
x, y = symbols('x y', real=True)
Trang 14from scipy import integrate
f = lambda x: math.sqrt(3*x+5) / x
result = integrate.quad(f, 3, 12)
print(result)
6.7758332587582695,1.1351502655320377e-10
Giá trị tích phân là:
0.0817030676982187
Trang 170 0
Trang 18NHIỆM VỤ 3: Vẽ đồ thị nghiệm riêng của các PTVP trong bài 9.4 với điều kiện ban đầu y(3) = 2 trên miền [3, 12].
y3 = 2
x = np.linspace(3,12)
y = odeint(model,y3,x)plt.plot(x,y)
plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
2 y ' + =4 x y
import math
Trang 19x = np.linspace(3,12)
y = odeint(model,y3,x)plt.plot(x,y)
plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
Trang 20return dydxy3 = 2
x = np.linspace(3,12)
y = odeint(model,y3,x)plt.plot(x,y)
plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
Trang 21plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
ValueError: math domain error - lỗi vì tập xác định của arcsin(x) là x € (-1;1)
Trang 22y3 = 2
x = np.linspace(3,12)
y = odeint(model,y3,x)plt.plot(x,y)
plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
Trang 23def model(y,x):
dydx = (2*x**4 + 2*y)/x return dydx
y3 = 2
x = np.linspace(3,12)
y = odeint(model,y3,x)plt.plot(x,y)
plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
Trang 24y = odeint(model,y3,x)plt.plot(x,y)
plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
ValueError: math domain error - lỗi vì tập xác định của arcsin(x) là x € (-1;1)
plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
Trang 250 0
Trang 26NHIỆM VỤ 4: Vẽ đồ thị nghiệm riêng của các PTVP trong bài 9.7 với điều kiện y(3)
= 2, y’(3) = 1 trên miền [3, 12].
Trang 350 0
Trang 36NHIỆM VỤ 5: Vẽ đồ thị nghiệm riêng của các PTSP trong bài 10.4 với điều kiện ban đầu x(0) = 3 trên miền [0, 9].
for n in range(1, N):
x[n] = n*x[n-1]+math.factorial(n)plt.plot(x)
plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
2 x(n+1) + nx(n) = n !
Trang 37for n in range(1, N):
x[n] = math.factorial(n-1)-(n-1)*x[n-1]plt.plot(x)
plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
Trang 38N = 9
x = np.zeros(N, float)x[0] = 3
for n in range(1, N):
x[n] = math.pow(9,n-1)*x[n-1]+math.pow(3,(n-1)**2+2*(n-1))plt.plot(x)
plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
4 x(n+1) + n3 x(n) = n!
3nimport math
import numpy as np
import matplotlib.pyplot as plt
N = 9
x = np.zeros(N, int)x[0] = 3
Trang 39for n in range(1, N):
x[n] = math.factorial(n-1)/math.pow(3,n-1)-(n-1)/3*x[n-1]plt.plot(x)
plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
5 x(n+1) = n2 x(n) + n !
2n+1import math
import numpy as np
import matplotlib.pyplot as plt
N = 9
x = np.zeros(N, int)x[0] = 3
for n in range(1, N):
x[n] =(n-1)/2*x[n-1]+math.factorial(n-1)/math.pow(2,n)plt.plot(x)
Trang 40plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
for n in range(1, N):
x[n] = math.pow(4,-(n-1))*x[n-1]+(n-1)*math.pow(2,-(n-1)**2-(n-1))plt.plot(x)
plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()
Trang 420 0
Trang 430 0
Trang 44NHIỆM VỤ 4: Vẽ đồ thị nghiệm riêng của các PTVP trong bài 9.7 với điều kiện y(3)
= 2, y’(3) = 1 trên miền [3, 12].
Trang 527 y ' ' − y =2 sin −4 cos x x
i t th