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.