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

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à

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à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

DANH SÁCH THÀNH VIÊN

Nhóm trưởng: Lê Thị Giang

ST

NHIỆM VỤ CỦA NHÓM

NV1 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

from 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)

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)

import numpy as np

from 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)

from sympy import *

x, y = symbols('x y', real=True)

from sympy import *

x, y = symbols ('x y', real=True)

z_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)

z = exp(atan(y/x))

Màn hình hiển thị:

- 0.010133348599340853

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)

z_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)

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)

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)

import 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.sin12

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 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

NHIỆ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 ' + y=4 x

import math

x = np.linspace(3,12)

y = odeint(model,y3,x)plt.plot(x,y)

plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()

return 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()

plt.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)

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()

def 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()

y = 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()

NHIỆ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].

import math

7 y ' ' − y=2 sin x−4 cos x

10 y ' ' − y=2 e x

−x2import math

NHIỆ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 !

for 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()

N = 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) + 3n x(n) = n!

3nimport 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] = 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)

plt.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()

for n in range(1, N):

x[n] = (n-1)*x[n-1]+math.factorial(n-1)*math.log(n)plt.plot(x)

plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()

for n in range(1, N):

math.exp(2*(n-1))*x[n-1]+math.exp((n-1)**2+(n-1))*math.log(n)plt.plot(x)

plt.xlabel('Truc x')plt.ylabel('Truc y')plt.show()

NHIỆM VỤ 6: Vẽ đồ thị nghiệm riêng của các PTSP trong bài 10.7 với điều kiện: x(0) =

plt.xlabel('Truc n')

plt.ylabel('Truc x')

plt.show()

plt.xlabel('Truc n')

plt.ylabel('Truc x')

plt.show()

plt.xlabel('Truc n')

plt.ylabel('Truc x')

plt.show()

BẢNG ĐÁNH GIÁ ĐIỂM 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

8

23 Lê Thị Giang Làm nhiệm vụ 6 và tổng hợp đáp án 8