Gợi ý giải Đề thi kết thúc học phần môn PHƯƠNG PHÁP SỐ VÀ LẬP TRÌNH - Học kỳ II - NH: 2013-2014

Câu 1: Ta có lưu đồ sau:

Câu 2:

Đáp án:

• t=0.2
  wk=(0.2, 0.6, 1.3, 2.15)
  sum=2.15
  
• t=0.3
  wk=(0.3, 0.65, 1.325, 2.1625)
  sum=2.1625
  
• t=0.4
  wk=(0.4, 0.7, 1.35, 2.175)
  sum=2.175
  
• t=0.5
  wk=(0.5, 0.75, 1.375, 2.1875)
  sum=2.1875

Câu 3:

Chương trình sai:

1     PROGRAM Noi suy Lagrange

2     IMPLICIT NONE

3     INTEGER::j,n

4     REAL::i,P,xnhap

5     PRINT*,'Moi ban nhap bac cua da thuc noi suy ='

6     READ*,n

7     REAL, ALLOCATABLE::x(:),y(:)

8     Do i=0,n

9         PRINT*,'moi ban nhap x_',i,'='

10         READ*,x(i)

11         PRINT*,'moi ban nhap y_',i,'='

12         READ*,y(i)

13    END DO

14    PRINT*,'Moi ban nhap diem x ma ban muon xac dinh gia tri ham'

15    READ*,xnhap

16    P=0

17    DO i=0,n

18        L=1.

19             DO j=0,n

20                 IF i/=j THEN

21             L=L*(xnhap-x(j))/(x(i)-x(j))

22                 END IF

23             END DO

24             P=P+L*y(i)

25    END DO

26    PRINT*, 'y(',xnhap,')=', P

27    END PROGRAM

Lỗi 1: Dòng 1: "PROGRAM Noi suy Lagrange", Tên chương trình không được chứa dấu cách.

  Sửa: Program NoisuyLagrange

Lỗi 2: Dòng 4: "REAL::i,P,xnhap", Biến i là biến đếm, phải là kiểu nguyên chứ không phải là kiểu thực.

  Sửa: Đưa biến i lên khai báo ở dòng 3 chung với các biến kiểu nguyên: "INTEGER::j,n,i"

Lỗi 3: Dòng 7: "REAL, ALLOCATABLE::x(:),y(:)", Phần khai báo biến phải nằm trong phần khai báo chứ không thể nằm ở phần thân chương trình.

  Sửa: Đưa dòng 7 lên ngay sau dòng 4.

Lỗi 4: Biến x, y là biến mảng động, phải được gán giá trị cho số phần tử của mảng.

  Sửa: Thêm vào ngay sau dòng 6 của chương trình dòng lệnh sau: Allocate(x(0:n),y(0:n))

Lỗi 5: Biến L chưa được khai báo biến.

  Sửa: Khai báo thêm vào dòng 4: REAL::P,xnhap, L

Lỗi 6: Dòng 20:"IF i/=j THEN", trong cú pháp câu lệnh IF thì điều kiện phải được đặt ở trong dấu ().

  Sửa: IF (i/=j) THEN

Sau đây là chương trình đã chỉnh sửa hoàn chỉnh:

1     PROGRAM NoisuyLagrange

2     IMPLICIT NONE

3     INTEGER::j,n,i

4     REAL::P,xnhap,L

5     REAL, ALLOCATABLE::x(:),y(:)

6     PRINT*,'Moi ban nhap bac cua da thuc noi suy ='

7     READ*,n

8     ALLOCATE(x(0:n),y(0:n))

9     Do i=0,n

10         PRINT*,'moi ban nhap x_',i,'='

11         READ*,x(i)

12         PRINT*,'moi ban nhap y_',i,'='

13         READ*,y(i)

14    END DO

15    PRINT*,'Moi ban nhap diem x ma ban muon xac dinh gia tri ham'

16    READ*,xnhap

17    P=0

18    DO i=0,n

19        L=1.

20             DO j=0,n

21                 IF (i/=j) THEN

22             L=L*(xnhap-x(j))/(x(i)-x(j))

23                 END IF

24             END DO

25             P=P+L*y(i)

26    END DO

