Đá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ưởngNguyệt Anh Lê Thị Nguyệt Anh MỤC LỤC... TIEU LUAN MOI download : skknchat
Trang 1BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC THƯƠNG MẠI
KHOA KINH TẾ - LUẬT
-BÀI THẢO LUẬN
MÔN: TOÁN CAO CẤP 2
Giảng viên hướng dẫn: Lê Văn Tuấn
Mã lớp học phần: 2040FMAT0211 Sinh viên thực hiện: Nhóm 1
Hà Nội, tháng 5 năm 2020 CỘNG HOÀ XÃ HỘI CHỦ
NGHĨA VIỆT NAM Độc lập - Tự do – Hạnh phúc
Trang 2II 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:
Phân công nhiệm vụ:
Trang 36) NV6: Hoàng Đức Anh, Bùi Phương 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:
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ưởngNguyệt Anh
Lê Thị Nguyệt Anh
MỤC LỤC
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 đầuy(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)
2 z =)
Giải :
>deriv(expression(log(x+sqrt(x^2+y^2))),c("x","y"),function.arg = TRUE,hessian
= TRUE)(1,10)
Trang 8[1,] NaN NaN
Trang 10TIEU LUAN MOI download : skknchat123@gmail.com
Trang 11Error in integrate(f, lower = 1, upper = 10) :
the integral is probably divergent
4.
Trang 12> f<-function(x){sin(4*x)/(cos(2*x))^2+4}
Error in integrate(f, lower = 1, upper = 10) :
the integral is probably divergent
Trang 13Error 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:
Trang 14> ham <- function (x, y,parms) {list(c((2*x/(1+x^2))*y))}
> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
2.
Giải:
Trang 15> ham <- function (x, y,parms) {list(c(4*x-y))}
> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
Trang 16> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
Trang 17Giải:
>plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
Trang 18Giải:
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]
Trang 19Giải:
>plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
Trang 21>plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
8.
Giải:
Trang 23> 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 2510
Trang 26> 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")
Trang 27Giải:
Trang 28> 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 29Giải:
Trang 30> 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 31Giải:
> ham <- function (x, y,parms) {list(c(2*x*(x^2+y)))}
Trang 33> 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 34> plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")
8.
Giải:
Trang 35> ham <- function (x, y,parms) {list(c((asin(x)-y)/(1-x^2)^1/2))}
> dieukien <- c(y = 2,y_next=0)
16
Trang 36> 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 3717
Trang 38> 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 41> 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")
19
Trang 42TIEU LUAN MOI download : skknchat123@gmail.com
Trang 43PTVP cấp 2
Trang 442.
Trang 45>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")
Trang 46TIEU LUAN MOI download : skknchat123@gmail.com
Trang 47> ham <- function (n, x,parms) {list(c(x[2], 3*x[2]-2*x[1]+4*n*3^n))}
"iteration")
> plot(nghiem, type = "l", which = "y",lwd = 2, ylab = "Truc y",main =
"PTSP")
Trang 48TIEU LUAN MOI download : skknchat123@gmail.com
Trang 49"iteration")
> plot(nghiem, type = "l", which = "y",lwd = 2, ylab = "Truc y",main =
"PTSP")
Trang 50TIEU LUAN MOI download : skknchat123@gmail.com
Trang 515.
Trang 52>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")
6.
Trang 53TIEU LUAN MOI download : skknchat123@gmail.com
Trang 54>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)))}
"iteration")
> plot(nghiem, type = "l", which = "y",lwd = 2, ylab = "Truc y",main =
"PTSP")
Trang 55TIEU LUAN MOI download : skknchat123@gmail.com
Trang 56"iteration")
> plot(nghiem, type = "l", which = "y",lwd = 2, ylab = "Truc y",main =
"PTSP")
Trang 57TIEU LUAN MOI download : skknchat123@gmail.com
Trang 58x[1]+6*cos((2*n+pi)/2)-4*sin((2*n+pi)/2)))}
"iteration")
> plot(nghiem, type = "l", which = "y",lwd = 2, ylab = "Truc y",main =
"PTSP")
Trang 59TIEU LUAN MOI download : skknchat123@gmail.com
Trang 60"iteration")
Trang 61TIEU LUAN MOI download : skknchat123@gmail.com