BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC THƯƠNG MẠI KHOA KINH TẾ - LUẬT... Bùi Phương Anh3.. Lê Thị Nguyệt Anh 8.. Nguyễn Thị Ngọc Anh 9.. Nguyễn Thị Vân Anh II.. Đánh giá: Các thành viên h
Trang 1BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC THƯƠNG MẠI
KHOA KINH TẾ - LUẬT
Trang 22 Bùi Phương Anh
3 Chu Hải Anh
4 Đặng Vân Anh
5 Hoàng Đức Anh
6 Lê Hoàng Quỳnh Anh
7 Lê Thị Nguyệt Anh
8 Nguyễn Thị Ngọc Anh
9 Nguyễn Thị Vân Anh
II Nội dung cuộc họp:
1 Họp nhóm lần 1
Thời gian: Ngày 21/04/2020 từ 19h đến 20h
Địa điểm: phần mềm Zoom
Số thành viên: 9/9
Nội dung:
- Nhóm trưởng thông báo nội dung đề tài thảo luận
- Lên kế hoạch thảo luận và phân công nhiệm vụ cho từng thành viên.Phân công nhiệm vụ:
1) NV1: Ngô Phúc An
2) NV2: Nguyễn Thị Vân Anh
3) NV3: Chu Hải Anh
4) NV4: Đặng Vân Anh, Lê Hoàng Quỳnh Anh
5) NV5: Nguyễn Thị Ngọc Anh
Trang 36) NV6: Hoàng Đức Anh, Bùi Phương Anh
7) Tổng hợp: Lê Thị Nguyệt Anh
2 Họp nhóm lần 2
Thời gian: ngày 10/5/2020 từ 19h đến 20h
Địa điểm: phầm mềm Zoom
Số thành viên tham gia: 9/9
Nội dung cuộc họp:
- Nộp nhiệm vụ thảo luận
- Nhóm trưởng đánh giá phần làm của mọi người
III Đánh giá:
Các thành viên hoàn thành nhiệm vụ, nhiệt tình, hăng hái
Hà Nội, ngày 10 tháng 5 năm 2020 Nhóm trưởng
Trang 4NV2: Tính các tích phân trong bài 8.1 với cận 1 đến 10 6
NV3: Vẽ phương trình nghiệm riêng của PTVP trong bài 9.4 với điều kiện ban đầu y(1)=2 trên miền [1;10] 8
NV4: Vẽ đồ thị nghiệm riêng của các PTVP trong bài 9.7 với điều kiện y(1) = 2, y’(1) = 1 trên miền [1,10] 18
NV5: Vẽ đồ thị nghiệm riêng của PTSP trong bài 10.4 với điều kiện ban đầu y(1) =
2 trên miền [1,10] 28
NV6: Vẽ đồ thị nghiệm riêng của các PTSP trong bài 10.7 với đk: y(1) = 2, y(10) =
5 trên miền [1,10] 36
Trang 5NV1: Tính đạo hàm riêng cấp 1 và 2 các hàm số trong bài 7.2 tại (1,10)
Trang 6> deriv(expression(exp(atan(y/x))),c("x","y"),function.arg = TRUE,hessian = TRUE)(1,10)
Trang 10Z’’xx(1,10) =0.5355604; Z’’xy(1,10) =0.0402399; Z’’yx(1,10) =0.0402399; Z”yy(1,10) = 1.328594
NV2: Tính các tích phân trong bài 8.1 với cận 1 đến 10
Error in integrate(f, lower = 1, upper = 10) :
the integral is probably divergent
4.
Giải :
6
Trang 12Error in integrate(f, lower = 1, upper = 10) : non-finite function value
In addition: Warning message:
In log(tan(x)) : NaNs produced
NV3: Vẽ phương trình nghiệm riêng của PTVP trong bài 9.4 với điều kiện ban đầu y(1)=2 trên miền [1;10]
1.
Giải:
8
Trang 13> ham <- function (x, y,parms) {list(c((2*x/(1+x^2))*y))}
> dieukien <- c(y = 2)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0)
> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
2.
Giải:
9
Trang 14> ham <- function (x, y,parms) {list(c(4*x-y))}
> dieukien <- c(y = 2)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0)
> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
3.
10
Trang 15> ham <- function (x, y,parms) {list(c(x*exp(-x^2)-2*x*y))}
> dieukien <- c(y = 2)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0)
> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
11
Trang 16Giải:
> ham <- function (x, y,parms) {list(c(sin(x)*cos(x)+y*sin(x)))}
> dieukien <- c(y = 2)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0)
> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
12
Trang 175
Giải:
> ham <- function (x, y,parms) {list(c(asin(x)+x-(x*y)/(1-x^2)))}
> dieukien <- c(y = 2)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0)
Chú thích: Trong phương trình, điều kiện của arcsinx xác định trên [-1;1], mà điều
kiện ban đầu xét phương trình chạy trên [1;10]
Không vẽ được phương trình vi phân trên [1;10]
13
Trang 18Giải:
> ham <- function (x, y,parms) {list(c(2*x*(x^2+y)))}
> dieukien <- c(y = 2)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0)
> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
14
Trang 20> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
8.
Giải:
> ham <- function (x, y,parms) {list(c((asin(x)-y)/(1-x^2)^1/2))}
> dieukien <- c(y = 2,y_next=0)
16
Trang 210 0
Trang 22> ham <- function (x, y,parms) {list(c((2*x/(1+x^2))*y))}
> dieukien <- c(y = 2)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0)
> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
2.
Giải:
9
Trang 23> ham <- function (x, y,parms) {list(c(4*x-y))}
> dieukien <- c(y = 2)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0)
> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
Trang 2410
:
Trang 25> ham <- function (x, y,parms) {list(c(x*exp(-x^2)-2*x*y))}
> dieukien <- c(y = 2)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0)
> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
11
Trang 26Giải:
Trang 27> ham <- function (x, y,parms) {list(c(sin(x)*cos(x)+y*sin(x)))}
> dieukien <- c(y = 2)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0)
> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
12
Trang 285
Giải:
Trang 29> ham <- function (x, y,parms) {list(c(asin(x)+x-(x*y)/(1-x^2)))}
> dieukien <- c(y = 2)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0)
Chú thích: Trong phương trình, điều kiện của arcsinx xác định trên [-1;1], mà điều
kiện ban đầu xét phương trình chạy trên [1;10]
Không vẽ được phương trình vi phân trên [1;10]
13
Trang 30Giải:
> ham <- function (x, y,parms) {list(c(2*x*(x^2+y)))}
> dieukien <- c(y = 2)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0)
> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
14
Trang 31Giải:
Trang 32> ham <- function (x, y,parms) {list(c((2*x^4+2*y)/x))}
> dieukien <- c(y = 2)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0)
15
Trang 33> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
8.
Giải:
Trang 34> ham <- function (x, y,parms) {list(c((asin(x)-y)/(1-x^2)^1/2))}
> dieukien <- c(y = 2,y_next=0)
16
Trang 35> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0)
Chú thích: Trong phương trình, điều kiện arcsinx xác định trên [-1;1], mà điều kiện
ban đầu xét phương trình chạy trên [1;10]
Không vẽ được phương trình vi phân trên [1;10]
9.
Giải:
Trang 36> ham <- function (x, y,parms) {list(c(x*log(x)+y/(x*log(x))))}
> dieukien < c(y = 2 y next=exp(2)/2)
Trang 37> dieukien <- c(y = 2,y_next=exp(2)/2)
> nghiem <- ode(y = dieukien, func = ham,times = 2:10, parms = 0)
> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
Chú thích: Xét phương trình ta có:
Tại x=1 thì xlnx = 0, phương trình không xác định sinh ra lỗi
x=1 không thỏa mãn
Ta xét phương trình chạy trên [2;10]
NV4: Vẽ đồ thị nghiệm riêng của các PTVP trong bài 9.7 với điều kiện y(1) =
2, y’(1) = 1 trên miền [1,10]
1.
> ham <- function (x, y,parms) {list(c(y[2], (2*y[2]+y[1])))}
> dieukien <- c(y = 2, y_phay = 1)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 1)
> plot(nghiem, type = "l", which = "y",lwd = 2, ylab = "Truc y",main = "PTVP cấp 2")
18
Trang 380 5.0e+08 1.0e+09 1.5e+09
Trang 39> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 1)
> plot(nghiem, type = "l", which = "y",lwd = 2, ylab = "Truc y",main = "PTVP cấp 2")
19
Trang 41> ham <- function (n, x,parms) {list(c(x[2], 3*x[2]-2*x[1]+2*n+3))}
> dieukien <- c(y = 2, y_next = 5)
37
Trang 42> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0, method =
Trang 43> ham <- function (n, x,parms) {list(c(x[2], 3*x[2]-2*x[1]+4*n*3^n))}
> dieukien <- c(y = 2, y_next = 5)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0, method =
"iteration")
> plot(nghiem, type = "l", which = "y",lwd = 2, ylab = "Truc y",main = "PTSP")
39
Trang 44> ham <- function (n, x,parms) {list(c(x[2], 4*x[2]-4*x[1]+(5^n)*(9*n+3)))}
> dieukien <- c(y = 2, y_next = 5)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0, method =
"iteration")
> plot(nghiem, type = "l", which = "y",lwd = 2, ylab = "Truc y",main = "PTSP")
40
Trang 45> ham <- function (n, x,parms) {list(c(x[2], 4*x[2]-4*x[1]+16*2^n+18*5^n))}
> dieukien <- c(y = 2, y_next = 5)
41
Trang 46> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0, method =
Trang 47> ham <- function (n, x,parms) {list(c(x[2],
5*x[2]-6*x[1]+cos((2*n+pi)/2)+12*sin((2*n+pi)/2)))}
> dieukien <- c(y = 2, y_next = 5)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0, method =
"iteration")
> plot(nghiem, type = "l", which = "y",lwd = 2, ylab = "Truc y",main = "PTSP")
43
Trang 48> ham <- function (n, x,parms) {list(c(x[2], 2*x[2]-4*x[1]+14*3^n))}
> dieukien <- c(y = 2, y_next = 5)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0, method =
"iteration")
> plot(nghiem, type = "l", which = "y",lwd = 2, ylab = "Truc y",main = "PTSP")
44
Trang 49> ham <- function (n, x,parms) {list(c(x[2],
x[1]+6*cos((2*n+pi)/2)-4*sin((2*n+pi)/2)))}
> dieukien <- c(y = 2, y_next = 5)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0, method =
"iteration")
> plot(nghiem, type = "l", which = "y",lwd = 2, ylab = "Truc y",main = "PTSP")
45
Trang 50> ham <- function (n, x,parms) {list(c(x[2], 3*x[2]-2*x[1]+2*n+1+2^(n+2)))}
> dieukien <- c(y = 2, y_next = 5)
> nghiem <- ode(y = dieukien, func = ham,times = 1:10, parms = 0, method =
"iteration")
> plot(nghiem, type = "l", which = "y",lwd = 2, ylab = "Truc y",main = "PTSP")
46
Trang 5147