27    PRINT*, 'y(',xnhap,')=', P

28    END PROGRAM

Câu 4:

Vậy phương trình có nghiệm y:

• Tại x=0.1 là y=0.18124
• Tại x=0.2 là y=0.26869

Bài làm trên chỉ mang tính chất tham khảo, nếu có gì sai sót mong được các bạn góp ý sửa chữa.

Read more

Đề thi kết thúc học phần môn PHƯƠNG PHÁP SỐ VÀ LẬP TRÌNH - Học kỳ II - NH: 2013-2014 - Thời gian: 90 phút

Câu 1: ( 3đ )

Vẽ lưu đồ cho chương trình EULSUM.F95 dưới đây:

PROGRAM eulsum

INTEGER::j,n

REAL::sum,t,dum,tmp

REAL,ALLOCATABLE::wk(:)

WRITE(*,*)'Input n'

READ(*,*)n

WRITE(*,*)'Input t'

READ(*,*)t

ALLOCATE(wk(n+1))

WRITE(*,*)'Input matrix wk'

READ(*,*)wk

sum=0

IF (n==1) THEN

  wk(1)=t

  sum=t/2

 ELSE

   tmp=wk(1)

   wk(1)=t

   DO j=1,n-1

     dum=wk(j+1)

     wk(j+1)=(wk(j)+tmp)/2

     tmp=dum

   END DO

   wk(n+1)=(wk(n)+tmp)/2

   IF (abs(wk(n+1))<=abs(wk(n))) THEN

     sum=sum+0.5*wk(n+1)

     ELSE

       sum=sum+wk(n+1)

   END IF

END IF

WRITE(*,*)wk,sum

END

Câu 2: ( 2đ )

Với chương trình ở câu 1, nếu các giá trị nhập vào ban đầu là: n=3; wk=(1, 2, 3, 4); t nhận một trong các giá trị sau: 0,2; 0,3; 0,4; 0,5 (SV chọn một trong các giá trị này để tính). Hãy cho biết giá trị của biến wk và sum sau khi kết thúc chương trình.

Câu 3: ( 3đ )

Hãy tìm 6 lỗi sai và đề xuất cách sửa cho chương trình dưới đây (không cần chép lại cả đề, chỉ viết lại dòng bị sai, gạch dưới chỗ bị lỗi và đề xuất cách sửa)

1     PROGRAM Noi suy Lagrange

2     IMPLICIT NONE

3     INTEGER::j,n

4     REAL::i,P,xnhap

5     PRINT*,'Moi ban nhap bac cua da thuc noi suy ='

6     READ*,n

7     REAL, ALLOCATABLE::x(:),y(:)

8     Do i=0,n

9         PRINT*,'moi ban nhap x_',i,'='

10         READ*,x(i)

11         PRINT*,'moi ban nhap y_',i,'='

12         READ*,y(i)

13    END DO

14    PRINT*,'Moi ban nhap diem x ma ban muon xac dinh gia tri ham'

15    READ*,xnhap

16    P=0

17    DO i=0,n

18        L=1.

19             DO j=0,n

20                 IF i/=j THEN

21             L=L*(xnhap-x(j))/(x(i)-x(j))

22                 END IF

23             END DO

24             P=P+L*y(i)

25    END DO

26    PRINT*, 'y(',xnhap,')=', P

27    END PROGRAM

Câu 4: ( 2đ )

Tìm nghiệm tại điểm x=0.1 và x=0.2 cho bài toán Cauchy dưới đây bằng phương pháp Euler cải tiến với bước nhảy là h = 0.1:
 \begin{cases} y'=e^{-2x}-2y+4x, 0\leq x\leq 2\\ y(0)=0.1\end{cases}

........................................Hết........................................

Read more

Phương pháp số và lập trình - Tính sai số

Read more

Phương pháp số và lập trình - Nội suy, Đạo hàm, Tích phân

Read more

Fortran cơ sở

Read more

Mở đâu môn Phương pháp số và lập trình

Read more

Phương pháp số và lập trình - Giải phương trình phi tuyến

Read more

Phương pháp số và lập trình - Giải phương trình vi phân

Read more

Đào vàng

Nhấn nút xuống để chơi game

Read more

Nông trại 2

Read more