68
Để thoát ra khỏi vòng lặp for hoặc while, ta sử dụng lệnh break. Khi chương trình
chạy đến dòng lệnh này thì nó tự động nhảy ra khỏi vòng lặp có chứa break. Ví dụ
for p = 7:8
for q = 3:5
for r = 1:2
fprintf('\n %3.0f, %3.0f, %3.0fn' , p, q, r)
end
if q == 4, break; end
end
end
fprintf('\n Xuong dong \n')
% Ket qua la
7, 3, 1n
7, 3, 2n
7, 4, 1n
7, 4, 2n
8, 3, 1n
8, 3, 2n
8, 4, 1n
8, 4, 2n
Xuong dong
Ví dụ sau đây chỉ ra sự tương tác giữa các hàm khi lập trình. Xét phương trình vi
phân
( , )
y
f t y
=
ɺ
với điều kiện đầu
0
0
( )
y t
y
=
.
Ta cần tìm nghiệm
( )
y t
trong khoảng thời gian
0
[ ,
]
f
t
t
t
=
.
Để giải bài toán bằng phương pháp Euler, ta sẽ chia khoảng thời gian trên thành
N
đoạn, với độ dài hay còn gọi là bước thời gian là
0
(
)/
f
h
t
t
N
=
−
. Gọi
( )
i
i
z
z t
=
là trị gần đúng của hàm
( )
y t
tại thời điểm
0
, (
1, 2,...,
)
i
t
t
ih
i
N
=
+
=
, theo
phương pháp Euler dãy các giá trị
i
z
được tính theo công thức sau
0
0
1
1
1
( ),
(
,
)
i
i
i
i
z
y t
z
z
h f t
z
−
−
−
=
=
+
⋅
Khi tiến hành thực hiện trong Matlab ta sẽ soạn ba m-file: một mô tả phương trình
vi phân, một mô tả phương pháp gần đúng Euler, và một chương trình chính.
Cụ thể để giải phương trình vi phân
1
0
0
(
), ( )
0
m
y
F
by y t
=
−
=
ɺ
,
0
0,
=10
f
t
t
=
.
Ta biết gnhiệm chính xác của phương trình này là
0
( /
)
( )
[1
]
F
b m t
b
y t
e
−
=
−
